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PREFACE 


The  objective  of  this  basic  research  program  is  to 

investigate  the  feasibility  of  new  storage  techniques  for 

large  archival  memories  using  ion  and  electron  beam  with 

semiconductor  targets.  The  goal  is  the  development  of  an 
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archival  memory  capable  of  storing  10  to  10  bits  with  rapid 
access  to  the  stored  information. 

The  currently  funded  effort  (Phase  I)  deals  with 
feasibility  studies  of  key  technical  areas.  The  program  covers 
experimental  and  analytical  investigations:  (1)  to  demonstrate 
the  feasibility  of  ion  implanted  and  alloy  junction  storage 
media,  (2)  to  determine  the  feasibility  of  beam  optics  design 
required  for  writing/reading  on  the  media,  and  (3)  to  select  the 
better  of  the  two  storage  methods  and  perform  a preliminary  paper 
design  of  the  concept  for  further  development. 

During  this  reporting  period,  effort  in  the  archival  memory 
program  included  experimental  work  and  theoretical  studies  on  a 
planar  diode  structure  (memory  target  substrate),  implantation 
of  ions  to  write  information  on  a planar  diode  structure,, 
feasibility  studies  of  the  alloy  junction  method  of  writing,  and 
the  specification  of  a high  performance  electron  beam  write  station 
for  alloy  junction  writing. 


A detailed  analytical  modeling  study  was  initiated  in 
order  to  properly  design  the  planar  diode  and  to  understand 
the  results  of  experimental  testing.  The  initial  model  design 
was  used  during  the  period  to  compare  test  results  with  theory 
and  define  new  processing  experiments. 

An  improved  method  of  ion  writing  using  inert  ions  to  form 
damage  sites  was  shown  to  be  feasible.  Preliminary  results 
indicate  a lower  writing  fluence  will  be  required. 

Alloy  junction  diodes  were  formed  using  a laser  beam. 

Diodes  as  small  as  5 microns  were  formed. 

An  electron  optical  system,  consisting  of  a high  performance 
field  emission  source  and  focusing  optics,  was  specified  and  a 
purchase  order  was  placed. 
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Section  I 
SUMMARY 
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Program  Objectives 

The  objective  of  this  basic  research  program  is  to  in- 
vestigate the  feasibility  of  new  storage  techniques  for  large 
archival  memories  using  ion  and  electron  beams  with  semi- 
conductor targets.  The  goal  is  the  development  of  an  archival 
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memory  capable  of  Storing  10  to  10  bits  with  rapid  access 
to  the  stored  infromation. 

The  currently  funded  effort  (Phase  I)  deals  with  feasi- 
bility studies  of  key  technical  areas.  The  program  covers 
experimental  and  analytical  investigations:  (1)  to  demonstrate 
the  feasibility  of  ion  implanted  and  alloy  junction  storage  media, 
(2)  to  determine  the  feasibility  of  beam  optics  design  required 
for  writing/reading  on  the  media,  and  (3)  to  select  the  better 
of  the  two  storage  methods  and  perform  a preliminary  paper  design 
of  the  concept  for  further  development. 

Accomplishments 

During  this  reporting  period,  effort  in  the  archival 
memory  program  included  experimental  work  and  theoretical 
studies  on  a planar  diode  structure  (memory  target  substrate), 
implantation  of  ions  to  write  information  on  a planar  diode 
structure,  feasibility  studies  of  the  alloy  junction  method 
of  writing,  and  the  specification  of  a high  performance  electron 
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beam  write  station  for  alloy  junction  writing. 

1.  Planar  Diode  Studies 

The  target  substrate  on  which  data  will  be  stored  is  a 
large  area  planar  diode.  The  diode  provides  a mechanism  for 
signal  amplification  within  the  target  structure.  The  top 
layer  of  the  diode  is  much  thinner  than  more  normal  planar 
diodes  and  constitutes  a departure  in  design.  In  order  to 
resolve  the  0.1  micron  stored  data  bits,  the  junction  must  be 
within  ~ 1000^  of  the  surface  and  the  electron  beam  energy 
must  be  in  the  2-3  keV  range.  This  dictates  a diode  with  little 
or  no  "dead"  region  near  the  surface  and  a junction  depth  ~ 
1000A.  A "dead  layer"  is  defined  as  the  region  in  which 
carriers  generated  by  the  electron  beam  are  not  collected  by 
the  diode. 

Various  combinations  of  materials  are  possible.  Both 
germanium  and  silicon  with  a wide  variety  of  dopants  are 
candidate  diode  structures.  In  order  to  properly  design  the 
diode  and  to  understand  the  results  of  experimental  testing, 
a detailed  analytical  modeling  study  was  initiated.  The  model 
is  designed  to  consider  the  probability  of  carrier  capture  at 
the  planar  junction,  as  affected  by  doping  gradients  and  other 
parameters.  The  initial  model  design  was  used  during  the  period 
to  compare  test  results  with  theory  and  define  new  processing 
experiments. 
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2.  Ion  Writing  Experiments 

One  of  the  two  storage  concepts  under  study  is  to  write 
bits  of  information  as  small  (0.1  micron  diameter)  surface 
diodes  by  ion  implantation.  During  the  period  an  interesting 
possibility  for  an  improved  writing  mechanism  was  discovered. 
This  new  method  involves  the  use  of  inert  ions  to  form  damage 
sites  rather  than  doped  sites.  To  test  the  concept,  helium 
implants  were  made  and  showed  excellent  readout  characteristics. 
More  detailed  work  is  planned  for  the  next  reporting  period. 
However,  preliminary  results  would  indicate  that  a much  lower 
fluence  will  be  required  for  this  form  of  writing. 

3.  Alloy  Junction  Writing  Experiments 

Alloyed  junctions  between  a thin  metal  coating  and  the 
surface  of  a planar  diode  will  be  made  with  a high  brightness, 
finely  focused  electron  beam.  Until  an  experimental  electron 
beam  test  facility  is  established,  experiments  were  conducted 
using  a laser  as  the  writing  beam.  Alloy  junction  diodes  were 
formed  at  the  surface  of  a planar  diode  and  demonstrated  basic 
feasibility  for  this  type  of  archival  memory  storage.  Diodes 
in  size  of  5-7  microns  were  formed  (limited  by  the  spot  size 
of  the  focused  laser  beam).  Additional  experiments  are 
anticipated.  However,  future  efforts  will  be  more  concentrated 
on  the  thermal  mechanisms  at  the  metal -semiconductor  interface 
and  the  selection  of  the  most  desirable  systems. 
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4.  Electron  Beam  Write  Column 

An  electron  optical  system,  consisting  of  a high 
performance  field  emission  source  and  focusing  optics,  vas 
specified  and  a purchase  order  was  placed.  Delivery  is 
expected  during  the  second  quarter  of  the  contract.  This 
equipment  will  be  used  to  conduct  alloy  junction  writing 
experiments . 

Plans  for  Next  Quarter 

The  overall  Program  Schedule  lists  the  major  areas  of 
development  envisioned  for  the  next  period.  The  current 
investigation  will  continue  with  refinements  to  the  analytical 
modeling  studies.  In  addition,  studies  of  both  electron  and 
ion  optics  will  be  initiated  as  beam  performance  parameters 
become  known  as  a result  of  the  theoretical  and  experimental 
work.  The  electron  optical  system  for  alloy  junction  writing 
is  expected  during  the  quarter. 

Data  Included  in  this  Report 

To  provide  a fuller  understanding  of  the  physics  of  these 

i 

recording  methods,  a large  body  of  substantiating  information 
is  included  in  this  report.  Specifically  the  two  Appendices, 
"Ion  Implantation"  and  "Electron  Penetration  into  Solids" 
contain  background  information  pertinent  to  the  methods  under 


study. 
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Section  II 


MEMORY  CONCEPT 

The  approach  selected  for  the  development  of  a large 

advanced  archival  memory  with  rapid  access  employs  beams  for 

the  writing  and  reading  processes.  Finely  focused  beams  can 

permanently  inscribe  a material,  producing  a record  of  archival 

character.  Beams  can  also  be  used  to  read  or  sense  the  stored 

information.  Data  can  be  written,  read,  and  accessed  rapidly, 

because  beams  can  be  deflected  at  extremely  high  rates. 

14  15 

Finally,  for  a memory  of  10  - 10  bits,  the  density  of  storage 

is  a critical  parameter . Beams  can  be  imaged  to  spot  diameters 
that  are  fractions  of  a micron  in  size. 

This  research  effort  is  directed  toward  the  writing  of 
data  by  the  formation  of  small  diodes  at  the  surface  of  a semi- 
conductor plane  (silicon  or  germanium)  by  the  action  of  a finely 
focused  beam  (ion  or  electron).  The  presence  or  absence  of 
diodes  represents  stored  zeros  or  ones. 

For  writing  data,  two  methods  for  forming  surface  diodes 
will  be  evaluated,  and  the  more  promising  one  will  be  selected 
for  further  development.  The  two  writing  methods  are: 

• Ion  Implanted  Surface  Diodes 

• Alloy  Junction  Surface  Diodes 


— 
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Data  Storage 

Ion  implantation  of  surface  diodes  is  a candidate  storage 
technique.  A finely  focused  ion  beam  is  scanned  over  a small 
area  of  a silicon  wafer  and  modulated  on-off  to  form  by 
implantation  small  surface  diodes.  Once  the  field  is  written, 
the  wafer  is  stepped  to  a new  location  for  additional  writing. 

This  method  is  shown  schematically  in  Figure  1. 

It  is  the  goal  of  this  study  to  show  that  it  is  feasible 
to  form  0.1  micron  size  diodes  in  the  surface  of  a silicon 
wafer  by  use  of  a finely  focused  ion  beam.  It  is  believed  that 
diodes  of  this  size  can  be  formed  and  that  this  high  resolution 
can  be  preserved  because  the  implanted  ions  will  have  very 
little  lateral  spread.  Diodes  as  small  as  0.5  microns  have 
been  formed  by  ion  implantation  through  a mask. 

Formation  of  surface  diodes  by  the  alloy  junction  method 
involves  the  use  of  an  electron  beam  to  alloy  two  materials. 

As  example,  a germanium  target  coated  with  a thin  indium  layer 
would  be  exposed  to  a finely  focused  electron  beam.  A thermal 
bias  would  maintain  the  target  within  a few  degrees  of  the 
alloying  temperature  of  160°C.  The  electron  beam  would  locally 
heat  the  surface  to  160°C,  causing  the  alloy  to  form  where  data 
are  to  be  recorded.  Figure  2 illustrates  this  method  schematically. 
The  finely  focused  electron  beam  would  be  electronically  positioned 
to  any  location  on  the  memory  storage  plane  by  use  of  a matrix 
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Figure  1.  Schematic  Diagram  - Ion  Implanted  Surface  Diodes  Write-Read 
Stations . 
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Figure  2.  Schematic  Design  - Alloy  Junction  Surface  Diode 
Write-Read  Station. 


lens.  This  lens,  sometimes  called  a fly's  eye  lens,  is  an 
array  of  small  electron  lenses  (lenselets),  each  with  the 
capability  of  beam  deflection.  In  this  case  the  complete 
memory  plane  could  be  written  without  the  need  for  mechanical 
motion. 

This  technique  has  very  high  resolution  because  the 
thermal  diffusion  of  heat  in  a germanium  substrate  is  very 
rapid.  The  heated  area  would  be  confined  to  the  size  of  the 
electron  spot  itself. 

Electron  Beam  Read-Out  of  Stored  Data 

The  surface  diode  memory  plane  is  shown  schematically 
in  Figure  3 (a  and  b) . The  unwritten  target  consists  of  a 
thin  layer  of  heavily  doped  N-type  silicon  (N+  layer),  as 
example,  on  a lightly  doped  P-type  substrate.  Stored  in- 
formation consists  of  small  heavily  doped  P-type  regions  (P+ 
regions)  in  the  top  surface  of  the  N+  layer.  These  diodes 
can  be  written  by  either  of  the  methods  previously  described. 

To  read  the  stored  information,  the  area  is  scanned 
with  a well  focused  electron  beam,  with  energy  —2.5  keV.  This 
produces  one  electron  hole  pair  for  every  3.6  eV  of  beam  energy 
(—700  pairs  for  a 2.5  keV  beam).  When  the  electron  beam 
impinges  on  an  unrecorded  area  (one),  the  minority  carriers 
(holes)  diffuse  to  the  N+P  diode  depletion  region  and  generate 
a signal  current,  Ig,  in  the  external  diode  biasing  circuit 
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Figure  3 -a.  Readout  of  One 


Figure  3-b.  Readout  of  Zero 
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(Figure  3a).  This  diffusion  of  holes  is  aided  by  the  N+ 
layer  if  there  is  a doping  gradient  across  the  N+  layer 

■f 

from  the  top  surface  to  the  N P diode.  The  signal  from 

a one  will  typically  be  a large  fraction  of  the  maximum 

possible,  which  for  a 2.5  keV  electron  beams  is  700.  When 

the  electron  beam  impinges  on  a recorded  zero,  the  electron 

hole  pairs  are  trapped  by  the  local  diode  depletion  field. 

This  results  in  a much  smaller  readout  signal  (Figure  3b). 

The  physics  of  the  carrier  motion  for  zero  readout  is 

essentially  the  same  as  for  the  one  readout,  only  now  the 

+ + + 

collecting  junction  (P  N instead  of  N P)  is  not  connected 
to  the  readout  circuit  and  the  generated  carriers  recombine 
locally.  In  effect  these  electron  hole  pairs  cause  a slight 
collapse  of  the  zero  bias  field  across  the  small  P+N+  junction, 
which  quickly  relaxes  because  these  diodes  are  likely  to  be  of 
very  poor  quality  (i.e.,  leaky). 

The  readout  signal  from  such  a surface  diode  target 
would  then  be  ~500  for  unrecorded  areas  (ones)  and  nearly  zero 
for  recorded  areas  (zeros). 


12 


I 

I 

Section  III 
PLANAR  DIODE  STUDIES 


Int roduction 

During  the  first  quarter  of  the  Archival  Memory  Program  much 
of  the  experimental  work  has  centered  upon  achieving  the  appropriate 
planar  diode  substrate  on  which  the  information  may  be  later  written. 

Both  schemes  of  writing,  the  ion  implantation  and  alloy  junction 
methods,  require  such  a planar  diode  structure.  Ion  implantation 
was  employed  for  fabrication  of  these  planar  diodes  because  of  the 
uniformity  and  control  that  can  be  achieved  by  this  method. 

This  section  first  summarizes  the  surface  diode  target  concept 
with  emphasis  on  those  aspects  which  relate  to  the  requirements  on 
the  planar  diode  structure.  Then  work  during  the  first  quarter 
directed  towards  meeting  these  goals  is  discussed. 

Target  Concept 

The  target  concept  is  shown  schematically  in  Figure  4.  The 
unwritten  target  consists  of  a thin  layer  of  heavily  doped  N-type 
silicon  (N+  layer)  on  a lightly  doped  P-type  substrate.  Stored 
information  consists  of  small  heavily  doped  P-type  regions  (P+ 
regions)  or  damaged  regions  in  the  top  surface  of  the  N+  layer. 

The  complementary  structure  with  N+  diodes  in  a P+  surface  on  N is 

also  possible.  For  clarity,  on  the  N+  on  P structure  is  discussed. 

+ + 

The  small  P regions  form  small  unbiased  diodes  with  the  N layer. 

These  diodes  can  be  written  by  several  methods,  such  as  focused 

ion  beams,  alloy  writing,  or  variations  on  electron  beam  lithography. 

I 

i 


-13- 


To  read  the  stored  information,  the  area  is  scanned  with  a well 

focused  electron  beam,  with  energy  ~2.5  keV.  This  produces  one 

electron  hole  pair  for  every  3.6  eV  of  beam  energy  (~700  pairs  for 

a 2.5-keV  beam).  When  the  electron  beam  impinges  on  an  unrecorded 

area  (one),  the  minority  carriers  (holes)  diffuse  to  the  N+P  diode 

depletion  region  and  generate  a signal  current,  I , in  the  external 

s 

diode  biasing  circuit  (Figure  5).  This  diffusion  of  holes  is  aided 
by  the  N+  layer  if  there  is  a doping  gradient  across  the  N+  layer 
from  the  top  surface  to  the  N+P  diode.  The  signal  from  a one  will 
typically  be  a large  fraction  of  the  maximum  possible,  which  for 
2.5-keV  electron  beams  is  700.  When  the  electron  beam  impinges  on 
a recorded  zero,  the  electron  hole  pairs  are  trapped  by  the  local 
diode  depletion  field.  This  results  in  a much  smaller  readout 
signal  (Figure  6).  The  physics  of  the  carrier  motion  for  zero 
readout  is  essentially  the  same  as  for  the  one  readout,  only  now 
the  collecting  junction  (P+N+  instead  of  N+P)  is  not  connected  to 
the  readout  circuit  and  the  generated  carriers  recombine  locally. 

In  effect  these  electron  hole  pairs  cause  a slight  collapse  of  the 

• 4*  ^ 

zero  bias  field  across  the  small  P N junction,  which  quickly 

relaxes  because  these  diodes  are  likely  to  be  of  very  poor  quality 

(i.e.,  leaky).  To  further  short  the  diodes  and  speed  the  relaxation, 

a thin  metal  top  layer  or  very  heavy  N+  surface  doping  can  be  used. 


PIP 


The  readout  signal  from  such  a surface  diode  target  would  then 
be  ~500  for  unrecorded  areas  (ones)  and  nearly  zero  for  recorded 
areas  (zeros). 

Frequency  Response 

The  frequency  response  of  the  target  can  be  extremely  high. 

The  intrinsic  target  limitations  on  frequency  response  are  the 
diffusion  or  drift  of  carriers  to  and  across  the  depletion  layer 
of  the  N+P  junction  and  the  RC  response  of  the  diode  capacitance. 

2 

The  diffusion  time  can  be  estimated  from  the  relation  - L / D, 

where  D is  the  diffusion  constant  for  the  minority  carriers  (D  = 5 

2 *f" 

cm  /sec  for  holes  in  0.1  ohm-cm  n silicon),  and  L is  the  N layer 

-11 

thickness  of  ~0.1  micron.  This  yields  2 x 10  seconds,  which  is 
much  too  short  to  be  of  any  practical  importance. 

The  drift  across  the  diode  depletion  field  can  also  be  shown 

to  be  negligible  for  this  application.  The  transit  time  is 

o 

rt  = 2e/qN^|i,  where  }i  is  the  minority  carrier  mobility  = 480  cm  /V  sec, 
Na  is  the  doping  level  of  the  P substrate,  q is  the  electronic  charge, 

and  e is  the  dielectric  constant  of  the  silicon.  For  a typical 

13  2 -9 

doping  level  of  10  /cm  , this  gives  2.7  x 10  seconds. 

The  dominant  response  limitation  is  therefore  the  RC  time 

constant  of  the  diode  capacitance  and  the  resistance  of  the  N and 

P layers.  The  resistance  of  the  P layer  (back  side  of  the  target) 

can  be  made  negligible  by  a heavily  doped  back  surface  and  a back 
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metalization.  These  same  methods  can  be  applied  with  caution  to 

the  top  surface,  but  they  may  interfere  with  the  target  operation 

and  in  particular  introduce  a dead  layer.  The  resistance  of  the 

top  target  surface  without  a metal  layer  would  be  (10^/N  ) 

s 

2 

ohms /square  where  N is  the  surface  doping  in  dopant  atoms /cm  . 

s 

For  the  range  contemplated  for  this  device,  N would  be  between 

s 

12  l^f  2 

10  and  10  /cm  . This  gives  a top  surface  layer  resistance  of 
# 2 

between  10  and  10  ohms /square.  The  diode  capacitance  is  limited 
ultimately  by  the  avalanche  breakdown  of  the  diode  and  by  the  sub- 
strate doping.  For  100  ohm-cm  p substrates  operated  near  avalanche, 

2 

this  minimum  capacitance  limit  is  50  pF/cm  . Diodes  operating 
near  this  limit  are  now  used  in  similiar  devices  such  as  particle 
detectors,  electron  beam  bombarded  microwave  amplifiers  and  BEAMOS 
devices . 

The  target  frequency  response  can  be  modeled  as  an  RC  transmission 
line  where  the  resistance  and  capacitance  are  proportional  to 
1/R  and  R,  respectively,  where  R is  the  distance  from  the  target 
center  (Figure  7).  The  target  is  assumed  circular  and  the  beam 
is  assumed  incident  at  the  target  center,  for  simplicity.  The 
solution  for  the  external  target  current  (iQ(t))  for  an  impulse 
current  input  beam  current  pulse  (i^(t)  = 6 (t) ) is  shown  in 

Figure  8.  a is  the  target  radius,  r is  the  resistance  per  square 

+ 2 
of  the  N -layer  and  c is  the  diode  capacitance/cm  . The  response 
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time  is  approximately  0.3a  rc  or  the  frequency  response  cutoff 

2 2 
fc~3/a  rc.  For  a 2-cm-diameter  target  with  100-pF/cm  diode 

capacitance  and  1000  ohms /square  N+  layer,  this  gives  33  MHz, 

which  is  acceptable.  To  reduce  the  response  time  still  further, 

several  approaches  can  be  taken.  By  fabricating  the  target  segments 

smaller  or  in  a strip  geometry,  distance  a can  be  reduced  with  strong 

reductions  in  the  response  time.  Alternatively,  a thin  metal  layer 

or  heavily  doped  region  can  be  fabricated  on  the  top  surface.  400 

angstroms  of  aluminum  would  reduce  the  resistance  to  less  than 

10  ohms/square,  which  would  reduce  the  time  constant  to  less  than 

3 x 10  ^ seconds,  and  the  response  would  be  limited  by  carrier 

transit  times. 

Resolution 

The  resolution  of  the  readout  is  limited  by: 

• Size  of  the  written  P regions 

• Depletion  field  spreading 

• Lateral  diffusion  of  minority  carriers  in  the  N+  layer 

• Scattering  of  the  reading  electron  beam 

• Electron  beam  spot  size 

+ -f 

Size  of  Written  P (or  N ) Regions 

This  size  will  depend  on  the  writing  process  used.  It  will  be 
assumed  in  this  section  that  doped  regions  of  any  size  can  be 
produced  by  the  writing  process. 


I 


Depletion  Field  Spreadin 


The  built-in  voltage  of  the  P N junction  spreads  in  both 
directions  from  the  metallurgical  junction.  The  depleted  region 
spreads  until  the  number  of  ionized  dopants  in  the  depleted  region 
balances  the  built-in  voltage.  The  width  of  the  depleted  region 
in  the  abrupt  junction  approximation  is  given  by 


2 Gsi  NA  + ND  ®T 


q n n 
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where  is  the  dielectric  constant  of  the  semiconductor,  q is  the 
electronic  charge,  and  are  the  doping  levels  of  the  acceptor 
and  donors  in  the  P+  and  N+  sides  of  the  junction,  respectively, 
and  ©T  is  the  built-in  voltage  of  the  junction.  ®T  depends  on 
doping  levels,  but  for  heavily  doped  junctions  such  as  those  con- 
sidered here,  cd^,  = 0.9  volt.  For  a symmetric  junction,  and 

the  electric  field  distribution  is  shown  in  Figure  9,  where  x is 
measured  relative  to  the  matallurgical  junction;  em  = for  a 

symmetric  junction.  The  dependence  of  W on  for  = NQ  is  shown 

in  Figure  10.  Hence  the  effective  size  of  the  P+  doped  region  is 
D + W,  where  D is  the  diameter  of  the  doped  region.  It  should  also 

be  noted  that  the  N+  layer  must  be  thick  enough  so  that  the  depth 

4.  4. 

of  the  doped  P region  + W/2  does  not  exceed  the  N region  thick- 

4 4* 

ness,  or  the  small  P diode  would  short  to  the  P substrate  and  give 
high  (not  low)  readout  gain,  as  desired. 
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Figure 


Lateral  Diffusion  of  Minority  Carriers  in  N+  Laver 

During  reading,  the  injected  minority  carriers  (holes)  will 
diffuse  laterally  so  that  some  sampling  of  the  adjacent  bit  site 
will  occur.  The  amount  of  lateral  diffusion  depends  on  the  boundary 
conditions,  in  particular,  the  written  data  pattern  on  the  top 
surface.  Typically  the  amount  of  radial  lateral  diffusion  is  approx- 
imately the  N+  layer  thickness.  In.  BEAMOS,  the  N layer  is  nearly 
uniformly  doped,  so  that  no  driving  fields  exist  in  the  N layer 
except  near  the  top  and  diode  interfaces.  This  absence  of  fields 
is  necessary  because  of  the  method  of  BEAMOS  operation.  However  in 
the  Archival  application  it  is  possible  to  have  a built-in  field  in 
the  N+  layer,  to  aid  carrier  motion  towards  the  planar  diode  depletion 
region  and  thereby  reduce  laterial  diffusion.  This  built  in  electric 
field  also  reduces  carrier  recombination  at  the  top  silicon  surface. 
This  will  dramatically  improve  the  readout  gain  and  makes  the 
doping  gradient  or  built  in  field  desirable  even  if  it  is  not  neces- 
sary to  minimize  laterial  diffusion.  The  effect  of  the  built  in 
field  on  readout  gain  and  the  problem  of  laterial  diffusion  are 
discussed  in  more  detail  in  subsequent  sections. 

The  built  in  field  is  achieved  by  incorporating  a doping 
gradient.  The  field  due  to  a doping  gradient  outside  of  any 
depleted  regions  is  given  by 


E - agd 


(2) 


where  kT  is  the  thermal  energy  and  N(x)  = Np(x)  - N^(x). 


An  idealized  N+  doping  profile  that  could  be  well  approximated 

by  conventional  processing  techniques  is  shown  in  Figure  11.  The 

4 » 

electric  field  distribution,  which  is  a constant  10  V/cm,  is  shown 

in  Figure  12.  Figure  13  shows  an  acceptor  distribution  for  a 

written  zero  that  is  taken  as  Gaussian.  This  distribution  could  be 

obtained  by  thermal  diffusion  or  by  ion  implantation.  In  the  case 

of  ion  implantation,  this  particular  distribution  would  be  obtained 

for  ARp  = 0.025  micron  and  with  a top  metal  layer,  so  that  Rp  falls 

at  the  top  surface  of  the  silicon.  A possible  choice  in  this  case 

would  be  boron  ions  at  20  keV  or  gallium  at  140  keV  with  a metal 

layer  about  0.08  micron  in  each  case.  Figure  14  shows  the  net  doping 

profile  Nd  - obtained,  and  Figure  15  shows  the  approximate  electric 

field  distribution.  The  sharp  negative  electric  field  spike  around 

+ + 

0.04  micron  is  due  to  the  depletion  field  of  the  P N diode.  This 
can  in  principle  be  calculated  exactly  for  nonuniform  doping  profiles; 
however,  this  is  an  estimate  based  on  the  abrupt  junction  approxim- 
ation. Note  that  for  z<0.04  micron,  the  field  is  negative,  so  that 
this  region  is  a potential  well  for  holes  (minority  carriers  in  the 
N+  region).  The  effect  of  such  built  in  fields  on  carrier  motion 
is  discussed  in  the  dead  layer  modeling  section. 

Scattering  of  Reading  Electron  Beam 

The  distribution  of  the  ionization  due  to  a low-energy  electron 
beam  is  shown  in  Figure  16.  (Appendix  II,  "Electron  Penetration  in 
Solids,"  discusses  this  subject  in  detail.)  Figure  16  is  a normalized 
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Figure  11. 


Idealized  Doping  Distribution  for  N+  Layer  of  Archival 
Target. 
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plot  of  equal  ionization  intensity  contours  for  a well  focused 
beam  incident  on  a single  point,  as  a function  of  penetration  z and 
lateral  distance  r.  This  distribution  applies  between  at  least  2 
and  10  keV  for  low-Z  solids  such  as  silicon  and  probably  germanium. 

The  distance  scaling  factor  is  R^,  the  so  called  Grun  range. 

The  dependance  of  R_.  on  beam  energy  for  silicon  and  germanium  is 
shown  in  Figure  17.  RG  scales  as  the  mass  density.  Note  that  most 
of  the  beam  energy  is  contained  within  a hemisphere  of  radius  RG, 
so  that  Rg  represents  a good  approximation  to  the  upper  limit  of 
the  effective  increase  in  diameter  of  the  electron  beam  as  it 
penetrates  the  material.  For  a finite  spot  size  that  is  of  the 
same  size  or  larger  than  the  Grun  range,  integration  of  the  ioniz- 
ation distribution  for  a small  spot  yields  an  ionization  distribution 
in  the  material  that  is  apparently  sharper  than  would  appear  from 
Figure  16.  This  is  the  case  because  electrons  in  the  center  of  the 
beam  cannot  scatter  far  enough  to  reach  the  edge  of  the  ionization 
distribution,  and  some  electrons  at  the  edge  scatter  into  the  center. 

Figure  18  shows  the  2-keV  ionization  distribution  at  the  edge  of  a 
spot  whose  diameter  is  larger  than  0.04  micron  for  three  depths,  z, 
in  silicon.  This  curve  is  for  the  case  in  which  the  beam  scans  in 
a line  rather  than  dwells  on  one  spot  and  is  based  on  the  measure- 
ments of  Possin  and  Norton.  Note  that  the  broadening  of  the  beam 
is  < 0.02  micron  in  radius.  Hence  for  a 0.4-micron  spot  at  2 keV, 
a cylinder  of  ~0. 08-micron  diameter  and  a 0.06-micron  depth  should 

ii 
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Figure  18.  Ionization  per  Unit  Volume  in  Arbitrary  Units  as  a 
Function  of  Lateral  Distance  x and  Penetration  z 
for  Sharp  Edged  Beam  of  Diameter  Greater  than  0.06 
Micron.  (Calculated  from  5-kV  Beam  Penetration  Data 
in  FMM;  Scaled  to  2-keV  and  Silicon  Mass  Density). 
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contain  80  percent  of  the  beam  energy  for  silicon.  For  other  low-Z 
materials,  both  distances  can  be  scaled  inversely  as  the  density. 

For  targets  constructed  of  germanium,  the  situation  is  there- 
fore even  more  favorable.  The  Grun  range  is  smaller  by  a factor 
of  two  at  the  same  voltage  (because  of  the  increased  mass  density), 
and  the  electron  hole  pair  production  rate  is  30  percent  larger 
(because  of  the  small  band  gap).  Hence  higher  energy  reading  beams 
can  be  used  to  give  the  same  size  ionization  region,  but  with  a 
target  gain  larger  by  a factor  of  two. 

The  above  discussion  has  indicated  the  factors  to  be  considered 

in  fabricating  surface  diode  targets.  Figure  19  shows  a possible 

target  with  a magnitude  of  all  of  these  quantities  indicated  to 

scale  for  bits  recorded  on  0.1-micron  centers.  The  two  P+  areas 

are  on  0.2 -micron  centers,  which  correspond  to  a 010  pattern  on 

0.1-micron  centers.  Figure  19  shows  a one  surrounded  by  zeros, 

which  is  the  most  difficult  pattern  for  this  device  to  read.  The 

and  Xp  are  the  penetration  of  the  diode  depletion  region  into  the 

N+  and  P layers,  respectively.  For  this  figure  it  is  assumed  that 

the  P substrate  is  400-ohm-cm  material ; the  N+  doping  profile  is 

shown  in  Figure  11.  The  main  N+P  diode  is  shown  back  biased  to  100 

volts.  As  can  be  seen,  nothing  but  lateral  carrier  diffusion  pre- 

•f  4* 

vents  achieving  0.1-micron  readout  resolution.  The  P N diode 
depletion  field  is  based  on  the  doping  levels  in  Figure  14.  The 
depletion  width  could  be  reduced  even  more  by  utilizing  higher  doping 
levels  in  the  P+  and  N+  regions. 
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Figure  19.  Cross  Section  of  101  Pattern  Recorded  on  Surface  Diode  Target  (All  Dimensions 
Shown  to  Scale) . 


Based  on  these  considerations  it  will  be  the  goal  of  the 

fabrication  effort  to  produce  planar  diodes  with  doping  levels 
18  -3 

~5  x 10  cm  near  the  surface  and  with  doping  profiles  like 
figure  11.  The  primary  reason  for  the  strong  doping  gradients 
is  to  produce  built  in  electric  fields  which  will  increase  target 
readout  gain  by  minimizing  carrier  recombination  at  the  top  surface. 
This  is  very  important  because  the  effect  of  writing  on  the  planar 
diode  is  to  increase  carrier  recombination  in  the  top  surface  region. 
To  maximize  the  readout  modulation  the  unwritten  target  state  must 
have  minimum  recombination  in  this  top  surface  layer.  High  doping 
levels  near  the  surface  are  also  necessary  to  prevent  depletion 
field  spreading  of  the  written  diodes  and  loss  of  spatial  resolution 
or  bit  packing  density.  In  addition,  heavily  doped  surface  layers 
increase  the  electrical  conductivity  of  the  n+  layer  and  thereby 
increase  the  frequency  response  of  the  readout  signal. 

Fabrication  of  Planar  Diodes 

The  planar  diode  studies  consisted  of  a number  of  implantations 
with  differing  processing  parameters.  The  diode  structures  were 
examined  in  an  effort  to  determine  the  optimum  conditions  for 
producing  the  desired  doping  with  a minimal  dead  layer  region  at 
the  surface.  The  influence  of  ion  species,  fluence,  and  energy 
as  well  as  substrate  crystallographic  orientation,  surface  oxide 
layers,  and  annealing  temperature  were  investigated  as  a part  of  a 
continuing  effort  in  this  area. 
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The  "dead"  layer  referred  to  is  a short -lifetime  region  near 
the  top  surface  of  the  target  which  appears  dead  to  the  reading 
electron  beam.  Carriers  generated  by  the  beam  in  this  area  recom- 
bine at  the  many  nearby  trapping  sites,  with  little  chance  to  diffuse 
down  to  the  collecting  junction  of  the  planar  diode.  This  results 
in  low  gain  for  the  device  at  low  beam  voltage  and  impairs  the 
device  resolution.  Therefore,  it  is  desirable  to  make  the  dead 
layer  as  thin  as  possible. 

The  dead  layer  phenomena  has  been  observed  by  other  researchers 
working  in  the  area  of  solar  cells  and  nuclear  particle  detectors. 

Dead  layers  are  quite  prevalent  for  very  heavy  phosphorus  diffusions. 
It  has  been  speculated  that  the  presence  of  large  numbers  of  oversize 
atoms  in  the  lattice  can  produce  sufficient  disorder  to  cause  the 
region  to  appear  dead.  In  the  case  of  nuclear  particle  detectors, 
for  example,  extreme  care  is  used  in  preparing  the  junctions.  Implan- 
tations at  very  low  energies,  5-10  kV,  directly  onto  the  crystallo- 
graphic axis  results  in  the  ions  being  gently  nudged  into  the 
crystal  to  produce  as  little  damage  as  possible  and  to  confine  the 
damage  to  a thin  surface  region.  (Ref.  1).  Since  the  archival 
memory  device  requires  relatively  heavy  dopant  concentrations  to  a 
depth  of  ~0.1u  this  kind  of  approach  is  not  a feasible  solution  for 
the  planar  diode  substrate. 

Solar  cell  researchers  have  had  some  success  in  producing 
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adequately  thin  dead  layers  for  solar  cell  applications  by  careful 
phosphorus  diffusions  and  special  surface  preparations.  However, 
even  with  these  procedures  the  dead  layers  are  often  unacceptable 
for  our  application  as  the  dead  region  is  around  500  8.  (Ref.  2). 

The  approach  pursued  in  these  studies  instead,  has  been  to 
achieve  the  necessary  doping  levels  and  gradients  through  the  use 
of  ion  implantation,  and  by  selecting  the  appropriate  combination 
of  ion,  energy,  fluence,  crystal  orientation,  surface  preparation, 
and  annealing  cycle,  to  restore  a high  quality  crystalline  structure 
with  a minimal  dead  layer  at  the  top  surface  of  the  device.  The 
influence  of  a number  of  these  parameters  has  already  been  investig- 
ated. Some  of  these  will  be  looked  into  further  while  work  on 
others  not  yet  begun  will  start  in  the  upcoming  quarter. 

For  this  initial  work  the  most  widely  accepted  implantation 
dopants,  phosphorus,  arsenic,  and  boron,  were  used  to  form  planar 
diode  structures  under  a number  of  implantation  conditions  in 
~5-20  A -cm  float  zone  silicon  using  the  400  KeV  Extrion  ion  implan- 
ter  at  the  Center.  The  popularity  of  these  column  III  and  IV 
elements  is  due  in  large  part  to  the  known  high  degree  of  electrical 
activity  of  these  ions  following  implantation  and  annealing.  The 
ready  availability  of  these  ions  due  to  their  ease  of  implantation 
made  them  ideal  first  candidates.  It  is  foreseen  that  the  suitabil- 
ity of  other  less  common  ions  will  be  evaluated  in  the  future  as  well. 
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All  of  these  Ions  weie  Implanted  at  room  temperature  at  a beam 
angle  7°  from  the  crystallographic.  It  was  felt  that  at  present, 
variation  of  the  implant  temperature  would  not  really  be  profitable. 
Hot  implants,  for  example,  would  result  in  less  damage  immediately 
after  implantation.  However,  higher  temperatures,  and  thus  spread- 
ing of  the  profile,  would  be  needed  to  anneal  out  residual  radiation 
damage.  Cold  implants,  on  the  other  hand,  would  result  in  a greater 
degree  of  channeling,  which  might  be  undesirable  in  view  of  the 
highly  uniform,  shallow  junction  being  sought.  For  similar  reasons 
no  attempts  at  angular  variation  were  attempted  to  utilize  channeling 
to  an  advantage,  since  variation  in  profile  and  changes  in  surface 
conditions  can  result. 

Energies  of  30  kV  and  100  kV  were  selected  for  the  first 
studies  to  offer  contrasting  situations  for  investigation  of  dead 
layer  formation.  At  30  kV,  as  shown  in  Figure  20,  the  depth  of 
ion  penetration,  and  the  resulting  damage  profile,  is  located  much 
closer  to  the  top  surface  of  the  device  than  at  100  kV.  Although 
the  lower  energy  implant  would  put  the  peak  of  the  dopant  distrib- 
ution near  the  surface  as  desired,  the  higher  energy  implant  could 
still  be  used  to  form  the  profile  shown  in  Figure  11,  by 
implanting  through  a layer  of  another  material  such  as  oxide  or 
metal,  or  by  layer  removal  at  the  surface  down  to  the  peak  of  the 
distribution.  Either  of  these  schemes  would  allow  the  production 
of  a top  highly  doped  region  with  a steep  gradient  through  the  upper 
n (or  p)  layer. 
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Figure  20.  Projected  Ranges  for  P+,  As  and  B+  at  30  and  100  kV. 
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Samples  with  200  A oxides  were  implanted  with  PI  , As  , and  B 

ions  for  comparison  with  polished  wafer  surfaces.  Several  oxide 
thicknesses  which  would  place  the  peak  of  the  implanted  ion  distrib- 
ution near  the  surface  were  also  calculated  for  each  ion  at  100  kV 
and  the  samples  were  implanted.  However  all  of  these  samples  have 
not  been  examined.  Yet  under  consideration  is  the  use  of  metal  top 
layers  in  a similar  way. 

A literature  survey  of  etches  which  remove  radiation  damaged 
material  in  a highly  controllable  and  uniform  manner  was  undertaken. 
Several  of  these  etching  techniques  will  be  tried  in  the  next  quar- 
ter. This  could  be  particularly  valuable  if  the  high  temperature 
anneals  prove  necessary  for  low  diode  noise  and  leakage  and  result 
in  thicker  dead  layers  as  mobile  defects  rise  to  the  surface.  Prelim- 
inary work  prior  to  the  proposal  indicated  that  some  of  the  dead  layer 
could  be  removed  by  etching.  However,  a residual  dead  layer  thick- 
ness, which  could  be  reduced  no  further,  remained.  Similar  results 
were  observed  by  Guldbery  & Schroder.  (Ref.  3).  Various  crystall- 
ographic orientations  were  used  for  the  work  to  determine  if  a 
thinner  dead  layer  would  result  for  any  particular  direction.  P-type 
wafers  with  <111>,  <110>  and  <110>  orientations  were  used  for  P+  and 
As  implants.  <110>  and  <111>  wafers  received  B+  implants. 

Other  research  has  produced  evidence  that  significant  differ- 
ences can  occur  in  the  regrowth  of  highly  damaged  regions  during 
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anneal  depending  on  the  orientation  selected.  (Ref.  4)  (Ref.  5). 

In  particular,  they  have  observed  more  residual  damage  in  <111>  Si. 

14 

However,  for  lower  doses  <10  , little  difference  has  been  noticed 

(Ref.  6)  in  the  annealing  of  damaged  regions.  Figure  21  shows  the 
regrowth  of  heavily  damaged  Si.  The  <111>  sample  did  not  completely 
anneal . 

The  sharpness  of  the  profile  can  also  be  altered  by  the  crystal 
orientation.  Supertailing,  thought  to  occur  as  a result  of  anomalous 
diffusion  effects  or  scattering  into  channels,  shows  a marked  varia- 
tion with  sample  direction  (Fig.  22).  The  effect  here,  makes  <111> 
Si  take  on  a sharper  dopant  gradient. 

Ion  fluence  is,  of  course,  also  known  to  have  considerable 
effect  on  the  amount  of  damage  a crystalline  material  such  as 

silicon.  For  the  doping  levels  required  for  the  planar  diode 

13  14 

structure,  fluences  of  10  or  more  likely  10  are  needed,  and  thus 
both  were  studied.  The  device  described  in  Section  IV  however, 
may  not  require  this  high  a doping  level.  These  fluences  are  not 
high  enough,  in  general,  to  drive  the  crystal  amorphous.  Amor- 
phousness, although  indicative  of  very  heavy  lattice  damage,  allows 
regrowth  with  high  electrical  activities  at  significantly  lower 
annealing  temperatures.  (Ref.  5).  Pre-implants  of  heavy  fluences 
(~10^)  of  inert  ions  are  planned  to  determine  if  this  will  produce 

higher  quality  diodes  for  low  temperature  anneals.  A preliminary 

15  -2  + 

experiment  using  10  cm  of  He  ions  for  post  damage  before  anneal 
was  not  successful  in  reducing  the  dead  layer. 
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Figure  22.  Comparison  of  depth  profiles  in  <111>  and  <110>  Silicon 
for  20  keV  P implants  at  10^  cm"^ . Curve  labeled  total 
is  the  total  P concentration  (electrically  active  plus 
active)  measured  by  the  radioactive  tracer  method  from 
Ref.  7 for  implantation  on  axis  into  <110>  after  anneal- 
ing at  400°C.  Other  curves  are  the  electrically  active 
P dopant  measured  by  stripping  and  four  point  probe, 

(Ref.  8),  after  annealing  at  420°C.  The  <110>  implant 
was  on  axis  and  the  <11 1>  implants  were  on  axis  and  15° 
off  axis.  Note  that  the  profiles  are  the  same  for  on  and 
off  axis  implants  into  <111>. 
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Probably  the  most  dominant  effect  on  the  amount  of  residual 
damage  in  an  ion  implanted  device  is  the  annealing  conditions  to 
which  the  structure  is  subjected.  Such  annealing  can  take  place 
during  as  well  as  after  the  implantation  process.  Heating  due  to 
high  dose  of  a sample  not  well  heat  sinked  can  cause  significant 
annealing.  Even  for  low  fluences,  lower  temperature  annealing 
stages  do  exist,  ~650°C  for  arsenic  and  phosphorous,  and  ~900°C 
for  boron.  Regions  near  these  stages  were  therefore  studied. 

Isochronal  annealing  studies  were  undertaken  to  determine  more 
precisely  the  temperature  needed  to  produce  the  desired  device  char- 
acteristics. In  addition,  more  detailed  isochronal  as  well  as 
perhaps  isothermal  annealing  studies  will  be  performed  to  determine 
the  optimum  temperatures  and  times  required  for  full  recovery  of 
radiation  damage  to  minimize  the  dead  layer  and  maximize  activation 
of  the  implanted  electrical  centers  for  satisfactory  diode  perfor- 
mance. 

The  lack  of  suitable  equipment  during  the  first  quarter  has 
precluded  the  investigation  of  the  device  potential  of  germanium. 
Particularly  due  to  its  possible  application  in  the  alloy  junction 
device,  this  material  requires  additional  study.  Studies  of  a type 
similar  to  those  aforementioned  can  be  undertaken  for  this  substrate 
material  during  the  next  quarter. 

Eva lua t ion  of  Planar  Diodes 


The  planar  diodes  fabricated  as  described  in  the  preceeding 
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pages  of  this  section  were  evaluated  in  the  Cambridge  Mark  II  Scanning 
Electron  Microscope  (SEM) . For  these  measurements  the  instrument 
is  used  in  a diode  current  imaging  mode  as  indicated  in  Figure  23. 

The  specimen  current  amplifier  (SCA)  of  the  SEM  is  used  to  amplify 
the  diode  current  signal  which  is  then  used  to  modulate  the  z-axis 
of  the  SEM  displays  and  is  also  displayed  on  monitor  oscillscope. 

The  diode  back  bias  (V  diode)  is  normally  1.5  volts.  The  electron 
beam  induced  gain  of  the  target  is  defined  as 

Gain  * G = Ig/IB  (3) 

where  I is  the  increase  in  diode  leakage  with  the  beam  on  and  I_ 

S D 

is  the  beam  current.  All  beam  currents  are  measured  using  a 
Faraday  cup  and  the  SCA  operated  in  the  electrometer  mode.  The 
beam  energy  is  measured  using  a calibrated  electrostatic  voltmeter 
attached  directly  to  the  filament  lead  of  the  SEM. 

The  electron  beam  induced  gain,  G,  is  measured  as  a function  of 
beam  energy.  The  complete  set  of  gain  measurements  on  all  targets 
studied  so  far  is  contained  in  Figures  24  through  59.  All  of  the 
gain  data  has  been  plotted  against  the  first  order  dead  layer  model 
to  be  described  later  in  this  section.  As  will  be  explained,  there 
is  no  reason  to  expect  the  experimental  results  to  follow  this  simple 
model.  However  it  does  help  in  visualizing  the  spatial  dimensions 
involved,  and  in  providing  a criteria  for  comparing  targets.  For 
example,  target  sample  P3A  in  Figure  25  can  be  said  to  behave  at 
2.5  keV  as  though  it  had  a dead  layer  of  thickness  T^  = 0.06  microns 
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Figure  23.  Specimen  Current  Amplifier  Operated  in  Diode  Current 
Imaging  Mode 
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voltage  measurements  for  12-18  ohm-cm,  <111>  p-Si,  implanted  , 
temperature  7°  off  axis  at  100  kV  to  a fluence  of  1013  ions/cm' 


ERG 


BEAM  ENERGY  (keV) 


Figure  28b.  Gain  vs.  voltage  measurements  for  12-18  ohm 
^Hl^j  p-Si,  implanted  at  room  temperature  ] 
axis  at  100  kv  to  a fluence  of  lO1^  ions /cm' 
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Figure  33.  Gain  vs.  voltage  measurements  for  12-18  ohm-cm, 
<111>  p-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10^  ions/cm^. 
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Figure  36.  Gain  vs.  voltage  measurements  for  7-10  ohm-cm, 

<111>  p-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10*^  ions/cm^. 
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Figure  38.  Gain  vs.  voltage  measurements  for  4. 

p-Si,  implanted  at  room  temperature 
100  kV  to  a fluence  of  10  ^ ions/cm^ 
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7°  off  axis  at 


1400 


GAIN 


Figure  44.  Gain  vs.  voltage  measurements  for  12-18  ohm-cm,  <111N 
p-Si,  implanted  at  room  temperature  7°  off  axis  at 
100  kV  to  a fluence  of  lO1^  ions /cm2,  through  1400  ! 
oxide . 
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Figure  45.  Gain  vs.  voltage  measurements  for  12-18  ohm-cm,  <111>  p-Si,  implanted  , 
at  room  temperature  7°  off  axis  at  100  kV  to  a fluence  of  10^  ions/cm' 
through  1400  X oxide. 


1400 


CM 

T3  O 
<U 


NIV9 


GAIN 


1200 


BEAM  ENERGY  (keV) 


Figure  51b.  Gain  vs.  voltage  measurements  for  12-18  ohm-cm, 
<111>  p-Si,  implanted  at  room  temperature  7° 
off  axis  at  30  kV  to  a fluence  of  lO1^  ions/cm2 
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Figure  52.  Gain  vs.  voltage  measurements  for  4. 8-7. 2 ohm-cm 
<100>  p-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10^  ions/cm^. 
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Figure  53.  Gain  vs.  voltage  measurements  for  4. 8-7. 2 ohm-cm,  <100'>  p-Si,  implanted 
at  room  temperature  7°  off  axis  at  100  kV  to  a fluence  of  10*  ions/cm 
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Figure  55.  Gain  vs.  voltage  measurements  for  12-18  ohm-cm, 
<111>  p-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  lo'^  ions/cm^, 
through  200  X oxide. 
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Figure  56.  Gain  vs.  voltage  measurements  for  11-13  ohm-cm, 
<111>  p-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10^  ions/cm^. 
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Figure  58.  Gain  vs.  voltage  measurements  for  11-13  ohm-cm, 
<111>  n-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10^  ions/cm^, 
through  200  K oxide. 
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Figure  59.  Gain  vs.  voltage  measurements  for  11-13  ohm-cm, 
<111>  n-Si,  implanted  at  room  temperature  7°  off 
axis  at  100  kV  to  a fluence  of  10^  ions /cm  . 


and  at  5 keV  as  though  = 0.12  microns.  This  provides  a convenient 
method  of  comparing  targets  while  at  the  same  time  giving  some  estim- 
ate of  the  thickness  of  the  dead  layers  involved.  To  avoid  confusion 
the  term  dead  layer  will  be  used  in  the  generic  sense  of  a phenomena 
which  causes  reduced  electron  beam  gain  from  planar  diodes.  T^  will 
be  used  to  describe  the  dead  layer  thickness  in  the  first  order 
model.  Tp  will  generally  require  the  measurement  voltage  to  also 
be  stated,  when  it  is  being  used  to  describe  experimental  data. 

In  addition,  many  of  the  targets  were  scanned  and  photographed 
to  search  for  gain  irregularities.  The  gain  irregularity  was  usually 
less  than  57,.  On  some  of  the  gain  curves,  the  error  bars  are  included 
indicating  the  amount  of  nonuniformity  observed.  Figure  60  shows 


such  a gain  photograph.  The  bright  oval  is  the  diode  area  and  the 
dark  central  area  is  the  area  of  the  diode  shadowed  by  the  contact 


finger  and  Faraday  cup.  The  upper  part  of  the  diode  has  57  higher 

§- 

gain.  This  area  had  an  oxide  layer  present  during  the  anneal.  Also 


note  the  gain  variations  over  the  entire  diode.  These  are  also 
less  than  57,.  These  small  variations  show  up  clearly  in  this  picture 
because  the  contrast  is  highly  expanded. 

Tables  1,  2 and  3 summarize  the  gain  measurements  on  the 
targets  examined  so  far.  All  samples  were  implanted  at  7°  off  the 
crystal  axis. 

The  best  results  were  obtained  for  samples  P11A  and  P11C 
(Figures  43,  44,  45)  and  A9C  and  A7B  (Figures  50  & 57).  For  all 
of  these  samples  the  implantation  peak  is  at  or  very  near  the 
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Figure  60.  SEM  Diode  Signal  Micrograph  Showing  Gain  Difference  For 
Substrate  A9B  After  Anneal  at  650°  With  Oxide  Left  on 
upper  half  of  Mesa  During  Anneal  Cycle  and  Removed 
Before  Measurement.  Beam  Energy  10  kV  and  Magnification 
~ 17x. 
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Table  1 

DEAD  LAYER  THICKNESSES  MEASURED  AT  2.5  and  5 keV  FOR  P+  IMPLANTED 
p-SI.  ALL  ANNEALS  1 HR.  IN  Ar  ATMOSPHERE.  SUBSTRATE  RESISTIVITIES 
PRIOR  TO  IMPLANT  5-18a  - CM. 


Sample 

Number 

Ion 

Energy 

(keV) 

Ion 

Fluence 

(0/cm2) 

Crystal 

Orient. 

4. 

Surface 

Prep. 

Anneal 

Temp. 

<°C) 

Dead  Layer 
Thickness  (A) 

at  2.5  keV 

at  5 keV 

P1A 

100 

1013 

<111  > 

20  min. 
at  550 

800 

1400 

P1B 

100 

io13 

<111  > 

770 

650 

1550 

P3A 

100 

io14 

<111  > 

770 

600 

1200 

P3B 

100 

1014 

<111  > 

575 

800 

1700 

675 

660 

1200 

770 

660 

1500 

io14 

860 

660 

1500 

P3C 

100 

<111  > 

650 

700 

1250 

P4C 

100 

10 14 

<111  > 

770 

660 

1300 

P4A 

100 

io14 

<111  > 

770 

220 

420  tr! 

P4B 

100 

io14 

<111  > 

~1050 

700 

~ 2000 

P5A 

100 

io14 

<110  > 

770 

630 

1300 

P5B 

100 

io14 

<no  > 

650 

700 

1400 

P6A 

100 

io14 

<100  > 

770 

630 

1400 

P6B 

100 

io14 

<100  > 

650 

NE 

NE 

P7A 

30 

io14 

<111  > 

770 

600 

1100 

P7B 

30 

io14 

<111  > 

650 

580 

1300 

P8C 

100 

io14 

<111  > 

200  A 

770 

620 

1150 

P8B 

100 

10 14 

<111  > 

650 

NE 

NE 

P9C 

100 

io14 

<111  > 

1000  A 

oxide 

770 

440 

900 

P9D 

100 

io14 

<111  > 

1000  X 

oxide 

650 

NE 

NE 

PIOC 

100 

io14 

<111  > 

1200  A 

oxide 

770 

430 

720 

P10A 

100 

10 14 

<111  > 

1200  A 

oxide 

650 

NE 

NE 

P11C 

100 

io14 

<111  > 

1400  A 
oxide 

770 

380 

600 
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Table  1 (con't) 


Sample 

Number 

Ion 

Energy 

(keV) 

Ion 

Fluence 

(#/cm2) 

Crystal 

Orient. 

Surface+ 

Prep. 

Anneal 

Temp. 

(°C) 

Dead  Layer 
Thickness  (X) 

P11A 

100 

io14 

<111  > 

1400  k 

650 

at  2.5  keV  at 

400 

5 keV 

800 

P12C 

100 

io14 

<111  > 

oxide 

770 

660 

1400 

P13B 

100 

io14 

<111> 

650 

800 

1800 

50  He 

io15 

<111> 

NE  = Not  Evaluated 

+A11  wafers  are  polished,  bright  etched  and  cleaned  in  BHF  before  implant 
(except  oxide  wafers) 
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Table  2 


DEAD  LAYER  THICKNESSES  MEASURED  AT  2 . 5 AND  5 keV  FOR  As+  IMPLANTED 
p-Si.  ALL  ANNEALS  1 HR.  IN  Ar  ATMOSPHERE.  SUBSTRATE  RESISTIVITIES 
PRIOR  TO  IMPLANT  5 - 15 A - CM. 


Sample 

Number 

Ion 

Energy 

(keV) 

Ion 

Fluence 

(#/cm2) 

Crystal 
Orient . 

Surface+ 

Prep. 

Anneal 

Temp. 

(°C) 

Dead  Layer 
Thickness  (1) 

at 

2.5 

keV  at  5 keV 

A3A 

100 

1014 

<111> 

770 

520 

700 

A3B 

100 

1014 

<111> 

650 

NE 

NE 

A4B 

100 

1014 

<110> 

770 

480 

700 

A4C 

100 

1014 

<110> 

650 

NE 

NE 

A7B 

30 

1014 

<111> 

770 

320 

400 

A7C 

30 

1014 

<111> 

650 

380 

700 

A8C 

100 

1014 

<100> 

770 

570 

1100 

A8A 

100 

1014 

<100> 

650 

700 

1400 

A9C 

100 

1014 

<111> 

2001 

oxide 

770 

320 

400 

A9B 

100 

1014 

<111> 

2001 

oxide 

650 

540 

1200 

NE  * Not  Evaluated 

+A11  wafers  are  polished,  bright  etched  and  cleaned  in 
BHF  before  implant  (except  oxide  wafers) . 


Table  3 


DEAD  LAYER  THICKNESSES  MEASURED  AT  2.5  AND  5 keV  FOR  B+  IMPLANTED 
n-Si.  ALL  ANNEALS  1 HR.  IN  Ar  ATMOSPHERE.  SUBSTRATE  RESISTIVITIES 
PRIOR  TO  IMPLANT  4 - 13A  - CM. 


Sample 

Number 

Ion 

Energy 

(keV) 

Ion 

Fluence 

(#/cm2) 

Crystal 

Orient. 

Surface* 

Prep. 

Anneal 

Temp. 

(°C) 

Dead  Layer 
Thickness  (l) 

B3A 

100 

10 14 

<111> 

200 X 

at 

950 

2.5  keV  at  5 keV 

2300 

B4A 

100 

10 14 

<111> 

oxide 

950 

1200  1900 

B5C 

100 

io14 

<110> 

950 

(3  keV) 

NE  NE 

NE  = Not  Evaluates 


+A11  wafers  are  polished,  bright  etched  and  cleaned 
in  BHF  before  implant  (except  oxide  wafers). 
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surface.  Pll  and  A9  used  implantation  through  oxide  layers  to 

achieve  a implant  peak  near  the  surface.  A7  used  30  keV  As+ 

implantation  to  achieve  a very  shallow  penetration  of  the  ions  (see 

Figure  20).  Note  also  that  650°C  annealing  for  the  As  samples  gives 

o 

consistently  poorer  results  than  760  C anneals.  Higher  temperature 
anneals  are  planned  for  the  As  samples.  For  the  P samples  650  and 
760°C  anneals  gave  about  the  same  results.  The  Boron  implanted 
samples  had  much  too  deep  implant  peaks  for  good  results.  Boron 
implants  through  thick  oxides  to  put  the  implant  peak  near  the 
surface  are  in  progress.  It  is  desirable  to  have  processing  tech- 
niques  for  both  n and  p surface  planar  diodes  so  that  both  n and 
p dopants  can  be  used  for  writing.  Also  damage  writing  may  give 
better  results  in  n or  p surfaces  depending  upon  the  type  of  deep 
traps  introduced  by  the  damaging  ion  beam.  The  extremely  good 
result  for  sample  P4A  (Figure  32)  could  not  be  reproduced  and  does 
not  fit  the  pattern  of  the  rest  of  the  samples.  The  only  identifiable 
difference  is  that  it  was  annealed  in  pre -purified  Ar  without  any 
trapping.  All  other  anneals  used  a dry  ice  and  acetone  trap  to 
additionally  clean  the  Ar.  One  other  sample,  P3A,  was  also  annealed 
without  a trap  in  an  attempt  to  reproduce  P4A.  P3  and  P4  are  from 
different  but  nominally  identical  silicon  wafers  implanted  sequentially 
in  the  same  implantation  run.  Note  also  that  P4B  and  P4C,  which 
received  similar  annealing  treatments  do  not  have  dead  layers  nearly 
as  small  as  P4A.  P4A  was  anamolous  also  in  that  it  occassionally 


1 


gave  much  lower  gain  indicating  a thicker  dead  layer.  This  occurred 
several  times  during  separate  measurements  of  the  sample  by  different 
operators  on  different  days.  One  of  these  data  points  is  indicated 
in  Figure  32a.  These  low  gain  data  points  could  never  be  consistently 
reproduced . 

The  good  results  for  the  Pll  and  A7  and  A9  samples  show  that 
acceptable  dead  layers  can  be  achieved  for  n+  surface  planar  diodes. 
These  targets  achieved  gains  in  excess  of  200  at  2 keV  which  is  507, 
of  the  maximum  possible.  Since  it  is  expected  that  beam  energies 
~2  keV  will  be  optimum  for  0.1  micron  bit  spacing  readout,  these 
targets  are  very  acceptable  for  ion  implantation  writing  and  readout. 
For  the  alloy  junction  writing  method  thinner  dead  layers  may  be 
necessary  as  discussed  in  the  alloy  junction  section. 

Figure  61  shows  an  interesting  effect  called  a bright  edge. 

This  is  a diode  gain  picture  as  in  Figure  60.  The  high  gain  strip 
along  the  left  edge  of  the  picture  corresponds  to  the  edge  of  the 


Figure  61. 


S KM  Piode  Signal  Micrograph  of  Sample  P10C  After 
0n(-  mul-  Anneal  at  770°  Showing  Bright  Ldge  (High 

Gain  Area)  at  Left. 
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however,  usually  it  extends  only  50  microns  as  in  Figure  61.  If  no 
inversion  layer  is  present  the  extent  of  the  bright  edge  is  a measure 
of  the  minority  carrier  lifetime  in  the  substrate.  A minority 
carrier  diffusion  length  of  100  microns  is  very  typical  for 
average  quality  silicon.  In  several  of  the  Figures  (32,  35,  45, 

50,  51,  57)  the  maximum  gain  in  the  bright  edge  is  shown  as  open 
circles.  Note  that  in  most  cases  this  gain  exactly  follows  the 
TD  = 0 dead  layer  curve.  This  means  that  all  the  carriers  generated 
by  the  electron  beam  are  being  collected  and  none  are  lost  in  a 
dead  layer.  This  indicates  that  the  dead  layer  phenomena  is  not 
an  intrinsic  property  of  silicon  surfaces  but  is  related  to  the 
introduction  of  dopant  near  the  surface.  Unfortunately  high  doping 
levels  near  the  surface  are  necessary  to  achieve  high  resolution  as 
discussed  earlier  so  that  ways  to  minimize  the  effect  in  the  presence 
of  high  doping  must  be  found. 

Dead  Laver  Models 

For  optimum  performance  of  the  surface  diode  targets  it  is 
necessary  to  fabricate  planar  diodes  with  very  thin  dead  layers. 

This  has  been  discussed  earlier  in  this  section.  By  a dead  layer  is 
meant  region  in  which  carriers  generated  by  the  electron  beam  are 
not  collected  by  the  diode.  Figure  50  shows  a typical  electron  beam 
induced  diode  gain  versus  beam  energy  curve  for  a planar  diode  target. 
The  straight  line  labeled  0 microns  is  what  would  be  expected  if  all 
the  carriers  generated  by  the  electron  beam  were  collected  by  the 


Pn  junction.  The  simplest  explanation,  for  deviations  from  this 

ideal  curve  is  that  a layer  of  the  silicon  near  the  surface  has 

a large  amount  of  damage  so  that  the  minority  carrier  lifetime  is 

very  short  and  the  generated  carriers  recombine  before  reacting 

the  junction.  Other  explanations  are  possible  as  will  be  seen  from 

the  following  discussion. 

1.  First  Order  Dead  Layer  Model. 

Figure  63  shows  a first  order  model  for  the  dead  layer.  We 

will  consider  n on  p diodes  for  def inetiness.  The  n layer  is 

assumed  to  be  uniformly  doped  and  have  a minority  carrier  lifetime 

much  larger  than  the  duffusion  time  across  the  n-layer  T^/D^  where 

T__  is  the  n-layer  thickness  and  D the  hole  diffusion  constant 
N p 

~ 2 

(D^  = io  cm  /sec).  All  carriers  generated  at  x<TD  are  not  collected 
by  the  diode  and  all  carriers  generated  for  z>TD  are  collected. 

This  includes  carriers  generated  at  z>T^.  No  assumptions  need  be 
made  about  the  reasons  for  these  collection  characteristics.  The 
generation  of  carriers  by  the  penetrating  electron  beam  is  well 


described  by: 


= »-(*&>) 


where  g(z,E)  is  the  number  of  electron-hole  pairs  generated  per 
3 2 

cm  per  incident  electron/cm  , E is  the  incident  electron  beam 
energy,  f the  fraction  of  electron  beam  energy  backscattered  by  the 
target,  Mq)  is  the  Grun  generation  function  and  R(E)  is  the  Griin 
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range  (see  Appendix  II).  For  A1  or  Silicon  f = 0.08.  This  should 
be  distinguished  from  the  number  backscatter  coefficient  which  for 
A1  or  Si  is  0.13.  For  silicon  or  silicon  dioxide: 

R(E)  « 0.019  microns  x E1'75  (5) 

where  E is  in  keV,  e^  is  the  average  electron-hole  pair  generation 
energy  which  is  3.65  ev/pair  for  silicon.  The  diode  gain  (G)  or 
number  of  holes  collected  per  primary  is  given  by: 

e.-jfd^E)  (6) 

til 

For  A(q)  the  4 order  polynomial  determined  by  Everhart  and  Hoff 
is  used: 

A (q)  = 0.6  + 6 . 2lq  - 12.40q2  + 5.69q3  (7) 

The  integration  is  terminated  at  z = 1.1  R(E)  where  the 
polynomial  approximation  to  A breaks  down.  Figures  64  and  65 
show  the  results  of  the  integration  of  equation  (6) . The  data 
shown  in  Figure  50  has  been  plotted  against  this  model.  Note  that 
the  fit  is  not  very  good.  This  is  typical  of  most  of  the  targets 
discussed  in  this  section,  and  is  not  surprising.  Much  of  the 
experimental  data  has  been  plotted  against  this  model  because  it 
gives  some  physical  insight  into  the  thickness  of  the  dead  like 
target  layer.  For  example  in  Figure  50  above  about  3 keV  the  data 
reasonably  well  follows  the  0.04  micron  dead  layer  curve.  This 
indicates  that  the  thickness  of  the  damaged  region  is  about  0.04 
microns.  The  deviations  at  smaller  energy  indicate  that  some  but 
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not  all  the  carriers  generated  below  0.04  microns  and  even  down  to 
0.01  microns  are  collected  by  the  diode. 

2.  Derivation  of  the  Collection  Probability  from  gain  versus  energy 
data 


The  first  order  model  can  be  improved  by  generalizing  the 
concept  of  a collection  probability.  The  target  gain  can  be 


represented  by: 


0° 

<5  r Jjz  P(Z)  £ (*,  e) 


where  g(z,E)  is  the  carrier  generation  function  in  equation  (4),  and 
P(z)  is  the  collection  probability,  for  carriers  generated  at  a depth 
z.  The  first  order  model  (equation  (4))  is  a special  case  of 
equation  (8)  with  P(z)  = 1 for  z>TQ  and  0 for  z<TD.  Calculation 
of  P(z)  from  models  of  the  doping  profiles  and  recombination  rates 
will  be  discussed  in  the  following  section.  In  this  section  the 
problem  of  calculating  P(z)  from  measurements  of  gain  versus  energy 
will  be  addressed.  A least  squares  method  will  be  used.  Equation  (8) 


is  rewritten  as: 


h«0) 


where  h = 1.1  so  that  the  integral  is  cut  off  at  the  depth  where 
g (z,E)  becomes  small.  A (a)  and  P(z)  are  written  as: 


(11) 


L,  i 

mbtu.  . I 


/7 


P(*)r  ^ 


K 


where  the  coefficients  b^  are  given  in  equation  (7).  Using  this 


notation  equation  (9)  becomes: 


\ - jt  K~i  r 1 


J"6 


(12) 


Performing  the  integrations  and  interchanging  summations  gives: 


&(£) 


e)  ' /L  a* C*  ^ 


K*o 


where 


(13) 


f'  ^ L-  L 


ck(e)-  ('-o©)CR^rs  &> 


K-0, 


4 *o 


(14) 


The  experimental  data  is  represented  by  a set  of  measurements  of 

gain  and  energy  and  weight  factors  as  0^,  E^,  W^)  i = 1,2, IM. 

W.  are  normally  assigned  according  to  the  degree  of  fit  desired  at 


energy  E^.  A logical  choice  is  = If ' % where  <7^  is  the  estimated 


error  or  standard  deviation  of  data  point  G^.  The  least  squares 
criteria  is: 


JM 


^[Gt-GCeS pW£  =■ 


l- 1 


(15) 
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Differentiating  with  respect  to  the  N coefficients  in  the  stan- 
dard fashion  and  assuming  a unique  minimum  exists,  yields  N 


simultaneous  equations  in  the  N unknowns  a (r  » 1,2,....N).: 
N r 


where 


J/A 

grK 

l-f 

ZM 

Dr  --  ^ <S;  Cr  CEi ) W; 
i=  | 


(17) 

(18) 


A computer  program  written  in  ALGOL  for  performing  this  least 
squares  fit  and  typical  output  is  shown  in  Figure  66.  This  program 
calculates  the  best  fit  P(z)  and  then  recalculates  the  gain  versus 
energy  curve  from  the  best  fife  P(z)  using  direct  integration  of 
equation  (8  ) and  summation  of  equation  (13).  The  simultaneous 
equation  set  (13)  is  solved  by  Crout  Reduction  (Ref.  11).  Typical 
results  of  this  program  are  shown  in  Figures  67  and  68.  These 
figures  show  the  effect  of  different  orders  of  fit  N.  The  order  of 
fit  must  be  chosen  carefully,  and  be  compatible  with  the  number  of 
data  points.  An  order  of  fit  comparable  to  the  number  of  data 
points  (e.g.  N = 6 in  this  case)  will  result  in  large  oscillations 
in  P(z)  because  it  is  possible  to  achieve  almost  an  exact  fit  to 
the  points  while  allowing  large  deviations  between  data  points  from 
the  smooth  curve  that  is  physically  expected.  This  is  particularly 
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uOOIO  BEGIN  COMMENT  CALCULATES  THE  COLLECTION  PROBABILITY  FROM 

00020  MEASURED  SAIN  VERSUS  ENERGY • USES  MODIFIED  E AND  H 

UG030  THIRD  ORDER  FIT  TO  GRUN  FUNCTION  AND  LEAST  SQUARES  CRITERIA; 

00G40  REAL  H, £u, SUM/ DET, EMAX/ Z, DEPI / DEPTH/ DEPTHM/ RR/ ENI / EMAXP/ ENERGY; 
00050  INTEGER  JM/  I / OPT/ N/ R/ K/  Q/  IM/* 

OOGoO  REAL  ARRAY  BLC  0 s 3 1 , D/  AC  0 : 1 0 3 / LC  1 : 1 00/  1 : 1 00 3 / V/ X C 1:1  0 3 /* 

00070  REAL  ARRAY  CACO: 100/ l:100]/BC0:10/0: 10  3/G/E/W/ ERROR! 1 1 1001; 

OOOoO  SWITCH  S1-b1/L2/L3/L4;  SWITCH  S2-LL1  / LL2/ LL3/ LL4I 
00062  SWITCH  S3-LX1/LX2; 

000*0  STRING  FILEl/DESC; 

00100  REAl.  PROCEDURE  PROSUi-iCI/M/N/P); 

oouo  value  M/N;  integer  i/M/N;  real  p; 

00120  BEGIN 

00130  REAL  tsum;  tsum-o; 

00140  FOR  I«M  STEP  1 UNTIL  N DO 
00150  TSUM-TSUM  *Pi 

00100  prosuh-tsum; 

00170  END  PROC  PROSUH/* 

00160  REAL  PROCEDURE  I NTEGRAL CF/ X/X1/X2/N); 

00190  ! INTEGRATES  BY  WEDDLES  RULE/  N MUST  BE  MOD  6! 

00200  VALUE  X1/X2/N; 

00210  INTEGER  l\i; 

00220  REAl.  X/X1/X2/F; 

00230  BEGIN 

00240  REAL,  SUM/DEE/ DEL2/XI; 

GGeoO  INTEGER  It 

002o0  IF  M 0 D ( M / 6 ) / s 0 THEN  PRINT!”  INTEGRATION  NOT  MOD  6”>; 
u027  0 DEL-(X2-X1 )/N;  DELa-E^OEL; 

ou2oo  xi»xi ; sum-o; 

00290  AGAIN:  FOR  1-0  STEP  2 UNTIL  6 DO 
00300  BEGIN  X-I*DEL+Xi;  SUM-SUM+F;  END! 

003 1 o x-del+xi;  sum-sum+s*f; 

00320  X-X+DEL2;  SUM-SUM+6*F; 

00330  X-X+UEL2;  sum-sum+s*f; 

00340  XI-b*DEL+Xi; 

003S0  IF  XI <X2-DEL  THEN  GO  TO  AGAIN; 

00360  INTEGRAL-SUM«‘UEL*>0.3z’ 

00370  END  WEDDLE  INTEGRAL; 

00360  PROCEDURE  SOLVECROUT UM ) ; VALUE  N;  INTEGER  N/' 

00390  BEGIN  COMMENT  NEEDS  MATRIX  L/  COEFF  VECTOR  V AND  SOLUTION  VECTOR  X 
00400  DECLARED  EXTERNAL  TO  PROCEDURE; 

00410  INT EG ER  US/l/U/US SI 

00420  REAL  ARRAY  ALC1:N,1:N3/  AVCl:NJ; 

00430  ALC  1/ n-LC  1/ n;  ALC2/  1I-LC2/  1 11 
00440  ALC  1/21-LC  1 /2I/ALC  1/  1 I/* 

00450  FOR  OS-1  STEP  1 UNTIi.  n do 
00460  BEGIN 

u047u  OSS-IF  US/*1  THeN  US  ELSE  31 
u04o0  FOR  I-USS  STEP  1 UNTIL  N DO 

00490  ALC  1/ JSJ-bC  I/USJ-PR0SU1-UQ/  1 / US-  1 / ALC  I / U3  *ALC  Q/ US  J ) /* 

00500  JSS-IF  OS/  = l THEN  US  ELSE  21 
uOS 1 0 FOR  D-oSS* 1 STEP  1 UNTIL  N DO 

ou520  begin  i-os; 

0u530  AuC 1/ JJ-!Lt I/DJ-PR05UM(U/ 1/ 1- 1/ ALC I / Q ] *ALC Q/ J J > > /ALC I / I J 1 

00540  end; 

ocsso  end; 


Figure  66a.  Computer  program  for  calculation  of  best  fit  collection 
probability  from  gain  versus  energy  data. 
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UuaeO  AVC 1 3-VC 1 3/ALC 1/ 1 37 

UOS  7 «j  FO H 1-2  jTZP  l UNTIL  N DO 

00a60  AVC  I 3-(.VC  I J-PROSUM! «/  1 / I * 1 > ALC  Ij  Ql+AVCQ]  ) ) /ALC  1/  I 37 
005*0  XCW3-AVCN37 

00600  FOR  1-N-l  STEP  -1  UNTIL  1 00 

006  10  XC  1 3 -AO  C 1 3 - PRO  SUN  (0/  I + 1 / N/  ALC  I / 03  *XC  Ci)  ) 7 
00620  END  PROCEDURE  CROUT; 

00o30  REAL  PROCEDURE  RANGECE37  REAL  £7 
i/OoAO  RANGE- 1 . 1 * ! 0 • 1 *Z) t 1 .7 57 

00o50  REAL  PROCEDURE  CCK/E37  REAL  £7  INTEGER  Ki 
OOooO  oEGIw 

00670  INTEJER  UJ  RtAi.  SUM 7 SUN-07 

Ouo60  FOR  u-0  SToP  1 UNTIL  UN  DO 

u0690  SUm-SUn+DdC U 3*ri*<J-fK*l  )/lii*0+\)} 

00700  C-E*RAW3£Cc.)  *R*SUn/*G7 

007  10  lifvD  PROC  CCH/E37 

00720  REAL  PROCEDURE  PR0b!£)7  REAL  i 7 
007  30  PROo-PROSUMCG/  0/ N/  t,C  03  *2.  t 03  7 

00731  REAL  PROCEDURE  3h1NIivT!E)7  RtAL  E7 
o07  32  DEO  IN  REAL  I IL  RAN7  RAW-RANGc! ENERGY) 7 

007  33  1 N- 1 NTE3RAL  ( PROSUM!  0/0/  UM/  dLC  03 1> < Z/RAN ) « 0)  fPROBC  Z ) / t,  0/H*RAN/  36  3 

0073A  GAININT-IN*£/!EG*RAN)7 

00735  END  PRO  GAIN  faY  INTE3RAT low! 

00736  REAL  PROCEDURE  GAINSUrKE)  7 REAL  E7 

00737  GAINS UN -PRO SUN  (U/0/N/ACQ3*C!0/£) ) 7 
007 AO  UN-37 

00750  H- 1.057 
007o0  EG-3.o5E-37 

00765  1-07  FOR  DLC I 3 -0 .6/ 6 .2  I / - 1 2 . A/ 5 . 69  DO  1-1+17 

00766  LAJLL3:  PRINT! "OPTI ONS s 1-TTY  , 2-FILE  INPUT")! 

00771  IwPUT ! OPT ) 7 GO  TO  S3COPT37 

U0772  uX2 : PRINT ("F ILE  NAN£"37  INPUTf < " ! A) "/ FILE 1 ) 7 
007  7 3 READCF  ILi  I * mCa50>  ",  DEF/  DESC3  7 
007  7a  REhDCF  IL£l/,,tI37">DEF/  11)3  7 
uG7  76  FOR  1-1  STEP  1 UNTIL  IN  DO 

0U776  READ!FIL£1/ "! 3Fo .2)"/  uc.F  , LC I 3/GC I 3/  ERRORC  13)7 
00777  RE'J 1NDCF ILE  1 ) 7 
00776  GO  TO  La! 

00760  LX  1 : F I LE 1 -"TTY" 7 ULSC-"  "7 

Gu7*0  PHI wT ! " I NPUT  ENERGY  IN  REV / G A I N/ ERROR  - EN£R3Y=0  TERMINATES")! 

00600  FOR  IN-1  STEP  1 UNTIL  1000  DO 

00601  bEG  I N I NPUTCEC  IN  3 / 3 C 1M 3 / ERRORC  IN  3 ) 7 

00602  IF  EC  I h 3 = 0 THEN  GO  TO  L57 

00603  END! 

00620  OEF:  PRINT ! "READI NG  ERROR  OF  FILE")!  GO  TO  1.L37 

00920  L5:  PRINT! "OPTIONS:  1-PRINT  DATA/  2 -SAVE  DATA/  3-RUN/  A - NEW  DATA")! 
00930  I.'jPUTCOPT)!  GO  TO  S I COPT  3 7 

00*36  Ll:  PRINT  ("F  li.£  NAME*"/ F ILE  1 ) 7 PRIhT!  D£SC)  7 

GO*AO  PRINT!"  ENERGY  GAIN  ERROR")!  PRINTC"  REV") 7 

00960  FOR  1-1  STEP  1 UNTIL  In  DO  PRINTr ! "! 3F6 . 2 7 "/ £C I 3/ G C I 3 / ERRORC I 3 ) 7 

00*ou  GO  TO  La! 

uO*  70  u2S  PR I wT ! "F ILE  NAME”)! 

U096D  INPUTF  !"!«o)"/FILiil  )7 

UD*o2  PRINT!  "FIi.ii  DESCRIPTION")!  INPUTF  !"  !A)"/D£SC)7 
uo  * 6 A 7/UTc,!r  IL*il/"!H»0)"  / DESd ) 7 
ou**u  JR1 TE!F IdEI / " l 1 3)"/ IN ) 7 
jluuO  F on  i-I  oTt,P  1 UNTIL  In  DO 

UiOlo  > n I T — ! r 1 ..i.  I /"!  3F  e .c  )"/  EC  1 3 / G C I 3 / dRRORC  I 3 ) 7 

ulula  ENUr  ILE!F  li.El  )7  HEW  1 No!  F ILo  1 3 7 
OlOcO  GO  To  067 


Figure  66a  (cont.) 
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L3:  PRINTC'PULYNOMIAL  ORDER  FUR  FIT  TO  PROoCE)"); 

INPUTCN); 

FOR  1-1  STEP  1 UNTIL  IM  DO  C I 3 - 1 E-4*<  ERRURC  I]  » ( -2  ) ) i 
FOR  K-0  STEP  1 UNTIL  N DO  FOR  I - 1 STEP  1 UNTIL  IM  DO 
CaCK, I 3-C<K,EC I ] ) j 
FOR  R-0  STEP  1 UNTIL  N DO 
FOR  K-R  STEP  1 UNTIL  N DO 

BCKj R3-BCR,K3-PROSUM<a, 1 , IM, CAC R, Q 3 *CAC K, Q 3 mV C Q3 ) i 

for  R-o  step  i Until  n do 

DC  hi -PROSUM (Q, 1 , im,gcu3*cacr, q3*jcq3 ); 

COMMENT  NOW  READY  TO  SOLVE  FOR  COEFF.  VECTOR  A i 
FOR  R-0  STEP  1 UNTIL  N DO 
BEGIN  VCR+13-DCR3; 

FOR  R-0  STEP  . 1 UNTIb  N Do  bC  R+  1 , K+  1 3 -bC  R,  K 3 i 
end;  go  to  sob; 

PRINT ("CACK* 1 1 K I");  FOR  K-0  STEF  1 UNTIL  N DO  FOR  1-1 
1 UNTIL  IM  DO  PKINTF<"<E13.3,2I2)",CACK,I3,K>I); 

PHI NTl "DC K> R3  «\  R"  ) ; FOR  R-0  STEP  1 UNTIL  N DO  FOR  R-0 
1 UNTIb  N DO  PRINTF < " < El  3 • 3> 2 1 2 ) "» BC R> R ] > K* h) i 
FOR  R-0  STbP  1 UNT I u N DO  PRINT C DC R3 , R) J 
SOLs  SOLVECROUTCN-l )i 

FOR  R-0  STEP  1 UNTIL  N DO  AC  R3 -XC  R-*- 1 3 1 
'.DISPLAY  HESULTS! 

PRINT  ( "COEFF  I C I ENTS  ARE"); 

PRINT ( " ACR3  R"); 

FOR  R-0  STEP  1 UNTIL  N DO 

FRINTF("('E14.7,  1 2 ) ",  AC  R3  , R)  ; 

EM AX -EC  1 3; 

FOR  1-2  STEP  1 UNTIL  IM  DO 

IF  £MAX<ECI3  THEN  EMAX-EC 13; 

LL5 : PR I NT  < "OPT  I ONS  s 1-COLL  PROEs, 2-GAIN  VS  E>  3-NEJ  DATA*  4- 

inputcopti;  go  to  s2Copt3; 

LLls  PRINT < "DEPTH  INCREMENT  AND  DEFTri  IN  MICRONS"); 

inputcdepi* depth) ; 

PRInTC "COBLECTION  DEPTH"); 

PRIimT  ( "PROD  MICRONS"); 

DEPT  HM -H  + RANGE ( EKAX ) i 
IF  DEPTH<DEPTHM  THEN  DEPTHM-DEPTH; 

FOR  dEPTH-0  STEP  DEPI  UNTIL  DEPTHM  DO 

PRINTF  <"(Fv  .5>F9.3)"*  PRO DC  DEPTH )* DEPTH ) ; 

GO  To  LLS; 

LL2 : PRINT ("ENERGY  INCREMENT  AND  MAX  ENERGY"); 

I NPU T < EN I * EMAXP ) i 

IF  EMAXP >£MAX  THEN  EMAXP-EMAX; 

PRINT ( "GAI N- I GAIN-S  ENERGY-KEV") ; 

FOR  ENERGY-ENI  STEP  ENI  UNTIL  EMAXP  DO 

PRI NTF  ( " < 3F9 • 2 ) "*  GAI NINTC  ENERG Y) * GAI WSUM ( ENERGY) » ENERG Y) 
GO  TO  LL5; 

LL4 : END 
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STEP 


QUIT") 


Figure  66a  (cont . ) 
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Figure  66b.  Typical  output  of  program  showing  data  access  from  a stored 
file. 
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COLLECTION  PROBABILITY 


a problem  at  the  larger  depths  z where  the  gain  versus  energy  curve 
is  insensitive  to  the  exact  shape  of  the  P(z)  curve.  Since  the 
generation  function  is  a very  broad  smooth  function,  the  gain  at  an 
energy  E is  a weighted  average  of  P(z)  over  all  z less  than  R(E) . 

For  small  z the  fits  obtained  are  essentially  independent  of  the 
order  N as  can  be  seen  in  Figure  67.  Collection  probabilities 
greater  than  one  are  not  physical  unless  avalanche  is  occurring 
in  the  diode  pn  junction  which  is  not  possible  for  the  low  diode 
bias  of  1.5  volts  used  for  these  measurements.  Hence  the  fit  for 
N = 3 or  4 is  the  most  reasonable  in  this  case.  The  solid  points 
shown  in  the  gain  versus  energy  curve  are  actual  measured  points. 

To  improve  the  fit,  extra  interpolated  points  indicated  by  crosses 
are  sometimes  added.  Figure  69  and  70  demonstrate  the  sensitivity 
of  the  fit  to  the  shape  of  P(z)  at  small  z.  Figure  69  shows  the 
collection  probability  from  Figure  67  for  the  3rd  order  fit 
and  the  dotted  line  is  an  assumed  collection  probability  which  is 
zero  for  z<0.04  microns.  This  corresponds  to  the  first  order  dead 
layer  model  with  a 0.04  micron  dead  layer.  Figure  70  shows  the  gain 
versus  energy  curves  determined  by  these  two  collection  probability 
functions  by  integration  of  equation  (9).  This  clearly  demonstrates 
that  the  shape  of  the  P(z)  curve  in  this  small  z region  is  not  an 
artifact  of  the  least  squares  fit.  The  finite  intercept  at  z = 0 
for  the  N = 3 fit  is  probably  an  artifact  since  there  are  no  data 
below  0.5  keV  to  determine  the  shape  of  P(z)  in  this  region.  Since 
the  derived  gain  for  the  N = 3 fit  deviates  from  the  assumed  data 
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DEPTH  Z (MICRONS) 

Figure  69.  Collection  Probability  P(z)  from  Figure 
67  for  N = 3 fit  compared  to  Tp  * 0.04u 

dead  layer  approximation  (dotted  curve) . 
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Figure  70.  Gain  vs.  energy  curves  derived  by  direct 
integration  of  equation  (9)  for  the  P(z) 
curves  in  Figure  69. 


point  at  .5  keV  significantly  the  true  P(z)  curve  probably  goes  to 
0 at  z = 0.  However,  as  will  be  seen  in  the  following  section  a 
finite  intercept  is  possible  if  the  surface  layer  of  high  recombin- 
ation is  of  near  zero  thickness  and  the  motion  of  generated  carriers 
is  controlled  only  by  diffusion. 

3.  Carrier  Diffusion  in  Planar  Diodes  with  a Doping  Gradient 
Guildberg  and  Schroder  (Ref.  3)  have  proposed  a model  for  the 
diffusion  of  carriers  in  the  presence  of  a field  due  to  a doping 
gradient  with  a dead  layer  at  one  surface.  The  model  is  shown  in 
Figure  71.  The  dead  layer  is  assumed  to  have  such  a high  electron 
hole  recombination  rate  that  no  carriers  generated  in  this  region 
escape  and  contribute  to  the  diode  signal. 

The  doping  gradient  extends  from  z = 0 to  z = S.  In  the  dead 
layer,  any  doping  is  ignored,  but  in  the  field  region  the  gradient 
is  assumed  to  be  such  that  the  electric  field,  e,  due  to  the  gradient, 
is  a constant.  In  the  bulk  region  the  doping  gradient  is  assumed  to 
be  so  small  that  the  field  is  approximately  zero.  This  region  may 
not  exist,  expecially  in  archival  targets.  The  space  charge  region 
is  associated  with  the  diode  depletion  field.  In  the  field  region 
and  bulk  region,  the  excess  carrier  density,  6p,  is  governed  by 
the  continuity  equation: 


DV*8p-v"£*  ZSp-SP/'r+S  * 0 


(19) 
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Figure  71.  Model  for  Diffusion  of  Carriers  in  Presence  of  Field 
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where  D is  the  minority  carrier  diffusion  constant,  u is  the  mobility, 
t is  the  recombination  time,  and  g is  the  generation  rate  of  the 
carriers.  It  is  assumed  that  t = « in  the  bulk  and  field  regions. 

The  boundary  condition  at  z = 6 is: 

~ $SP  (20) 

where  h is  the  unit  vector  normal  to  the  surface  at  z = 6 and  S is 
the  surface  recombination  velocity.  This  is  the  standard  boundary 
condition  for  an  idealized  surface  with  recombination  centers  (inter- 
face states).  The  boundary  condition  at  the  depletion  region  boundary 
is : 

sp(^  = w)=  0 

(21) 

because  the  very  high  fields  that  exist  in  this  region  sweep  out 
carriers  very  rapidly.  General  solutions  to  Equation  19  in  the  field 
region  are  of  the  form  Cq  + exp(Az)  and  in  the  field  free  region 
az  + 0.  The  additional  boundary  conditions  are  carrier  and  current 
continuity  across  boundaries  where  the  minority  current  is: 


j = -dS  8f 


(22) 


Using  Equation  19  and  these  boundary  conditions,  the  solutions  in 
each  region  can  be  found.  In  fact,  any  series  of  field  regions,  each 
with  different  uniform  fields,  can  be  solved  by  this  method,  which 
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reduces  to  the  solution  of  simultaneous  linear  equations.  The 
electric  field  due  to  a doping  gradient  is: 


(23) 


where  N = ND  - NA>  This  assumes  that  charge  neutrality  exists,  which 
is  similarly  required  for  Equation  19. 

Because  N(x)  <*  exp  (-dqp/kT),  where  cp  is  the  electron  potential, 

voltage  drop  between  Z^  and  Z ^ is 

Ai-'/f  /?*•  & f*  ^ •*  ($w) 


xi 


The  model  can  be  solved  if  the  carrier  generation  g = 

6 (z  - z^);  that  is  all  carriers  are  generated  in  a plane  at  a depth 
Zq.  A specific  case  solution  of  the  model  is  for  Z^  < 8^.  This  is 
of  most  interest  for  the  surface  diode  target,  because  the  electron 
beam  penetration  will  be  shallow  and  will  generate  carriers  only 
within  the  field  region.  The  solution  for  the  diode  gain  (signal 
current/beam  current)  normalized  to  the  maximum  possible  gain,  G^, 
which  occurs  for  S =0  is : 

-§•  6 " e***)  ft 


For  the  special  case  W 


6^  this  reduces  to: 


(24) 


(25) 
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The  constant  A = (q/kT)  e.  is  proportional  to  the  electric  field. 
Because  of  the  presence  of  at  least  a thin  damaged  surface  layer, 

S is  very  large.  The  limiting  maximum  value  for  S is  the  thermal 
velocity,  ~10^  cm/sec,  which  is  probably  approached  in  most  cases. 

To  calculate  the  normalized  gain  for  an  arbitrary  generation  function, 
g(z),  Equation  24  or  25  can  be  integrated  over  g(z)  i.e. : 

rSi  jx  $ 

G fw-apj  ~ RSi 

6m  (%(z0  + 6') 

J r (26) 

T-S, 


assuming  g (zq  >w  + s^)  =0. 

Figure  72  shows  the  dependence  of  the  normalized  gain  or  collec- 
tion probability  on  A or  the  electric  field  for  several  values  of 
surface  recombination  velocity.  This  figure  is  for  zq  = 0 and  is 
valid  if  AW>3  or  for  A<10^  cm  ^ and  W>0.03  microns,  which  covers  the 
range  of  values  of  A and  W of  interest. 

Figure  73  shows  the  dependence  of  normalized  gain  or  collection 
probability  on  generation  depth  zq . This  figure  is  for  S/D  = 10 
cm  ^ and  again  is  valid  for  AW>3.  For  AW<3  it  is  specific  to  W = 

0.5  microns.  Note  the  condition  AW<3  includes  all  the  curves  for 
A<0 . 

To  apply  this  model  to  prediction  of  collection  probabilities 
and  gain  versus  energy  curves  it  is  necessary  to  know  the  doping 
profiles  and  recombination  characteristics  of  the  targets.  As 
discussed  in  a following  section,  the  second  quarter  will  include 
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work  directed  towards  this  end.  However,  based  on  the  expected 
implantation  profiles  coupled  with  some  intuition  it  is  possible 
to  make  some  predictions  which  agree  qualitatively  with  the  exper- 
imental results. 

Figure  74  shows  the  model  to  be  used.  Figure  74a  is  the 
expected  ion  implantation  profile  which  in  LSS  theory  is  a simple 
gaussian.  Assuming  a gaussian  profile: 


and 


r 2f/ Cfff  ARr) 


Njj  • Uf 


(27) 


(28) 


where  F is  the  implantation  fluence  (and  A.Rp  the  standard  deviation 
of  the  expected  range  (see  Appendix  I).  The  position  of  the  center 
profile  is 


z = R - d 
p p ox 


(29) 


where  R^  is  the  projected  range  (see  Appendix  I)  and  dQx  is  the 
thickness  of  the  oxide  layer  present  during  implant.  For  the  low 
mass  ions  used  in  these  studies  so  far  the  stopping  powers  of  SiC>2 
and  Si  are  identical.  The  implantation  tail  shown  will  have  to  be 
guessed  at  for  now  until  the  profiling  measurements  become  available. 
Wj  is  the  ideal  junction  depth  which  is  determined  by: 


l^j  — hi  A -0s  Nptfif  (r  /t±  *r  ) ~ ^3 


(30) 
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Figure  74.  Model  Approximation  of:  (a)  implanted  doping  profile 
(b)  Electric  field  distribution. 
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COLLECTION  PROBABILITY 


Figur-  74c. 


Model  Approximation  of  Collection  Probability. 


where  Ng  is  the  acceptor  doping  of  the  starting  wafer.  For  most 

~ 15  -3 

of  the  samples  used  so  far  Ng  =10  cm  The  effect  of  tailing 

can  be  quite  severe  (see  Figure  22),  and  we  will  estimate  that 

W.  = WT  + 2(W_  - z ). 

A I I p 

Assuming  that  the  doping  profiles  profile  is  given  by  equation 
27  and  28  at  least  in  the  vicinity  of  the  doping  peak.  The  electric 
field  is  given  by  equation  23  is: 

r~  z JS.  . £ Cz  ~ 

% (31) 


In  order  to  fit  this  field  distribution  into  the  solution  represented 
by  equation  (25)  the  dotted  approximation  shown  in  Figure  74b  will 
be  used  with: 


ET  (z^Zf  =- 


(32) 


Equation  24  could  also  be  used  and  this  would  allow  a field  free 
region  near  the  junction.  In  general,  equation  19  can  be  solved 
for  any  number  of  uniform  field  regions  to  better  approximate  the 
actual  field  distribution.  It  is  also  possible  to  solve  equation 
19  by  finite  difference  relaxation  for  any  field  and  recombination 
distribution.  This  will  be  considered  during  the  second  quarter  if 
it  seems  desirable. 

Using  figures  72,  73,  and  74  predictions  about  the  collection 
probability  can  be  made.  Table  4 shows  the  calculated  values  of 
A^  = (q/kT)  E^,  Zp,  etc.  for  several  fabricated  planar  diodes 


Table  4 

CALCULATED  VALUES  FOR  ARSENIC  AND  PHOSPHORUS  PLANAR  DIODES 

Energy  d R AR  F Z En  A„  W,  W, 

^ ox  p p p R R 1 A 

-3  -1 

Sample  Species  keV  cm  v/cm  cm 

A7B  As  30  — 210  47  1014  210  l.lxlO5  4.2xl06  375  700 

A8A  As  100  — 574  122  10 14  574  4.1xl04  1.6xl05  986  1810 

P12C  P 100  — 1230  350  1014  1230  1.4xl04  5.7xl05  2350  4600 

P11C  P 100  1400  1230  350  1014  170  1.4xl04  5.7xl05  950  3200 
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discussed  in  this  section  (see  Tables  I & II) . Note  that  these  are 
calculated  only  from  the  implant  conditions  and  not  from  measured 
doping  profiles.  In  some  cases  the  values  for  are  extremely  high 
corresponding  to  electric  fields  ~105  v/cm.  Such  high  fields  are 
not  likely  to  actually  exist  in  these  samples  because  the  doping 
profile  is  not  as  sharp  as  assumed,  equation  23  does  not  apply  to 
very  steep  doping  gradients  where  the  Debye  length  X = 
is  long  compared  to  the  distance  over  which  the  doping  ischanging. 

For  ND  5=  10^  cm~^,  ^ = 100  8.  Hence  we  will  assume  A = 10**  cm  ^ 
is  a more  realistic  estimate  of  the  built-in  fields.  From  Figure  73 
note  that  for  z<z^  (Region  I)  A = -10"*  cm  ^ effectively  constrains 
all  the  carriers  within  that  region  and  few  escape  past  z = z^. 

That  is,  this  region  acts  as  a potential  well  which  holds  the  carriers 
until  they  recombine  at  the  surface  z = 0 which  has  high  surface 
recombination.  In  applying  Figure  73  to  region  I it  can  be  assumed 
that  region  II  (z>z^)  acts  as  a collecting  diode  (i.e.  zp  = W) . 

Hence  the  collection  probability  for  carriers  generated  in  Region 
I is  near  zero  as  indicated  in  Figure  74c.  Carriers  on  the  other 
side  of  Zp  will  have  a collection  probability  somewhat  like  the 
curve  A = 10^  cm”^  in  Figure  73.  For  carriers  generated  in  Region  II 
the  boundary  z = z^  is  a high  recombination  surface,  since  carriers 
which  pass  this  point  to  z<z^  are  effectively  lost  and  recombine 
as  Region  I.  The  shape  of  the  curve  in  Region  II  (Figure  74c)  does 
not  turn  up  as  sharply  as  in  figure  73  because  the  electric  field 
is  zero  ai  z and  increases  gradually  rather  than  having  a constant 
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large  value  from  z to  W.  as  assumed  in  the  calculation.  The  exact 

P A 


shape  of  the  curve  in  region  II  will  depend  on  the  electric  field 
distribution.  As  mentioned  earlier  it  is  possible  to  calculate 
exactly  the  collection  probability  for  any  doping  profile  and  this 
will  be  considered  during  the  following  quarters  if  it  seems  approp- 
riate. The  curve  shown  in  Figure  74c  is  a reasonable  approximation 
to  the  result  of  such  an  exact  calculation. 

Figure  75  shows  the  calculated  collection  probability  curves 
from  the  data  described  earlier  in  this  section  compared  to  the 
estimated  collection  probability  based  on  the  model  described  by 
Figure  74  and  equations  27  through  32.  As  can  be  seen  the  predic- 
tions are  reasonably  consistent.  In  particular  the  simple  model 
just  discussed  seems  to  predict  the  general  trends  of  the  probability 
curves.  This  indicates  that  the  basic  concepts  of  a dead  like  region 
due  to  a negative  electric  field  and  back  diffusion  of  carriers 
against  the  positive  electric  field  can  explain  most  of  the  dead 
layer  effects  observed.  It  is  possible  that  in  carefully  prepared 
samples  there  is  no  dead  layer  effect  in  the  sense  of  a high  recom- 
bination region.  Further  investigation  will  be  necessary  to  resolve 
this . 

For  reference  the  best  fit  collection  probability  curves  and 
gain  energy  curves  for  several  additional  samples  discussed  in  the 
planar  diode  fabrication  and  measurements  sections  are  included  in 
the  following  figures  76  through  80.  In  each  case  the  order  of  the 
fit  is  indicated.  Figure  78  again  compares  the  first  order  dead  layer 
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Figure  75 . Comparisons  of  calculated  collection  probability  by 

least  squares  fit  with  the  model  calculation  based  on 
figure  74  and  Table  IV.  (a)  Sample  A8A,  lOOkeV  As  at 
10*^  cm"2.  <100''  orientation  annealed  for  1 hour  in  Ar 
at  650°C.  (b)  Sample  A7B,  30  keV  As  at  10^  cm“^t  <111^ 

orientation,  annealed  for  1 hour  in  Ar  at  770°C. 
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Figure  75  (Cont.)  Comparisons  of  calculated  collection  probability 
by  least  squares  fit  with  the  model  calculation  based  on 
figure  74  and  Table  IV.  (c)  Sample  P12C,  100  keV  P at 
10™  cm  , <111>  orientation  annealed  for  1 hour  in  Ar  at 
760°C.  (d)  Sample  P11C,  100  keV  P at  10i4  cm-2.  <111  > 
orientation  annealed  for  1 hour  in  Ar  at  760°C. 
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Figure  76b  Gain  vs.  Energy  for  Sample  A8C,  lO1^  - As  - 100  keV 
<100>  after  760°C. 
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Figure  77a.  Best  Fit  Collection  Probability  for  Sample  A8A 
10 14  - As  - lOOkeV  after  650°C  anneal. 
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Figure  78a.  Best  Fit  Collection  Probability  for  Sample  P11C  - 
10 14  - P - lOOkeV  - <111>  through  1400  A oxide. 
Curve  labeled  = 0.06u  compare  first  order  model 
approximation  to  Data. 


14 

Figure  78b.  Gain  vs  Energy  for  Sample  P11C  - 10  - P-lOOkeV- 

<111>  through  1400  8 oxide.  Curve  labeled  TD  = 
compare  first  order  model  approximation  to  Data. 
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Figure  79a.  Best  Fit  Collection  Probability  for  Sample  P12C 
P - lOOkeV  - <111>  after  760°C  anneal. 
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Figure  80a. 


Best  Collection  Probability  for  two  implanted  sample 
annealed  at  different  temperatures.  Both  A7B  and  A7 
arj|  sections  of  the  same  wafer,  implanted  with  30  ke 
As  at  10^  cm~2  into  <111>  p type  silicon.  A7B  was 
annealed  at  760°C  and  A7C  at  650°C  for  1 hour  in  Ar. 
The  difference  for  z<0.02  microns  may  not  be  real  si 
there  is  no  gain  data  below  1 keV.  The  electron  ran 
at  1 keV  is  200  8.  The  differences  for  0.02  microns 
microns  are  real.  Comparison  of  these  samples  also 
some  idea  of  the  sensitivity  of  the  gain  curves  to  t 
collection  probability. 


BEAM  ENERGY  ( keV) 

Figure  80b.  The  Gain  versus  Energy  Data  used  for  the  fit  and  also  the 
best  fit  gain  curves  calculated.  Both  A^B  and  A7C  are  sections  of 
the  same  wafer,  implanted  with  30  keV  As  at  10^  cm"^  into  <111> 
p type  silicon.  A7B  was  annealed  at  760°C  and  A7C  at  650°C  for  1 hour 
in  Ar.  The  differences  for  z<0.02  microns  may  not  be  real  since  there 
is  no  gain  data  below  1 keV.  The  electron  range  at  1 keV  is  200  8. 

The  differences  for  0.02  microns<z<0. 2 microns  are  real.  Comparison 
of  these  samples  also  gives  some  idea  of  the  sensitivity  of  the  gain 
curves  to  the  collection  probability. 


model  to  the  best  fit  collection  probability  showing  again  that  the 
first  order  model  is  inadequate  to  fit  the  data.  Figure  80  compares 
the  collection  probability  for  two  As  implanted  samples  annealed 
at  different  temperatures.  Note  that  the  collection  probability 
appears  to  be  shifted  to  the  left  by  a distance  of  about  200  8. 

This  could  be  explained  by  a decrease  in  the  junction  depth  at  the 
higher  temperature  anneal.  This  could  occur  if  the  amount  of  active 
dopant  decreased  at  the  higher  temperature  anneal.  This  type  of 
behavior  can  occur  for  Boron  implants  but  does  not  usually  occur 
for  As.  Another  explanation  is  that  the  amount  of  active  dopant 
increased  only  in  a small  region  near  the  peak  thereby  increases 
the  built-in  field  in  this  region.  Note  that  if  the  amount  on 
active  dopant  increases  by  a constant  ratio  over  the  entire  n+  region 
then  the  built-in  field  does  not  increase  (see  equation  23).  Profil- 
ing of  these  samples  may  resolve  this  question. 

Determination  of  the  implanted  doping  profiles  is  planned  for 
next  quarter.  One  approach  will  use  sheet  resistivity  measurements 
coupled  with  controlled  anodic  oxidation  and  oxide  stripping.  (Ref. 

9).  This  technique  measures  the  electrically  active  dopant  or  carrier 
concentration  profile  assuming  the  mobility  is  constant  and  known. 

The  sheet  resistivity  will  be  determined  by  four  point  probe.  Calib- 
ration of  the  oxide  thickness  will  be  carried  out  by  ellipsometry 
(Ref.  10).  Computer  programs  for  all  the  necessary  calculations 
and  corrections  (Ref.  10)  have  been  acquired  and  developed. 


The  second  technique  will  be  secondary  ion  mass  spectrometry 
(SIMS) . This  measures  the  total  dopant  profile  (active  plus 
inactive)  and  is  complementary  to  the  resistivity  technique. 
Unless  a SIMS  unit  is  acquired  by  the  Center,  this  will  be  done 
on  a contract  basis  with  outside  laboratories. 
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Section  IV 


RESOLUTION  MODEL  I 


Figure  19  shows  a somewhat  idealized  model  of  the  surface  diode 
target.  This  can  be  even  more  idealized  as  shown  in  Figure  81. 

Radial  symmetry  is  implied.  The  areas  labeled  P+  are  the  written 
areas  and  should  be  thought  of  as  the  outer  edges  of  the  P+N+  diode 
depletion  field  in  the  case  of  surface  diode  writing  or  the  boundary 
of  the  high  recombination  region  in  the  case  of  damage  writing.  In 
either  case  the  boundary  condition  along  these  surfaces  is  V = 0 
where  V is  the  excess  minority  carrier  density. 

The  motion  of  generated  excess  carriers  is  governed  by  the 
continuity  equation  as  discussed  earlier  in  this  section. 

DVaV  = 0 

(33) 

where  D is  the  minority  carrier  diffusion  constant,  p.  is  the  mobility, 
r the  minority  carrier  recombination  time  and  g is  the  generation 
rate  of  carriers.  The  boundary  conditions  can  very  generally  be 
expressed  as: 


(34) 
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where  £ is  the  unit  vector  normal  to  the  surface  and  S is  the  sur- 
face recombination  velocity.  For  the  edge  of  the  written  regions 
S-*«  or  equivalently  V = 0 as  previously  mentioned.  The  boundary 
condition  at  the  z = 0 surface  for  r = a/2  is  assumed  to  be  S = 0. 

This  is  equivalent  to  no  carrier  recombinations  at  the  z = 0 surface, 
or  to  no  dead  layer.  The  goal  of  the  dead  layer  studies  is  to 
produce  targets  approximating  this  boundary  condition. 

Carriers  are  assumed  to  be  generated  in  the  region  z<d  and 
r<a/2  during  reading.  The  readout  gain  is  then  the  number  of  carriers 
collected  by  the  diode  (z  = W surface)  normalized  to  the  number  gen- 
erated. 

While  it  is  the  goal  to  eventually  develop  solutions  for  this 
general  resolution  problem  just  described,  a more  restricted  problem 
was  considered  during  the  first  quarter.  First,  the  effects  of 


electric  fields  and  carrier  recombination  are  ignored. 

Second,  the  geometry  shown  in  Figure  82  is  used.  This  can  be  thought 
of  as  two  extreme  limits  of  Figure  81  as  indicated  in  Figure  83.  The 


that  carriers  which  cross  the  dotted  boundary  in  Figure  83a  are 
assumed  to  be  collected  by  the  N+P  diode  (planar  diode)  and  do  not 
back  diffuse  into  the  z<P  region  or  recombine  at  the  z - p surface 
as  indicated  in  Figure  83a.  This  approximation  tends  to  give  a higher 
readout  gain  than  would  actually  be  observed.  If  T * W (Figure  83d) 
then  the  approximation  being  made  is  that  the  carriers  that  cross  the 
dotted  boundary  in  Figure  83c  do  not  back  diffuse  into  the  r<a/2 
region  or  to  the  collection  junction  at  z * W but  recombine  at  the 
z = P surface.  If  the  readout  gain  is  taken  as  the  carrier  flux  to 
the  z = T * W surface  in  Figure  83d  then  some  carriers  which  might 
have  reached  the  diode  will  not  be  included  and  the  readout  gain  will 
be  lower  for  the  Figure  83d  approximation  than  for  Figure  83c.  Here 
the  calculated  gain  for  T = P & T = W brackets  the  actual  gain  for 
the  model  of  Figure  81. 

With  the  assumptions  that  E * 0 and  r = 0 equation  (33)  reduces 
to  the  form  of  Poissons  equation.  This  is  th«  reason  for  the  intro- 
duction of  the  notation  V for  fiP  with  g equivalent  to  electric  charge 

density.  The  boundary  conditions  are  V - 0 for  r = a/2,  V = 0,  for 

z = T and  = 0 for  z = 0.  Hence  the  problem  is  equivalent  to  a 

conductive  cylindrical  box  with  an  open  end  where  the  boundary  condit- 
ion at  the  open  end  is  aV  = 0.  Since  closed  form  solutions  are 
available  for  charges  in  cylindrical  conductive  boxes,  the  method  of 
images  can  be  used  to  yield  the  desired  solution.  This  is  summarized 
in  Figure  84  along  with  the  notation  to  be  used,  q is  the  strength 
of  the  point  injection  of  carriers  at  points  z * d and  z =-d.  This 


3 3 

is  equivalent  to  g = q 6 (r-r)  where  5 (r)  is  the 

three  dimensional  delta  function.  The  net  rate  of  carrier  generation 
is  then  q.  By  symmetry  —■  * 0 at  z =0.  That  is,  the  net  carrier 
flux  across  the  z - 0 plane  is  zero. 

Solutions  of  the  potential  problem  shown  in  Figure  84b  can  be 
obtained  by  a series  expansion  in  Bessel  functions.  See  for  example 


Smythe  (Ref.  12)  p.  189-90.  The  solution  to  Figure  84b  is: 
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These  are  also  solutions  to  Figure  84a  since  they  satisfy  the 
boundary  conditions  and  are  solutions  of  equation  33.  Only  the 
solution  for  z>d  is  needed  for  the  following  calculations,  the 
-d<z<d  is  included  for  reference.  Note  that  A is  equivalent  to  a/2. 

The  normal  derivative  of  the  carrier  density  at  z = T is: 


a V | 
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Integrating  equation  38  according  to  equation  35  gives  the 
normalized  readout  gain  (with  q = 1) 
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Equation  39  can  be  easily  summed  using  a computer.  Typically  only 
5 or  10  terms  are  necessary  to  give  convergence  <10  unless 

■j 

T - d/T  < 0.1. 

Figure  85  shows  the  result  of  evaluation  of  equation  39.  As  an 
example  of  the  application  of  these  results  consider  the  model  of 
Figure  81  with  P/w  =0.5,  P = 0.05  microns  and  a generation  depth 
d = P/2.  Figure  86  shows  the  calculated  normalized  gain  as  a function 
of  2a  which  can  be  identified  as  the  bit  spacing  in  Figure  81  which 
is  an  isolated  zero  pattern.  The  actual  gain  for  the  model  of  Figure 
81  will  be  between  these  two  limiting  curves  as  discussed  earlier. 

Note  that  2a /W  should  be  larger  than  2 for  large  readout  gain.  For 
bit  spacings  ~0.1  micron  a junction  depth  W of  ^0.05  microns  is 
therefore  desired.  Inclusion  of  electric  fields  in  equation  33  will 
modify  this  result  and  should  increase  the  gain  for  a given  bit 
spacing  relaxing  the  requirement  on  junction  depth.  Inclusion  of 
these  electric  field  effects  will  be  considered  during  the  following 
quarters.  This  preliminary  model  will  then  be  useful  as  a limiting 
case  check  on  the  more  sophisticated  models. 
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Figure  85a. 


Normalized  readout  gain  form  equation  (39).  D is  the 
generation  depth  of  the  carriers,  and  T the  thickness 
of  the  T and  A the  dimensions  of  the  diffusion  region 
(see  Figure  82). 
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Figure  85b.  Normalized  readout  gain  form  equation  (39).  D is  the 
generation  depth  of  the  carriers,  and  T the  thickness 
of  the  T and  A the  dimensions  of  the  diffusion  region 
(see  Figure  82). 
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.ized  gain  from  equation  (39)  for  P/w  = 
0.25  (see  Figure  81  and  82)  for  the  li 
described  by  Figure  83  of  T = WJ  T = P 


Section  V 


ION  IMPLANTATION  WRITING 

Since  focused  ion  beam  systems  are  not  available,  all  ion  beam 
writing  experiments  during  the  first  phase  will  use  aperature  masks 
on  the  surface  of  planar  diodes  to  define  the  areas  written  by  a 
defocused  ion  beam.  This  also  has  the  advantage  of  providing  com- 
plete flexibility  in  choice  of  ion  species  and  energy.  Figure  87 
shows  schematically  the  experimental  method.  The  aperature  mask 
will  normally  be  a resist  (photo  resist  or  electron  beam  resist)  of 
thickness  :»0.2  microns  sufficient  to  stop  the  ion  beam.  With  care- 
fully prepared  and  exposed  electron  beam  resists,  it  will  be  possible 
to  implant  areas  0.1  microns  in  diameter  and  demonstrate  the  bit 
packing  capability  of  the  target.  For  other  writing  evaluation  tests 
photo  masks  will  be  used.  The  edges  of  carefully  prepared  photo 
masks  can  be  reasonably  sharp  (possibly  as  good  as  0.1  microns)  so 
that  something  about  the  bit  packing  density  can  still  be  learned  . 
The  photo  mask  technique  has  the  advantage  of  simplicity  for  prelim- 
inary evaluation  experiments. 

An  initial  writing  test  was  performed  prior  to  the  submission  of 
the  proposal.  These  early  tests  were  performed  using  40  ohm-cm, 

Boron  doped  <100>  silicon.  The  N+  surface  was  formed  using  900°C 

POCl^  gas  diffusion.  The  surface  doping  level  was  calculated  to  be 

19  -3  + 

~10  cm  . The  sheet  resistance  of  the  N layer  after  deposition 
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Figure  87.  Aperature  Mask  Ion  Beam  Writing 
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and  drive-in  was  100  ohms /square.  The  depth  of  the  N surface  is 
estimated  as  0.5  microns. 

A mask  was  formed  on  the  target  surface  using  electron  beam 
exposure  of  polymethyl  methacrylate  (FMM)  and  development  in  isopropyl 
alcohol  and  methyl  isobutyl  ketone.  This  is  a standard  electron 
beam  resist  system,  which  is  capable  of  submicron  resolution  (Refs. 

13  and  14).  An  optical  photograph  of  an  area  of  the  generated  mask 
is  shown  in  Figure  88.  The  resist  thickness  was  0.3  micron.  The 
mask  exposures  were  performed  in  the  SEM.  After  development  of  the 
resist  pattern,  the  test  wafer  was  exposed  with  20-keV  boron  ions. 

20  keV  boron  penetrates  0.07  micron  with  a gaussian  spread  in  the 

12 

penetration  depth  (ARp)  of  0.028  micron.  Implant  fluences  of  7 x 12 
13  2 

and  5 x 10  ions /cm  were  used.  This  should  produce  peak  acceptor 

18  18  3 

levels  of  1.4  x 10  and  10  boron/cm  , respectively,  if  all  of  the 
boron  is  active.  Very  little  difference  was  noted  in  the  behavior 
of  the  two  implant  fluences. 

Figures  89  and  90  show  diode  current  readout  signals  from  this 
target,  using  a 3 keV  electron  beam.  The  modulation  in  the  video 
pictures  is  the  diode  signal  current. 

Figure  91  shows  the  signal  levels  from  recorded  bit  and  bar  pat- 
terns of  various  dimensions  as  a function  of  the  reading  beam  current. 
Ones  are  the  signals  between  the  implanted  regions  in  Figure  89. 

Zeros  are  the  signals  from  the  implanted  region.  Note  that  the 
resolution  is  relatively  independent  of  beam  current  at  least  up  to 
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Figure  88. 
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Optical  Micrograph  Resist  Pattern  Used  as 
Ion  Mask. 
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Figure  89.  SEM  Micrograph  Showing  Implanted  Diode  Patterns 
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SKM  Micrograph  Showing 
Implanted  Diode  Pattern 
Diodes  on  1.3  Micron 
Centers . 


Readout  Waveform 
1.1  and  1.3  Micron 
Spaced  Implanted  Diodes 


Figure  90.  Implanted  Diode  Archival  Memory  Readout 

("The  largest  gain  modulation  is  100  times 
the  reading  beam  current.) 
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Figure  91.  Effect  of  Reading  Beam  Current  and  Bit  Spacing 
on  Readout  Signal  Levels  (The  upper  set  of  data 
points  (ones)  are  for  the  unimplanted  area 
between  the  implanted  spots.  The  low  sets  of 
points  (zeros)  are  for  implanted  spots  of  the 
indicated  size. ) 
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-0.5  nA.  Without  more  careful  additional  work,  it  is  difficult  to 
separate  the  effects  of  increasing  the  electron  beam  spot  size,  as 
the  beam  current  is  increased,  from  any  saturation  effects  on  the 
resolution,  but  all  of  the  resolution  degradation  can  be  attributed 
to  increasing  the  electron  beam  spot  size.  Figure  92  shows  the 
normalized  signal  levels  at  a low  reading  current,  where  the  spot 
size  is  smallest  (~0.5  micron),  versus  bit  size.  As  can  be  seen,  the 
recorded  patterns  are  well  resolved  down  to  at  least  0.6-micron  bit 
size.  An  actual  signal  trace  at  the  0.6-micron  bit  size  is  shown 
in  Figure  90.  The  zeros  (implanted  areas)  show  resolution  degrada- 
tion at  the  smaller  sizes,  as  would  be  expected  if  most  the  resolution 
degradation  is  due  to  electron  beam  spot  size  effects.  Because  of 
depletion  field  spreading,  implanted  areas  should  always  appear 
larger  than  expected  in  the  readout  signal  and  should  not  show  this 
gain  increase  at  smaller  sizes.  However,  if  the  read  electron  beam 
spot  size  is  too  large  during  reading,  then  some  ionization  is  always 
created  outside  the  implanted  zero,  and  the  effect  seen  in  these 
readout  signals  occurs.  An  alternate  explanation  for  the  resolution 
degradation  at  0.5  micron  is  the  lack  of  edge  sharpness  in  the  PMM 
mask.  This  would  produce  a tapered  implant  distribution  that  could 
also  explain  the  observed  result. 

It  has  been  tacitly  assumed  that  this  writing  evaluation  demon- 
strates the  surface  diode  type  of  writing  mechanism.  However  it  is 
possible  that  the  readout  observed  depended  only  upon  the  damage  done 
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Figure  92.  Dependence  of  Readout  Resolution  on  Written  Bit  Size 
Upper  Set  Ones  and  Lower  Set  Zeros.  (Degradation  of 
zero  resolution  indicates  the  reading  resolution 
degradation  is  due  to  electron  beam  spot  size  effects.) 


to  the  silicon  by  the  P ion  beam.  Writing  studies  were  not  planned 
until  the  second  quarter  of  the  contract  but  it  seems  desirable  to 
settle  this  question  by  an  exploratory  writing  experiment.  He+  ions 
were  chosen  as  the  writing  ions.  Being  an  inert  gas  with  no  known 
doping  effects  in  silicon  they  seem  appropriate  as  a test  of  the 
damage  writing  method.  A planar  diode  target  (P12B)  was  fabricated 

from  <111>,  12-18  ohm-cm  P-Si  by  implantation  of  P+  ions  at  100  keV 

o 2 

at  room  temperature,  7 off  axis  to  a fluence  of  10  cm  . The 

target  was  annealed  for  1 hour  in  Argon  at  760°C.  (This  target  is 

identical  to  P12C).  The  target  was  then  coated  with  resist  and 

patterned  with  a relatively  coarse  bar  mask  (Air  Force  Resolution 

Chart).  The  patterned  wafer  was  then  exposed  to  a 30  keV  room  temp- 

erature  He  implant,  7°  off  axis  to  a fluence  of  10  cm  . Figure 

93  is  a secondary  emission  image  from  the  SEM  of  a part  of  the  exposed 

resist.  Figure  94  and  95  are  SEM  micrographs  using  diode  signal  modul 

ation  after  stripping  of  the  resist.  Figure  96  shows  scope  traces  of 

line  scans  across  the  edge  of  one  of  the  bars.  The  edge  sharpness  of 

the  diode  signal  is  about  0.6  microns  at  5 keV  beam  energy  and  1.2 

microns  at  7 keV.  Secondary  emission  pictures  of  indicated  that  the 

spot  size  of  the  SEM  was  ~0.4  microns  at  5 keV  and  ~0.2  microns  at  7 

keV.  In  addition  the  broadening  of  the  electron  beam  by  scattering 

in  the  substrate  is  ~0.3  microns  at  5 keV  and  ~0.6  microns  at  7 keV 

(Figure  17).  The  total  width  of  the  traces  in  Figure  96  will  be 

approximately  twice  the  effective  spot  size  in  the  target  including 


Figure  94.  Scanning  Electron  Micrograph  using  Diode  Current 
Modulation  Readout  showing  He  Implanted  Archival 
Target  P12D  Electron  Beam  Energy  3keV. 
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SEM  Micrograph  using  Diode  Signal  Modulation  showing 
He-Implanted  Bar  Pattern  in  Target  P12B.  Electron 
Beam  Energy  3keV.  Spacing  of  Finest  Bar  Patterns 
is  38  Microns . 


Figure  96 
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Single  line  scope  traces  across  the  edge  of  one  of 
the  wide  bais  on  He  implanted  target  P12B.  Multiple 
traces  at  random  positions  along  the  bar  are  shown 
to  indicate  uniformity.  Top  flat  line  is  zero  signal 
or  beam  off  trace.  Scope  sensitivity  0.2  uA/cm  target 
current.  (a)  Beam  energy  7 kcV,  horizontal  distance 
scale  0.  >2  microns/div  as  indicated,  beam  current  ~6  x 
lO'^^A(b)  Beam  energy  5 keV.  horizontal  distance  scale 
0.45  mi  crons ''div  as  indicated.  Beam  current  ~1.5  x 10'^A 
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scattering  and  incident  spot  size.  This  gives  an  estimated  edge 
sharpness  comparable  to  that  observed  so  all  the  resolution  degrad- 
ation can  be  accounted  for  by  spot  size  and  scattering  effects. 
Readout  of  this  pattern  with  a smaller  spot  size  at  lower  beam  energy 
will  be  very  interesting  when  the  new  Coates  & Welter  SEM  under 
order  by  General  Electric  becomes  available  later  in  the  year. 

Figure  97  shows  the  measured  gain  versus  beam  energy  energy 
curves  and  best  fit  collection  probability  for  the  He  written  target 
areas,  unwritten  areas  and  the  control  sample.  Note  the  very  large 
gain  differences  between  the  written  and  unwritten  areas  and  the 
good  reproducibility  between  the  two  identically  processed  targets 
P12C  and  P12B.  Readout  at  2.5  keV  would  be  possible  with  this 
target  and  writing  method.  Low  voltage  readout  is  important  because 
it  minimizes  beam  scattering  during  readout.  Figure  97a  shows  the 
collection  probability  calculated  from  the  data  of  figure  97b.  The 
expected  region  of  He  damage  between  z = R^  + aR^  is  indicated.  The 
collection  probability  is  significantly  reduced  in  this  area  and  at 
smaller  z.  It  is  reduced  at  smaller  z probably  because  of  the  block- 
ing effect  of  the  damaged  region.  This  also  shows  that  as  expected 
the  writing  was  very  deep  in  this  experiment.  Clearly  0.1  micron 
resolution  will  be  difficult  to  achieve  with  this  energy  He  ions. 
Lower  energy  or  heavier  mass  ions  are  preferred. 

This  preliminary  test  indicates  that  damage  writing  is  possible 
in  the  planar  diode  target.  Damage  writing  depends  on  a reduction 
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Lgure  97.  Comparison  of  the  calculated  collection  probability  and 
Gain  versus  energy  for  the  He  written  sample  and  control. 
Wafers  P12C  and  P12B  were  implanted  at  100  keV  at  10^  P 
into  <111>  12-18  ohm-cm  n type  silicon  and  annealed  for 
1 hour  in  Ar  at  760°C.  P12B  was  written  in  selected  areas 
with  He  at  10*^  cm”^  and  30  keV.  (a)  shows  the  best  fit 
collection  probability  for  the  written  area  of  P12B  and 
target  P12C  which  functions  as  a control.  (b)  shows  the 
gain  versus  beam  energy  data  and  the  best  fit  gain  curves 
calculated  from  the  collection  probability  in  Figure  (a). 
The  solid  circles  are  the  data  from  P12C  and  the  open 
circles  from  the  unwritten  area  of  P12B.  The  triangles 
are  the  data  from  the  written  areas  of  P12B.  In  (a) 
the  calculated  position  of  the  He  implant  peak  R and  the 
standard  deviation  in  the  projection  range  &R  P are 
indicated.  The  damaged  region  should  coincide  with  the 
implanted  region.  Note  that  the  region  of  reduced  collec- 
tion probability  agrees  well  with  the  predicted  region 
of  damage  generated  by  the  He. 


(See  pages  174  U 175) 
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Figure  97b.  Gain  curves  calculated  from  best  fit  collection 
probability  in  Figure  97(a)  and  gain  versus 
energy  data. 
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in  the  minority  carrier  lifetime  in  the  written  area  due  to 
damage  by  the  writing  ion  beam.  This  prevents  carrier  collection 
by  the  planar  diode.  In  effect,  this  is  a written  dead  layer. 

An  alternative  mechanism  may  depend  on  altering  the  doping  level 
of  the  written  area  by  the  introduction  of  deep  traps.  These 
traps  deplete  the  carrier  density  in  the  n+  layer  causing  a 
negative  field  which  effectively  traps  the  generated  carriers 
until  they  recombine.  That  is,  it  reduces  or  eliminates  the 
electric  field  generated  by  the  doping  gradient.  If  acceptor 
levels  are  introduced,  it  may  be  possible  to  cause  type 
conversion  although  this  is  not  known  to  happen  for  He  implants. 

The  success  of  damage  writing  will  require  rethinking  of 
some  of  the  specifications  placed  on  doping  levels  and  fabrication. 
In  particular  depletion  field  spreading  of  the  bits  may  not  be  as 
severe  a problem.  Higher  writing  rates  may  be  possible,  since 
potentially  less  fluence  may  be  needed  to  damage  the  semiconductor 
than  to  dope  it  to  an  electrical  activity  sufficient  to  produce 


diode  bits. 


Future  Plans:  The  writing  effects  of  other  ions,  particularly 
the  inert  gas  ions  will  be  explored,  as  well  as  the  limitations 
on  writing  fluence.  Degradation  mechanisms  will  also  be  explored 
including  readout  fatigue  of  the  stored  patterns  and  room  tem- 
perature annealing  of  the  damage  writing.  The  doping  levels  in 
the  planar  structure  will  be  changed  as  well. 
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Section  VI 

ALLOY  JUNCTION  WRITING 

Selection  of  Metal  Semiconductor  Systems  for  the  Alloy  Junction  Target 

An  alloy  junction  target  is  shown  in  cross  section  in  Figure  98. 
To  write,  an  area  is  heated  locally  by  an  electron  beam  (or  laser 
beam  in  evaluation  tests)  to  cause  local  alloying  of  the  metal- 
semiconductor.  In  the  idealized  concept,  the  heated  region  forms  a 
pn  junction.  This  can  then  be  read  by  the  electron  beam  because  of 
its  effect  on  the  generated  carriers.  This  section  discusses  the 
criteria  for  the  choice  of  the  best  semiconductor  metal  combinations 
and  identifies  the  most  promising  binary  systems. 

Alloy  Junction  Formation 

Figure  99  shows  the  motion  of  carriers  during  readout  of  a 

written  bit.  When  the  beam  impinges  on  the  written  area,  there 

should  be  only  a small  readout  signal  (i  ) generated  in  the  main 

s 

pn  junction  circuit.  For  this  to  occur  the  carriers  generated  by 
the  electron  beam  must  either  recombine  in  the  local  region  where 
they  are  generated  or  not  be  generated  at  all.  Local  recombination 
can  occur  if  the  written  area  is  actually  a pn  junction  as  shown  in 
Figure  99a.  In  this  embodiment  electron-hole  pairs  are  trapped  by 
the  diode  depletion  field  until  they  recombine  by  some  process, 
for  example,  recombination  at  an  ohmic  contact  across  the  top  of 
the  structure.  The  ohmic  contact  or  short  could  be  a metal  contact 
or  just  a heavily  doped  semiconductor  layer. 
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Figure  99.  Writing  Methods 
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Another  way  to  cause  local  recombination  is  to  cause  the  written 
region  to  have  a very  short  minority  carrier  lifetime  (Figure  99b). 
This  could  be  accomplished  by  damage  to  the  lattice  (introduced, 
for  example,  by  the  rapid  regrowth  of  the  semiconductor  from  the 
semiconductor-metal  melt  or  by  a heavy  concentration  of  dissolved 
metal. 

An  alternative  writing  and  readout  method  would  depend  on  preven- 
ting electron-hole  pair  generation  in  the  written  region  because  the 
area  is  actually  a metal  or  * metal  with  semiconductor  crystallites 
imbedded  (Figure  99c) . 

All  of  these  writing  methods  depend  upon  the  formation  of  a 
liquid  metal- semiconductor  melt.  The  solubility  of  the  semiconductor 
in  the  melt  must  be  large  for  good  writing.  Figure  100  shows  the 
formation  of  the  melt  above  some  critical  temperature  T^.  The 
solubility  by  weight  (w)  of  the  semiconductor  in  the  melt  is: 

fsfs  \ 

pMtM  + Dsts  J (41) 

Where  t is  the  thickness  of  the  dissolved  semiconductor,  tx.  the 
initial  metal  layer  thickness  and  pg  and  p^  the  mass  densities  of 
the  semiconductor  and  metal,  respectively.  Equation  41  can  be  solved 
for  the  ratio  thickness  of  dissolved  semiconductor 

ts/tM  = (pM/Ps)(l^) 

(42) 
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During  cooling  of  the  written  region  the  melt  will  solidify  leaving 
an  affected  region  of  depth  tg  in  the  semiconductor. 

During  cooling  any  of  the  situations  shown  in  Figure  99  could 
occur  depending  upon  the  material  properties  and  cooling  rates.  How- 
ever, for  all  three  cases  the  thickness  of  the  affected  region  tg 
should  be  large  to  produce  a large  reduction  in  readout  diode  gain. 

(It  is,  of  course,  possible  that  none  of  the  three  situations  occurs 
if  the  semiconductor  regrows  with  very  little  damage  or  doping.  This 
is  unlikely,  expecially  because  of  the  rapid  cooling  rates.)  It  is 
therefore  important  to  choose  metal-semiconductor  pairs  where  tg  is 
large.  Tables  5 and  6 summarize  the  properties  of  a large  number 
of  metal-Ge  and  metal-Si  pairs.  All  of  these  pairs  (except  As  which 
is  not  really  a metal)  have  simple  eutectic  phase  diagrams.  As  an 
example,  the  phase  diagram  for  Al-Si  is  reproduced  in  Figure  101. 

All  of  the  data  in  Tables  5 and  6 is  taken  from  Hansen  and  its 
supplements  (Refs.  15,  16,  17).  This  collection  is  believed  to 
represent  all  of  the  eutectic  pairs  for  Ge  and  Si  with  eutectic 
temperatures  below  700°C.  III-V  and  II-IV  semiconductors  might  also 
be  considered  but  these  are  less  well  characterized  and  are  of  course 
trinary  systems  which  are  more  complicated.  They  need  not  be  reviewed 
unless  none  of  the  Ge  and  Si  systems  are  successful. 

Note  that  only  a few  pairs  have  tg/t^  greater  than  even  1 °L. 

For  most  pairs  it  is  much  less  and  was  not  even  calculated.  Based 
on  the  criteria  of  large  tg/t^  only  8 pairs  are  candidates:  Al-Ge, 
Al-Si.  Au-Ge,  Au-Si,  Ag-Ge,  Ag-Si,  Sb-Ge,  and  Zn-Ge. 
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Thermal  bias  of  the  target  will  probably  be  necessary  because 
of  the  limited  power  of  the  electron  beam  and  the  high  thermal 
conductivity  of  the  semiconductors.  Since  a minimum  of  several 
hours  is  required  to  write  a 10H  bit  memory  chip,  the  target  will 
have  to  tolerate  extended  periods  (hours  or  days)  at  a temperature 
tens  of  degrees  below  the  eutectic  temperature  and  probably  several 
cycles  between  the  thermal  bias  temperature  and  room  temperature. 

Hence  any  pair  where  the  diffusion  or  solid-phase  reactions  can 
occur  is  less  attractive. 

Diffusion 

The  diffusion  constants  30°C  below  the  eutectic  temperature  are 

tabulated  in  Tables  V and  VI.  In  many  cases  these  are  very  extended 

extrapolations  of  much  higher  temperature  data  and  so  may  be  grossly 

in  error.  In  10^  seconds  (2.8  hours)  a diffusion  constant  of  10  ^ 

2 

cm  /sec  would  result  in  a 100  X diffusion  depth.  This  could  begin 

'"16  2 

to  cause  some  spontaneous  writing  so  D < 10  cm  /sec  is  unattractive. 
Based  on  this  criteria  Ag-Ge,  Ag-Si,  and  Sb-Ge  can  be  clearly  elim- 
inated. The  diffusion  constants  are  fairly  well  known  for  these 
pairs  because  the  eutectic  temperatures  are  so  high.  The  high  temp- 
eratures these  systems  have  also  makes  them  less  attractive.  Au-Si 
would  also  appear  to  be  ruled  out.  However,  because  the  diffusion 
constant  is  not  as  well  known  at  this  temperature,  it  should  not  be 
completely  eliminated. 

So lid -So lid  Reactions 

Solid  phase  reactions  are  another  potential  problem.  It  is  known 


1 
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that  Al-Ge  and  Al-Si  (Refs.  18  - 23)  films  will  dissolve  and  migrate 

at  temperatures  much  below  the  reaction  temperature  until  the  solid 

solubility  of  the  semiconductor  in  metal  is  reached.  For  Al-Ge 

a 3000  8 A1  film  reacts  completely  in  20  min.  at  less  than  300°C, 

even  though  the  eutectic  temperature  is  424°C.  For  Al-Si  the  reaction 

(21) 

is  affected  by  impurities  and  is  aided  by  the  presence  of  poly- 

(22) 

crystalline  or  amorphous  silicon  but  the  reaction  does  occur  on 

(22)  (19) 

single  crystals.  For  Sb-Ge  the  reaction  apparently  does  not 

occur  after  1 hour  at  280°C.  The  addition  of  an  amorphous  Ge  layer 

(19) 

inside  the  Sb  layer  causes  the  reaction  to  occur 

Solid-solid  reactions  are  generally  undesirable  for  this  applic- 
ation. The  solid  solubility  of  semiconductor  in  the  metal  is  gener- 
ally low  as  shown  in  Tables  V and  Vt.  Upon  cooling  a doped  semicon- 
ductor layer  will  regrow  or  the  semiconductor  may  precipitate  out 
in  the  metal.  If  the  solid  solubility  is  large,  the  regrown  doped 
region  can  be  large  and  produce  an  entirely  written  target.  If  the 
dissolution  and  regrowth  are  nonuniform,  each  cycle  between  room  temp- 
erature and  thermal  bias  temperature  will  cause  greater  nonuniformity 
and  eventually  permanent  damage  to  the  storage  surface.  Hence,  it 
will  be  necessary  to  find  methods  to  prevent  solid-solid  reactions 
or  at  least  ensure  that  they  occur  uniformly.  This  points  to  the 
desirability  of  a diffusion  barrier  which  breaks  down  or  melts  near 
the  eutectic  temperature  of  the  metal-semiconductor  pair.  It  will  be 
necessary  to  study  these  solid-solid  reactions  since  they  are  not  well 
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characterized  except  for  the  A1  systems.  Since  the  solubility  is 
low  (<3%)  for  all  the  candidate  systems,  there  is  no  reason  to  rule 
out  any  pair  until  more  is  known  about  the  reactions.  Solubilities 
<57,  should  not  prevent  feasibility  demonstrations. 

Electron  Beam  Penetration 

During  readout,  the  electron  beam  must  penetrate  through  the 
metallization.  The  metallization  could  be  stripped  before  readout, 
but  this  is  undesirable  because  verification  of  the  written  data 
would  be  impossible  until  the  entire  plate  was  written  and  stripped 
of  metal.  The  beam  energy  can  be  increased  to  penetrate  through  the 
metal  layer,  but  this  will  increase  the  lateral  scattering  of  the 
beam  in  the  semiconductor.  This  effect  is  sketched  schematically 
in  Figure  102.  To  a first  approximation,  the  electron  range  and 
lateral  scattering  is  proportional  to  the  mass  density  for  a given 
beam  energy.  If  the  beam  has  sufficient  energy  to  penetrate  into 
the  semiconductor  then  the  penetration  and  lateral  scattering  in 
the  semiconductor  will  be  approximately  t = R - t^  (45) 

P f>s 

as  shown  in  Figure  102,  where  R is  the  electron  range  in  the  semi- 

s 

conductor.  This  is  only  a reasonable  approximation  if  tp/Rg>  1 (46) 

so  that  lateral  scattering  in  the  metal  is  negligible.  To  achieve 
a bit  spacing  B,  the  lateral  scattering  in  the  semiconductor  must 
be  less  than  B/2  so  that: 

Ra  ' Wk  - S s 8/2  (47) 
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To  satisfy  equation  46  and  assure  a significant  amount,  beam  energy 
is  dissipated  in  the  semiconductor.  It  is  assumed  that: 

R = 2t  (48) 

s p 

Solving  equations  42,  45,  47,  and  48  for  t yields: 

s 

t s B/2  (_w ) (49) 

1-w 

Note  that  equation  49  is  independent  of  the  mass  density  of  the 
metal.  The  reason  is  that  for  the  same  solubility  a thinner  layer 
of  a dense  metal  can  dissolve  more  semiconductor  and  compensate  for 
the  thinner  metal  layer  required  to  permit  adequate  beam  penetration 
into  the  semiconductor. 

To  assure  a large  modulation,  t should  be  as  a significant 

s 

fraction  of  the  electron  penetration  into  the  semiconductor  and  must 
also  be  larger  than  the  dead  layer  thickness.  For  both  of  these 
reasons  tg  > 0.02  microns  is  a reasonable  requirement.  Hence  for  a 
bit  spacing  of  0.1  microns,  equation  49  can  be  solved  for  the  minimum 
required  solubility.  This  gives  w a 307«;  this  calculation  is  only 
a rough  approximation  but  it  points  to  solid  solubilities  much  less 
than  10%  as  very  undesirable.  Based  on  this  criteria  the  Al-Ge  is  by 
far  the  most  favorable  with  Au-Ge,  Al-Si  the  next  most  favorable 
choices,  of  the  remaining  uneliminated  pairs. 

Choice  of  the  Most  Favorable  Pairs 

Based  on  the  above  discussion  the  metal-semiconductor  pairs 
listed  in  Tables  V and  VI  can  be  reduced  to  Al-Ge,  Au-Ge  and  Al-Si 
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with  Al-Ge  by  far  the  most  attractive.  Au-Si  should  be  considered 
if  the  diffusion  constant  of  Au  in  Si  is  smaller  than  predicted.  The 
In-Ge  system  (Ref.  21,  22)  does  not  look  attractive.  Sb-Ge  is  also 
possible  although  the  diffusion  constant  of  Sb  in  Ge  is  too  large 
based  on  published  data. 

Alloy  Junction  Laser  Diode  Experiments 

An  established  method  for  producing  p-n  diodes  involves  alloying 
metal  dopant  material  on  an  oppositely  doped  semiconductor.  For 
example  heating  of  an  indium  (p-type  dopant)  on  n-type  germanium 
forms  a local  alloy  melt  as  the  temperature  reaches  approximately 
160°C,  in  which  the  proportion  of  indium  to  germanium  exceeds  the 
solid  solubility  of  indium  in  germanium  (Figure  103) . Upon  cooling, 
indium  doped  germanium  regrows  at  the  alloy -germanium  interface  by 
a process  similar  to  liquid  epitaxy,  resulting  in  a junction  between 
the  p-  and  n-type  germanium.  Such  junctions  are  often  produced  under 
conditions  that  inhibit  oxidation  of  the  semiconductor  and  the  alloy 
during  heating,  such  as  a reducing  atmosphere,  so  that  oxide  at  the 
interface  does  not  hinder  alloying. 

This  technique  has  been  proposed  as  a writing  method  for  the 
surface  diode  archival  memory  concept  to  form  the  tiny  p-n  diodes 
in  the  surface  of  a planar  n-p  junction  structure  as  shown  in  Figure 
104. 

The  writing  method  for  producing  alloy  surface  diodes  is  illus- 
trated in  Figure  105.  In  this  process  the  surface  of  the  n-p  planar 
structure  is  coated  with  a thin  layer  of  the  desired  metal,  and  a 
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Figure  103.  Formation  of  Alloy  Junction  with  Indium  on  N-Type 
Germanium 
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Figure  104.  Cross  Section  of  Alloy  Junction  Diode 
Archival  Memory  Project 
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Figure  105.  Focused  Electron  Beam  Method  for  Formation 
of  Alloy  Junctions 


focused  electron  beam  is  used  to  locally  heat  a bit  region  to  the 
temperature  required  for  alloy  action.  The  entire  structure  is 


thermally  biased  to  a temperature  near  the  eutectic  temperature, 
so  that  less  local  heating  by  the  electron  beam  would  be  required 
for  alloy  formation,  thus  reducing  the  writing  beam  requirements. 

The  metal  coated  n-p  diode  blanks  can  be  prepared  using  well 
established  processing  techniques  (e.g.  ion  implantation  or  diffusion 
of  an  n or  p layer)  prior  to  insertion  in  the  writing  station.  Clean- 
ing the  semiconductor  surface  in  vacuum  just  before  deposition  of 


the  metallic  coating  reduces  the  possibility  of  interface  oxidation 


inhibiting  the  alloy  phenomena.  Special  intermediate  layers  to 


reduce  oxide  formation  and  thereby  increase  the  rate  of  eutectic 
formation  might  also  be  used.  These  would  have  a function  similar 
to  that  of  flux  in  a soldering  process. 

A number  of  III-V  dopant -semiconductor  systems  exist  that  might 
be  applicable  to  this  alloy  writing  technique,  for  which  the  roles 
of  the  n and  p regions  could  be  interchanged.  The  most  promising 
of  these,  due  to  their  relatively  low  eutectic  temperatures  are: 

System  Eutectic  Temperature  (°C) 


Indium-Geranium  156 
Antimony-Germanium  590 
Thallim-germanium  302 
Aluminum- germanium  424 
Aluminum- silicon  577 
Antimony-silicon  630 
Indium-silicon  156 
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Germanium  systems  have  the  advantages  of  a thermal  diffusivity 
one-third  that  of  silicon  and  surfaces  less  prone  to  oxidation. 

To  estimate  the  amount  of  beam  current  necessary  to  heat  a local 
volume  of  the  alloy  junction  target,  the  simplified  model  shown 
in  Figure  106  is  used.  The  beam  energy  is  assumed  to  be  completely 
dissipated  within  a hemisphere  of  radius  a.  Radius  a is  determined 
by  the  electron  scattering  in  the  material  and  by  spot  size  b.  A 
reasonable  approximation  to  a is : 

2a  = b + Rg  (50) 

where  R^  is  the  Grun  range  which  is  a measure  of  electron  penetration 
into  solids.  Figure  107  shows  the  dependence  of  R^  on  beam  voltage 
for  several  materials  of  interest. 

The  solution  of  the  heaL.  diffusion  equation  for  a uniformly 
heated  hemisphere  at  the  surface  of  a semiinfinite  medium  has  been 
solved  by  Goldenberg.  Figure  108  shows  the  temperature  distribution 
in  a target  at  various  times  as  a function  of  radial  distance  r.  Temp- 
erature and  time  are  normalized  as  follows: 


AT0  = 3P/47TKa 

(51) 

r - a2/k 

(52) 

where  k = thermal  diffusivity  of  the  medium,  K « thermal  conductivity, 
and  P = total  thermal  power  into  the  hemisphere  of  radius  a.  The 
thermal  parameters  of  silicon,  germanium,  aluminum,  and  indium  are 
given  in  the  Table  VII.  Note  that  the  parameters  for  indium  and 
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Figure  106. 


Thermal  Writing  Model  for  Alloy  Junction  Writing 


Table  VII 

THERMAL  PARAMETERS 


Element  Specific  Heat  Thermal  Condition  Thermal9Dif fusivity 

(Cv  j/g°C) (K(W/cm°C) (K  cm  /sec) 

0.9 

0.36 

0.29 

0.38 

and  germanium  are  nearly  identical,  so  a target  sandwich  constructed 
of  these  materials,  can  be  treated  as  thermally  uniform.  For  other 
combinations,  such  as  aluminum-germanium,  the  situation  is  more 
complicated  and  the  solution  of  Figure  108  does  not  exactly  apply. 

From  Figure  108,  it  is  evident  that  for  a given  beam  power  the 
highest  temperature  can  be  achieved  for  large  t/r,  without  very  much 
loss  in  spatial  resolution,  where  t is  the  time  needed  to  write  one 
bit.  For  germanium  or  indium  the  thermal  time  constant,  t,  is 
3 x 10  ^ seconds  for  a = 0.1  micron.  Therefore,  using  a beam  dwell 
time  that  is  long  compared  to  t will  only  limit  the  writing  rate 
to  <300  Mbits /sec.  Assuming  a thermal  bias  such  that  the  desired 
temperature  rise,  AT,  is  >15°C,  Figure  109a  indicates  the  beam  power 
required  for  an  asymptotic  temperature  rise,  aTq  = 30°C,  at  the 
center  of  the  spot,  as  a function  of  hemisphere  radius  a.  Bit  spac- 
ing B as  identified  with  hemisphere  diameter  2a.  The  approximate 
size  of  the  heated  area  will  exceed  the  alloy  temperature  if  the 
thermal  bias  is  chosen  ATq/2  below  the  alloy  temperature  (Figure  108). 

-200-  ] 

„W,;a.W  ...  , II  m^irn  j 


Silicon 

0.7 
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Germanium 

0.31 
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Figure  109a.  Beam  Power  Required  for  Asymptotic  Temperature 
Rise  of  30°C  Versus  Bit  Spacing  for  Germanium 
Alloy  Diode  Target 
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The  required  beam  current  and  beam  voltage  can  then  be  determined 
from  equation  50  as  a function  of  spot  size  b.  This  is  shown  in 
Figure  109b  for  two  bit  spacings,  B.  Based  on  these  requirements 
an  electron  column  for  alloy  junction  diode  formation  has  been 
purchased.  The  details  of  the  column  and  its  selection  are  given 
in  Section  VII  of  this  report. 

To  develop  a feel  for  power  requirements  and  the  effect  of  ther- 
mal bias,  a laser  system  has  been  used  for  the  radiation  source 
required  to  fffrm  the  surface  diodes.  Although  there  are  significant 
differences  in  the  nature  of  the  laser  and  electron  beam  interaction 
with  the  target,  the  laser  experiments  can  provide  useful  data  on 
the  feasibility  of  the  alloy  junction  concept.  Questions  concerning 
the  time  of  formation  and  power  required  to  form  the  diodes  can  be 
answered  in  part  if  not  in  whole  regardless  of  the  radiation  differ- 
ences. The  details  of  the  laser  experiments  performed  thus  far  are 
discussed  in  the  following  paragraphs. 

Laser  System 

The  laser  system  used  for  the  alloy  junction  experiments  is  an 
Electro  Scientific  Industries  Model  25  laser  trimmer.  The  system 
uses  a continuously  pumped  Nd  doped  YAG  laser  with  an  acousto-optic 
Q-switch.  It  is  capable  of  delivering  at  le*st  3.0  watts  C.W.  in 
the  TEMqq  mode  and  up  to  30  watts  C.W.  multimode  in  Q-switched  oper- 
ation. The  Q switch  modulator  itself  is  powered  by  a 50  Mhz.  RF 
driver  unit.  The  Q switch  frequency  or  pulse  repetition  rate  is 
continuously  variable  from  1 to  20  KHz  with  a typical  pulse  time  of 
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150  ns.  The  laser  can  be  apertured  between  the  reflecting  surfaces 
to  limit  the  higher  order  modes.  These  apertures  also  influence 
the  final  spot  size  since  the  optics  are  diffraction  limited.  Various 
optical  configurations  are  possible: 

• 2 lens  system-condensor  and 
objective  « .8  mil  final  spot; 

*10  mil  depth  of  focus. 

• 3 lens  system-condensor,  projector 
and  objective  ss2-3  mil  final  spot; 

-30  mil  depth  of  field 

• 3 lens  system-condensor  and 
projector  lens  switched  - final 
spot  less  than  1 mil. 

For  the  alloy  junction  diode  experiments  the  optics  des- 

cribed last  above  were  used  exclusively.  Final  output  power  is  a 
function  of  lamp  current,  limiting  aperture,  optical  chain  and  Q 
switch  frequency. 

The  system  is  equipped  with  a manually  operated  two-position 
handler.  Each  position  has  vacuum  hold  down  for  securing  chips  up 
to  2 1/4  by  2 1/4  inches  in  size.  The  two  position  handler  allows 
the  operator  to  load  or  unload  substrates  while  a second  substrate 
is  being  exposed  to  the  laser  in  the  other  next.  The  handler  has 
positional  repeatability  of  *0.001  inch. 

The  laser  beam  is  positioned  on  the  sample  by  mirrors  and  lenses 
mounted  on  an  X-Y  table  which  is  actuated  by  linear  motors.  In 
moving  over  the  substrate  area  the  laser  beam  remains  perpendicular 
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to  the  work  surface  resulting  in  no  change  in  focus  or  beam  energy. 
Position  of  the  laser  beam  can  be  either  manual  or  under  the  process 
control  computer  (Digital  Equipment  Corporation  PDP8  series). 

Targets 

Alloy  junction  target  blanks  were  fabricated  from  BEAMOS  memory 
planes.  The  basic  BEAMOS  target  structure  that  was  used  is  shown  in 
crossection  in  Figure  110.  It  consists  of  a p substrate  with  an 
expitaxially  grown  n layer  covered  with  an  oxide  and  aluminum  layer. 
Typical  parameters  are: 

A1  layer  thickness 
Oxide  thickness 
n epitaxi  ~2.0  urn 
n layer  resistivity  . 9 A - cm 
p layer  resistivity  100  A - cm 

The  BEAMOS  target  structures  were  modified  slightly.  First,  the 
aluminum  layers  were  removed  with  aluminum  etch  and  the  oxide  strip- 
ped with  hydroflouric  acid.  The  silicon  target  surfaces  were  then 
plasma  etched  and  sputter  coated  with  a thin  aluminum  layer.  These 
last  two  steps  were  done  sequentially  under  vacuum  to  prevent  oxide 
growth  on  the  silicon  surface.  Targets  with  sputtered  aluminum 
layers  of  500  X,  1000  8,  and  2000  $ were  fabricated.  All  tests  to 
date  have  been  performed  on  the  targets  with  1000  X layers . 

Laser  Writing  Experiments 

As  a first  step  the  laser  experiments  were  aimed  at  determining 
the  operating  conditions  for  diode  formation.  The  process  control 
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Figure  110.  Cross  Section  of  BEAMOS  Target 
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computer  was  programmed  to  produce  the  test  patterns  shown  in  Figure 
111.  The  laser  is  held  stationary  at  each  spot  and  gated  on  for  the 
time  duration  indicated  above  each  dot  on  the  figure.  Lamp  drive 
current,  Q switch  frequency  and  laser  apertures  were  manually  changed 
variables.  The  number  of  pulses  delivered  to  each  spot  is  deter- 
mined by  the  gate  on  time  and  the  Q switch  frequency.  Power  delivered 
to  the  target  depends  on  the  lamp  drive  and  the  optical  aperturing. 
The  first  tests  were  conducted  with  a 20A  lamp  drive  and  a 10% 
transmitting  beam  splitter  in  the  optical  path.  These  conditions 
were  set  so  the  number  of  variables  would  be  limited  and  their 
values  were  based  on  laser  trimming  experience.  The  Q switch 
frequency  and  the  laser  aperture  size  were  varied  over  the  time 
durations  shown  in  Figure  111.  The  best  operating  conditions  were: 

L^  = lamp  current  = 20A 

107o  Beam  Splitter 

f q = Q switch  frequency  = 2kHz. 

d^  = Laser  aperture  diameter  = 60  mils 

Memory  readout  was  obtained  on  the  Cambridge  Stereoscan  Scanning 
Electron  Microscope.  The  target  planar  diode  signal  was  used  for 
readout  as  shown  in  Figure  112.  Memory  readout  from  a laser  irrad- 
iated spot  before  removal  of  the  1000  8 A1  layer  is  shown  in  Figure 
113.  Figure  114  shows  readout  of  the  same  area  after  the  aluminum 
layer  was  etched  off.  These  two  figures  verify  the  initial  belief 
that  readout  could  be  achieved  with  or  without  removal  of  the 
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Figure  113. 


Planar  Diode  Signal  Showing  Laser  Alloyed 
Region  on  Si  Substrate  Coated  with  1000  ° Al. 
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Figure  114.  Planar  Diode  Signal  Showing  Laser  Alloyed 

Region  on  Si  Substrate  After  Removal  of  1000  X 
A1  Layer. 
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aluminum  layer.  A secondary  electron  surface  scan  showing  the  top- 
ography of  this  irradiated  area  is  shown  in  Figure  115.  The  indented 
crater  with  raised  center  as  indicated  by  this  figure  seemed  to  be 
characteristic  of  areas  with  good  readout.  The  existence  of  aluminum 
within  the  laser  irradiated  area  was  verified  by  scanning  electron 
microscope  nondispersive  x-ray  analysis.  The  results  of  this  analysis 
within  the  irradiated  area  are  shown  in  Figure  116.  Here  the  marker 
at  1.74  keV  denotes  the  Silicon  peak.  At  the  lower  value  of  1.49 
keV  the  much  smaller  aluminum  peak  is  found.  The  same  type  of  analy- 
sis from  the  area  just  outside  the  irradiated  spot  is  shown  in  Figure 
117.  As  seen  here  only  the  marker  at  1.49  keV  remains  indicating 


the  absence  of  aluminum  outside  of  the  irradiated  area.  Although 
not  conclusive  proof  the  results  from  the  diode  read  and  x-ray  analy- 
sis give  strong  indications  that  an  alloy  junction  surface  diode  was 
formed.  The  actual  crater  area  from  which  memory  readout  was  obtained 
was  on  the  order  of  15  um  across.  For  this  particular  spot  the  laser 


was  gated  on  for  1.56  seconds  corresponding  to  3120  pulses  at  the 
Q switch  frequency  of  2kHz.  The  memory  readout  and  the  video  wave 
form  for  two  spots  written  with  fewer  pulses  is  shown  in  Figure  118. 
Here  the  laser  gate  on  time  was  0.083  and  0.162  sec  for  the  spot  on 
the  left  and  right  respectively.  At  the  2kHz.  Q switch  rate  these 


correspond  to  166  and  324  pulses.  The  video  waveform  corresponds  to 
1.0  volt/cm  and  demonstrates  good  discrimination  between  written  and 


unwritten  areas. 
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Figure  116.  SEM  Nondispersive  X-Ray  Analysis  of  laser  Alloyed 
Region  on  Si  Substrate  after  removal  of  1000  X 


A1  layer. 
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SEM  NONDISPERSIVE  X-RAY  ANALYSIS  OF  Si  SUBSTRATE 
NEAR  LASER  ALLOYED  REGION  AFTER  REMOVAL  OF  1000  A 
Al  LAYER 


Figure  117.  SEM  Nondispe rs ive  X-Ray  Analysis  of  Si  Substrate 
near  Laser  Alloyed  Region  after  removal  of  1000  8 
Al  layer. 


Figure  118.  Video  Waveform  avid  ' ier  , icde  Signal  from  two 
adjacent  Laser  M1<  v ’ i v-gions. 


One  of  the  purposes  of  the  laser  experiments  was  to  demonstrate 
that  alloy  diodes  could  be  formed  on  a short  time  scale  (i.e.  10  MHz 
writing  rate).  The  preceeding  experiments  did  not  do  this  because 
of  the  large  number  of  pulses  involved.  However,  the  model  used 
for  the  alloy  junction  formation  in  the  introduction  to  this  section 
indicates  that  the  diodes  should  form  in  one  pulse  if  they  form  at 
all.  Hope  is  added  to  this  concept  due  to  the  fact  that  diode  read- 
out was  just  as  good  from  an  area  irradiated  with  166  pulses  as  one 
irradiated  with  3120  pulses.  Therefore,  the  next  experiments  were 
dedicated  to  reducing  the  number  of  laser  pulses  per  spot  with  the 
final  goal  of  writing  with  one  pulse. 

As  a starting  point  for  this  goal  the  process  control  computer 
was  reprogrammed  to  write  an  NxN  array.  The  size  of  the  array,  in 
terms  of  spot  to  spot  center  distance  and  number  of  spot  N as  well 
as  the  laser  "gate  on"  time  per  pulse  were  programmable  inputs.  A 
new  laser  operating  mode  was  sought  using  the  40  mil  aperture  so 
that  the  spot  size  would  be  smaller.  With  this  aperture  it  was  found 
that  diodes  were  formed  with  a lamp  drive  current  of  18  A and  a Q 
switch  frequency  of  2kHz.  The  10%  transmissive  beam  splitter  was 

also  removed  from  the  optical  system. 

4 

Planar  diode  readout  of  a 10  spot  array  formed  with  these  oper- 
ating conditions  and  10  pulses  per  spot  is  shown  in  Figure  119. 
Considerable  variability  from  point  to  point  and  even  the  absence 


of  readout  at  some  points  can  be  seen.  A magnified  scan  of  part  of 
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Figure  119.  Diode  Current  Signal  10  Spot  Array  - SEM  Readout 
50X  - Laser  Irradiated  - 10  pulses  per  spot  - 1000 
8 A1  on  Si. 
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this  array  showing  a missing  location  and  the  variability  from 
point  to  point  more  clearly  is  given  in  Figure  120.  A secondary 
emission  surface  scan  of  this  same  area  is  shown  in  Figure  121.  With 
this  figure  the  variability  in  the  diode  signal  can  be  linked  directly 
to  surface  nature  variations  of  the  irradiated  spots.  Note  that  at 

the  missing  signal  location  there  was  apparently  no  laser  pulse  as 

S . * * '•  ' r * * / f> 

the  surface  scan  shows  no  radiation  damage.  All  of  these  variations 

are  attributed  to  variations  in  the  laser  output  due  to  a poorly 
regulated  power  supply  and  an  aging  lamp.  More  will  be  said  about 
this  when  the  laser  power  measurements  are  discussed. 

Figure  122  shows  the  planar  diode  signal  of  a portion  of  this 
array  at  still  higher  SEM  magnification.  The  spots  are  on  20  urn 
centers  indicating  an  irradicated  spot  size  which  is  typically  10 
urn  in  diameter.  The  secondary  emission  surface  scan  of  the  same 
3x3  spots  is  shown  in  Figure  123.  At  this  magnification  the  granu- 
larity of  the  sputtered  aluminum  layer  is  clearly  visible  in  this 
figure.  This  same  granularity  shows  up  but  at  a lesser  degree,  in 
the  planar  surface  diode  output  of  Figure  122.  The  reason  for  this 
is  that  the  beam  is  stopped  at  some  points  and  therefore  no  electron 
hole  pairs  are  generated  and  hence  no  diode  signal.  As  far  as  planar 
diode  signal  this  result  is  the  same  as  at  an  irradiated  spot  where 


the  excess  minority  carriers  experience  enhanced  recombination. 

Finally  a video  waveform,  or  the  actual  planar  diode  current  for  the 
center  three  spots  of  this  portion  of  the  array  is  shown  in  Figure  124. 
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Figure  120.  Diode  Current  Signal  from  10  Spot  Array  - SEM 
Readout  500X  - Laser  Irradiated  - 10  pulses  per 
spot  - 100  R A1  on  S:i 
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Figure  121.  Secondary  Emission  Surface  Scan  from  10  Spot 
Array  - SEM  Readout  500X  - Laser  Irradiated  - 
10  Pulses  per  Spot  - 1000  8 A1  on  Si. 
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Figure  122.  Diode  Current ^Signal  oi  i x 3 Laser  irradiated 
Spots  from  ] 0 Spot  Array  - SEM  Readout  2000  X 
- Laser  Irradiated  - 10  Pulses  per  Spot  - 1000 
A1  on  Si. 
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Figure  124.  Video  Waveform  from  3 Adjacent  Alloy  Junction  Diodes  - 
Laser  Irradiated  - 10  Pulses  per  Spot  - 1000  A A1  on 
Si. 
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The  noise  between  diode  signals  if  just  the  result  of  the  granul- 
arity of  the  sputtered  aluminum  in  this  presentation. 

The  10  pulses  per  spot  shown  in  the  preceeding  figures  was  the 
minimum  number  of  pulses  that  could  be  programmed.  The  next  step 
was  to  try  single  pulses  and  therefore  a program  change  had  to  be 
made.  The  process  control  computer  was  reprogrammed  so  that  the 
table  moved  with  a constant  velocity.  Based  on  the  10  um  spot  size, 
as  long  as  the  table  velocity  was  greater  than  20  mm/sec,  single 
pulses  at  a 2kHz  rep  rate  would  be  resolved.  Figure  125  shows 
readout  through  the  planar  diode  of  three  successive  single  pulse 
irradiated  areas.  The  laser  operating  conditions  were  the  same  as 
the  preceeding  experiment;  i.e.  40  mil  aperture,  18  A lamp  current, 
2kHz  rep  rate  and  the  table  velocity  was  50  mm/sec.  The  irradiated 
spots  are  on  25  um  centers  and  have  a typical  diameter  of  7.5  um. 

The  difference  in  spot  size  between  these  and  the  10  pulse  per  spot 
results  is  attributed  to  the  changing  surface  conditions  and  increased 
damage  effects  due  to  multi  pulses.  The  planar  diode  current  signal 
or  video  waveform  from  a scan  through  the  three  pulses  is  shown  in 
Figure  126. 

Laser  Power  Measurements 

The  beam  power  delivered  to  the  target  in  the  model  25  system 
was  measured  for  various  operating  conditions.  This  data  is  needed 
in  order  to  obtain  a calibration  on  the  temperature  rise  achieved 
in  the  preceeding  experiments  as  well  as  for  future  experiments.  Of 
course  the  full  power  delivered  to  the  target  is  not  utilized  and 
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Figure  125.  Readout  Through  Piniar  hi  ode  ot  Three  Single  Pulse 
Irradiated  Areas . 
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Figure  126.  Waveform  Scan  Through  Three  Irradiated  Areas. 
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must  be  modified  by  the  reflectance  of  the  surface.  The  reflectance 
data  is  still  under  investigation.  At  present  we  have  determined 
the  power  delivered  to  the  target  as  a function  of  the  laser  system 
parameters . 

Average  power  was  measured  at  the  target  plane  with  a Coherent 
Radiation  model  205  power  meter.  The  model  205  consists  of  a direct 
absorption,  convection-cooled  head  and  a control-indicator  unit.  The 
instrument  is  specified  for  CW  laser  measurements  covering  the  broad 
wavelength  range  of  .3  to  30  microns.  A relative  measure  of  peak 
pulse  power  was  obtained  from  a photo  diode  that  looks  at  the  laser 
beam  in  the  model  25  laser  system.  The  signal  from  this  photo  diode 
was  displayed  on  a high  speed  Tektronix  storage  scope,  model-7633, 
and  photographically  recorded. 

Power  measurements  were  made  for  the  full  range  of  lamp  drive 
currents  with  aperture  diameters  of  40  and  60  mils  at  2kHz  Q switch 
frequency.  Power  versus  Q switch  frequency  was  also  measured  at  18 
A lamp  drive  and  with  a 60  mil  aperture.  The  peak  power  per  pulse 
was  obtained  from  graphical  integration  of  the  normalized  to  unity 
photo  diode  signal  as : 

P . = P.„  T c 

peak  AV  (53) 

If (t)dt 


where  P 


AV 


T c 

f(t) 


average  power 
Q switch  cycle  time 
normalized  photo  diode  signal. 
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The  average  and  peak  powers  at  the  target  plane  as  a function  of 
lamp  current  at  2kHz  are  shown  in  Figure  127.  Curves  showing  the 
peak  and  average  power  as  a function  of  Q switch  frequency  are  shown 
in  Figure  128.  The  data  used  in  plotting  these  curves  is  also 
summarized  in  Table  VIII. 

TABLE  8 

POWER  SUMMARY  ND : YAG  LASER 


Lamp 

Current 

amps 

Q switch 
freq  - kHz 

Average  Power 
- watts 

Peak  Power 
kilowatts 

Laser 

Aperture 

12 

2 

.24 

.53 

60 

14 

fl 

.43 

1.175 

II 

16 

II 

.55 

1.52 

II 

18 

II 

.70 

1.95 

II 

18 

4 

1.03 

1.20 

II 

18 

6 

1.17 

.887 

II 

18 

8 

1.25 

.627 

II 

18 

10 

1.36 

.468 

II 

20 

2 

.83 

2.38 

II 

22 

II 

.96 

2.99 

II 

18. 

II 

.02 

.034 

40 

18.5 

II 

.04 

.097 

II 

20 

II 

.09 

.308 

II 

22 

II 

.20 

.966 

II 

The  photo  diode  output 

pulse  used  to  obtain  the  peak 

power  is 

shown 

for  the  two  operating  conditions  used 

in  the  laser 

diode  exper 

iments 

in  Figure  129  a and 

b.  Considerable 

variation  in 

the  output 

waveform  and  hence  the  output  power  is  noticeable  at  both  operating 
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Figure  129a.  Photo  diode  signal  - 20  A Lamp  drive  60  mil 
aperture  2kHz. 
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Figure  129b.  Photo  diode  signal  - 18  A Lamp  drive  40  mil 
aperture  2kHz. 
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points.  This  power  variation  is  believed  to  be  the  cause  of  the  laser 

4 

spot  variation  noted  previously  in  the  discussion  of  the  10  array 
experiments . 

This  power  data  will  be  used  in  the  near  future  to  get  a feel 

for  the  temperature  rise  required  for  diode  formation.  Before  this 

can  be  done  the  reflectance  of  the  sputtered  A1  film  must  be  known. 

However,  the  power  density  per  pulse  delivered  to  the  target  can  be 

determined  from  the  present  data.  Diodes  were  formed  with  20  A lamp 

current  and  a 60  mil  aperture  at  2kHz.  With  this  system  a 107o  beam 

splitter  was  also  used.  Power  was  measured  with  and  without  the 

beam  splitter  and  the  actual  transmission  was  found  to  be  6.677,. 

Therefore  using  this  value  and  the  peak  power  of  2383.35  watts  and 

a typical  irradiated  spot  size  of  15  ym  the  calculated  power  density 

per  pulse  is  8.99  x 10  watts /cm  . The  40  mil  aperture  18  A lamp 

current  operating  conditions  did  not  use  a beam  splitter  in  the  optical 

path.  Therefore,  with  a peak  power  of  33.87  watts  and  a typical  diode 

size  of  7.5  urn  the  calculated  power  density  per  pulse  is  7.67  x 10^ 

2 

watts /cm  . The  power  densities  used  to  form  the  diodes  under  different 
operating  conditions  are  therefore  in  reasonably  good  agreement. 

Lap  Stain  of  Written  Junctions 

In  order  to  further  investigate  the  diodes,  or  damaged  areas, 
formed  by  the  laser,  a sample  target  was  angle  lapped  and  stained. 
Figure  130  shows  a photomicrograph  of  the  lapped  action.  This  photo- 
graph reveals  that  the  aluminum  dissolved  in  the  silicon  (or  vice 
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Figure  130  Angle  Lapped  (5  40')  and  Stained 
Section  of  a Prototype  Archival 
Memory  on  a BLAMOS  Substrate. 
Patlern  is  10.jm  Spots  on  20um 
Centers . 


versa)  and  that  the  bit  pattern  went  to  a depth  of  5 to  6 microns. 
It  also  shows  that  the  "write-in"  was  nonuniform  across  the  wafer. 


This  method  investigation  probably  will  not  be  adequate  for  junctions 
of  0.2  micron  depth. 


Section  VII 


ALLOY  JUNCTION  EXPERIMENTAL 
WRITE  COLUMN 

Introduction 

Alloy  junction  surface  diodes  are  formed  by  locally  heating 
two  suitable  materials.  As  described  in  Section  VI,  junctions  of 
this  type  have  been  produced  using  an  imaged  laser  beam.  The  laser 
is  capable  of  delivering  far  more  energy  than  necessary  to  form  the 
diodes.  For  a memory  system,  however,  it  is  more  desirable  to 
use  an  electron  beam  for  writing,  since  it  can  be  deflected  more 
easily  over  the  desired  field  of  view.  Additionally,  an  electron 
beam  will  be  used  to  read-out  the  stored  information.  Part  of  this 
phase  of  the  advanced  Archival  Memory  Contract  will  determine  if  an 
electron  beam  system  can  be  structured  with  the  beam  characteristics 
necessary  for  alloy-j unction  writing.  As  part  of  these  studies,  and 
as  an  aide  in  determining  the  beam  requirements,  an  electron  beam 
writing  column  will  be  set-up  and  experimental  studies  conducted. 

Unfortunately,  the  analysis  to  determine  the  optimum  configur- 
ation for  such  a column  is  dependent  in  large  measure  on  the  stor- 
age target  studies.  The  laser  beam  experiments  will  be  helpful 
in  determining  power  levels,  but  the  manner  in  which  the  beam  energy 
heats  the  localized  areas  will  be  different  when  an  electron  beam 
is  employed.  For  this  reason,  an  electron  beam  system,  that  will 
come  as  close  as  possible  to  meeting  the  theoretical  and  analytical 


estimates  available  at  this  time,  will  be  assembled  and  electron 
beam  studies  conducted.  The  method  of  selection  and  specification 
of  such  a column  will  be  discussed  in  the  following  paragraphs. 


Column  Specifications 

The  amount  of  beam  current  necessary  to  heat  a local  volume 
of  the  alloy  junction  target  was  estimated  for  the  advanced  Archival 
Memory  proposal.  For  these  purposes  the  simplified  model  shown  in 
Figure  131  was  used.  The  beam  energy  is  assumed  to  be  completely 
dissipated  within  a hemisphere  of  radius  a given  by 

a = \ (b  + Rg)  (54) 

ii  it 

where  b,  is  the  beam  diameter  and  Rg  is  the  Grun  range.  The  Grun 
range  depends  on  the  material  and  is  a measure  of  electron  penetrat- 
ion into  the  material  which  is  proportional  to  beam  energy  raised 
to  the  1.75  power.  The  beam  requirements  to  locally  raise  the 
temperature  of  an  Indium  Germanium  system  30°C  above  a thermal  bias 
for  selective  alloying  have  been  determined  using  this  model.  The 
resulting  beam  current  and  voltage  as  a function  of  incident  beam 
spot  size  for  1000  and  2000  8 bit  sizes  are  shown  in  Figure  132. 

The  possibility  of  meeting  these  beam  requirements  with  an 
electron  column  as  shown  schematically  in  Figure  133,  has  also 
been  estimated.  A magnetic  final  lens  is  used  because  of  the  ext- 
remely short  focal  lengths  achievable.  Since  the  target  is  non- 
magnetic it  can  be  brought  close  to  the  lens  and  even  inside  the 
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Figure  132.  Predicted  Beam  Current  and  Beam  Voltage  Requirements 
for  Alloy  Junction  Writing  on  In-Ge  Target  at  0.1 
micron  and  0.2  micron  Bit  Spacing. 


BEAM  VOLTAGE  (kV) 


Figure  133.  Optics  for  Alloy  Junction  Experimental  Write  Station 


pole  pieces  without  distorting  the  focusing  field.  Hence  the  lens 
can  operate  with  an  extremely  short  focus  and  small  spherical 
aberration  which  results  in  larger  beam  currents.  However,  oper- 
ation of  the  final  lens  at  such  a short  focal  length  leaves  no  room 
for  deflection  between  the  lens  and  target.  Therefore,  a double 
deflector  preceding  the  final  lens  similar  to  that  used  in  SEM  op- 
tics is  used.  The  two  stage  deflector  is  arranged  so  that  the 
deflected  beam  always  passes  through  the  optical  center  of  the  lens. 
Final  deflection  over  the  field-of-view  results  from  varying  the 
angle  with  respect  to  the  optical  axis  of  the  principal  ray  as  it 
passes  through  lens  center. 

The  expected  beam  current  from  this  type  of  column  was  estimated 

. 2 

assuming  a thermionic  emission  Radley  Pierce  source  with  a 10  A/cm 
cathode  loading  and  a much  higher  birghtness  temperature  aided 
field  emission  (TFE)  source.  High  quality  magnetic  lenses  with  a 
spherical  aberration  coefficient  of  .13  cm  of  a type  previously 
built  at  General  Electric's  Research  and  Development  Center  were 
also  assumed.  The  results  of  this  preliminary  analysis  showing 
the  beam  current  capabilities  of  the  column  and  the  beam  current 
requirements  (from  Figure  132)  are  shown  in  Figure  134.  From  this 
figure  it  appears  as  though  a field  emission  type  source  is  required. 


Column  Selection 

Making  use  of  the  above  background,  the  specifications  for  a 
column  were  determined.  The  specifications  are: 
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Figure  134.  Predicted  Beam  Current  vs.  Spot  Size  for 
Experimental  Alloy  Junction  Write  Optics 
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1^  = beam  current  s 150  nA 
dg  = spot  size  s 1000  8 

= Landing  potential  = 3Kv  < s 5kV 

In  order  to  make  use  of  experimental  data  gained  from  writing 
experiments,  a decision  on  how  to  quickly  set-up  such  an  equipment 
was  necessary.  Several  possibilities  existed. 

• Design  and  construct  ourselves 

• Purchase  component  parts  from  various  vendors  and 
assemble  in  house. 

• Locate  a manufacturer  who  could  build  a system  capable 
of  meeting  the  specifications. 

This  latter  approach  was  selected  since  it  was  found  that  the  Coates 

and  Welter  Company  could  assemble  a complete  column  that  would 

meet  these  specifications.  After  obtaining  Contracting  Offices'* 

approval  for  the  purchase  of  the  equipment,  it  was  ordered  with 

delivery  promised  during  the  second  quarter. 

The  purchased  Coates  and  Welter  column  is  conceptually  similar 

to  the  column  of  Figure  133  as  shown  in  Figure  135.  The  column  has 

a field  emission  electron  source  housed  in  an  ultra  high  vacuum 

region  that  is  both  internally  and  externally  pumped  to  a pressure 
-10  -9 

of  10  to  10  Torr.  The  gun  chamber  is  differentially  pumped 
through  a conductance  limited  aperture  so  that  the  main  column 
pressure  can  be  approximately  1000  times  poorer  than  the  gun  chamber. 
The  gun  chamber  can  be  closed  off  with  an  isolation  valve  to  allow 
rapid  changing  of  target  samples.  The  first  lens  is  a specially 
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Figure  135.  Schematic  of  Purchased  Alloy  Junction  Write  Optics 
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designed  low  aberration  electrostatic  lens  and  the  final  lens  is 
a high  quality  magnetic  lens  as  assumed  in  the  preliminary  analysis. 
Magnetic  prelens  deflection,  similar  in  concept  to  the  electrostatic 
prelens  deflection  shown  in  the  column  of  Figure  133,  is  used. 
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Appendix  A 
ION  IMPLANTATION 

Ion  implantation  is  a technique  for  placing  atoms  in  the  surface 
layer  of  a solid  material  by  bombarding  the  substrate  with  energetic 
(keV  to  MeV)  ions . This  is  a particularly  valuable  tool  in  aemi- 
conductor  work,  because  a number  of  solid-state  properties  are  sen- 
sitive to  small  concentrations  of  foreign  atoms.  With  implantation, 
a wide  range  of  ions  can  be  introduced  in  specified  concentrations 
in  the  desired  location. 

Because  of  its  potential  in  device  fabrication,  there  is  consid- 
erable interest  in  the  ion  implantation  process.  This  process 
provides  an  attractive  alternative  to  thermal  diffusion  for  intro- 
ducing atoms  into  the  lattice.  Dopant  ions  can  be  implanted  at 
temperatures  considerably  lower  than  those  required  for  typical 
diffusions.  Impurity  concentrations  are  not  limited  by  the  usual 
thermodynamic  solubility  rules,  and  dopants  that  cannot  be  added  by 
conventional  means  can  be  used  with  implantation.  By  adjustment  of 
parameters  such  as  ion  energy  and  target  crystal  orientation,  the 
desired  profile  of  the  dopant  in  the  substrate  can  be  achieved  as 
well. 

Range  of  Implanted  Ions 

As  the  ion  being  implanted  penetrates  a solid,  it  loses  energy 
through  collisions  with  the  nuclei  and  electrons  in  the  target  mater- 
ial. The  actual  total  path  length  the  ion  travels  as  it  slows  down 


is  called  the  range  (R) , and  the  projection  of  the  range  along  the 
incident  direction  of  the  ion  beam  is  referred  to  as  the  projected 
range  (R^)  (Figure  A-l). 

The  amount  of  energy  lost  per  collision  and  the  number  of  collis- 
ions are  random  variables,  so  all  of  the  ions  in  a given  implant  do 
not  come  to  rest  in  the  same  spot.  Instead,  there  is  a distribution 
or  standard  deviation  in  the  ion  range,  which  can  be  calculated  and 
measured.  The  result  is  a Gaussian  shaped  profile  of  implanted  ions 
in  the  substrate  (Figure  A-2). 

For  an  amorphous  target,  this  range  distribution  depends  on  such 
parameters  as  the  ion  energy,  ion  mass  and  atomic  number,  target 
atomic  mass  and  atomic  number,  target  atomic  density,  target  temper- 
ature during  implantation,  and  ion  dose.  In  single  crystal  targets, 
surface  conditions  and  crystal  orientation  relative  to  the  ion  beam 
are  also  important  influences  on  the  range. 

As  the  ions  come  to  rest  they  lose  energy  by  two  major  mechanisms 
(Figure  A-3): 

• Inelastic  collisions  with  the  target  atoms,  which 
result  in  excitation  of  the  orbital  electrons  of  the 
ion  and  substrate  atoms. 

• Elastic  collisions,  involving  the  screened  nuclei  of 
the  ion  and  target  atoms. 

As  an  approximation,  these  processes  can  be  considered  independent 
of  each  other  for  calculation  of  energy  losses  and  the  ion  range. 
Thus,  for  an  incident  ion,  the  rate  at  which  energy  is  lost  with 
distance  can  be  given  as 
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Total  Path  Length  R Traversed  by  Implanted  Ion  and  Projected 
Range  R , Actual  Penetration  Depth  into  Substrate  (Source:  Ref.  24) 


2 


3 


X/( 

TARGET 

SURFACE 


Figure  A-2.  Gaussian  Distribution  Profile  for  Ions 

Implanted  Into  Amorphous  Target  (Source: 


Ref  .24) 


dE/dx 


Figure  A-3.  Nuclear  and  Electronic  Energy  Loss  dE/dx  as 
Function  of  Ion  Velocity  V (Source:  Ref. 25) 


af  ' N CS„  (E)  + S„  (E)J 


n 


(A-l) 


whe  re : 


N = Average  atomic  density 

S^(E)  = Nuclear  stopping  power 

S (E)=  Electronic  stopping  power 
e 


Integrating  Equation  A-l,  range  R for  an  ion  with  initial  energy  E 


can  be  written 


(Ref. 25)  (A-2) 


Nuclear  Stopping 

Considering  the  collisions  between  the  incoming  ion  and  substrate 
nuclei  as  elastic  interactions,  the  energy  loss  due  to  nuclear  stop- 
ping can  be  estimated.  Inelastic  collisions  resulting  from  very  high 
ion  energies  or  very  low  energies,  the  ones  that  involve  chemical 
binding,  are  not  relevant  to  most  implantation  conditions.  Thus, 
assuming  that  each  target  nucleus  acts  independently,  ions  are  slowed 
by  two  body  scattering  events,  as  shown  in  Figure  A-4.  The  energy 
loss  occurs  as  the  ion  travels  from  x to  x + Ax  through  the  target 
material,  passing  within  distance  p of  Nax  2irp  dp  target  nuclei, 
giving  up  an  energy  T(p^,E)  to  each  (Figure  A-5). 

The  energy  lost,  AE,  due  to  nuclear  interactions  can  then  be 
found  by  integrating  over  all  possible  impact  parameters 

(p,e)2^df>  (A.3) 
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Figure  A-5.  Model  for  Energy  Loss  T (p,E)  to  each 

Substrate  Atom  in  Annular  Ring  with  Volume 
27rp^dp  (Source:  Ref .25) 
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where  /\E  «E  (Ref  .25).  Dividing  by  Ax  and  taking  the  limit  as 
Ax  0,  the  average  rate  of  energy  loss  is 


g 


90  « 

n Crn  C r,£) 

So 


(A -4) 


= N S (E) 
n 


(4-5) 


For  this  type  of  scattering  interaction,  as  the  incident  ion 
passes  within  distance  b of  a target  atom,  it  is  deflected  through 
an  angle  0 as  it  gives  up  energy  T (b,E),  which  can  vary  from 


to 


T (p  = »,E)  = 0 
r (f*  0/  ' (M,  t rti.} 


(A -6) 

(Ref  .24)  ( A-7  ) 


Therefore,  knowing  T (p,E),  can  be  calculated.  The  solution 
of  this  elastic  collision  problem  is  obtained  using  classical 
mechanics  and  a change  of  coordinates,  switching  from  the  laboratory 
coordinates  of  Figure  A- 4 to  a system  in  which  the  origin  of  the 
coordinates  travels  with  the  center  of  mass  (Figure  A-6).  The 
incident  ion  energy,  E,  then  becomes 


E 

cm 


(A-8) 


where 


u = = reduced  ™ass. 

(Mx  + M2) 


(A-9) 


From  conservation  of  energy  and  angular  momentum  considerations 


E (l-p2u2  ) = V (u  ) 

cm  r max  max 


(A-10) 
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where  u = 1 /r,  and 


U *V  0* 


(A-ll) 


Using  this  basis,  the  center  of  the  mass  deflection  angle  can 


be  derived: 


* -*p 


t . z1*- 


(Ref.26)  (A- 12) 


For  a given  E , P,  and  interaction  potential,  V (u) , this  integral 
can  be  evaluated.  A general  form  for  this  potential  is 


V(r)  = 


ziz2  fM 

e r \a  / 


(A-13) 


where 


Z^  = Atomic  number  of  implanted  atom, 

Z ^ = Atomic  number  of  target  atom, 

and 

f^rj  = Screening  function  with  screening  parameter  a 

For  reasonably  accurate  numerical  calculations,  an  approximation  to 
The  Thomas-Fermi  potential  (Figure  A-7)  is  useful 


where 


a = 0.885  a (Z  2/3  + Z 2/3)"1/2 
o 1 2 


(A-14) 


8 


(r/a) 

Figure  A-7  Thomas -Fermi  Screening  Function 
Approximation  (Dashed  Line)  (Sou 
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and  C = /3  (Ref.27).  After  evaluation  of  the  integral,  0 can  be 


transformed  back  to  the  laboratory  frame  with 


tan  9, 


sin  Qcm 


cos  0 + M,  /M 

cm  12 


(A- 15) 


and  T calculated  using 


COS  6,  = 


1+  Ax.  d±\  JZ 

- C £ J e 


(A -16) 


Thus  T can  be  estimated  from  p using  the  sequency  p 0cm  9 T, 
and  therefore  S . Assuming  energy  loss  due  to  a series  of  small 
angle  scattering  processes,  9cm  is  small,  and  using  Equation  A- 15 
and  A -16 


S^— *T  2irp  dp 

(A -17) 

= T d-7 

(A -18) 

- V 

**  i /n 

(A -19) 

“ t1/2 


As  a crude  approximation,  the  nuclear  stopping  power  in  silicon 
can  be  given  in  the  form 


S = 2.8  x 10 
n 


-15 


Z1Z2 


M, 


,7  2/3  , _ 2/3. 1/2  , 

(Z1  + Z2  ) + M2 


eV  • cm  (Ref. 24)  (A-20) 


Electronic  Stopping 

Several  models  can  be  used  to  calculate  the  energy  lost  through 
electronic  interactions,  the  Lindhard,  Scharff,  and  Schi^tt  (LSS) 
and  Firsov  models  being  the  most  prominent. 


In  the  LSS  model  the  ion  is  considered  to  be  stopped  by  a free 
electron  gas  through  close  electron  collisions  and  plasma  resonance 


A^  = Atomic  weight  Ox  implanted  ion 

and 

A2  = Atomic  weight  of  substrate  atom  (Ref.27). 

The  Firsov  model  considers  the  overlap  of  electronic  orbitals 
during  the  collision  between  an  ion  and  a substrate  atom.  The 
transfer  of  electrons  from  one  to  the  other  results  in  a momentum 
transfer,  which  slows  the  ion  (Ref.28).  If  the  target  is  considered 
amorphous,  Sg  can  be  estimated  by  a scheme  rimilar  to  that  used  for 
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S , where  Tg,  the  energy  lost  due  to  electronic  interactions  in  a 
two  body  event  is  calculated  and  a random  distribution  function  is 


employed  to  result  in 

NS.Ce)  - "fn  (',«>  trrPJf 

o 


(A-26) 


Considering  the  situation  shown  in  Figure  A-8,  the  ion  and  target 
atom  form  a sort  of  molecule,  and  as  the  ion  comes  near,  the  outer 
less  tightly  bound  electrons  can  be  exchanged.  Electrons  going 
from  the  moving  ion  over  to  the  stationary  target  atom  lose  momen- 
tum, whereas  those  transferred  to  the  ion  must  gain,  with  these 
energy  exchanges  taking  a toll  from  the  ion's  energy.  Thus  as  the 
ion  travels  from  -®,  near  a traget  atom,  to  +»,  it  loses  an  amount 
of  energy 

where  v is  given  in  cm/sec  and  p is  given  in  cm.  providing  v<orbital 
electron  velocity.  Therefore  Sg  can  be  found  using  T in  Equation 

A-26 

S (V)  = 2.34  x 10'23  (Z.  + Z,)  v (eV  - cm)  (A -28) 

e i 4 

In  this  view,  all  of  the  electrons  are  considered  to  take  part  in 
energy  exchanges;  however,  in  actuality,  shielded  inner  electrons 
do  not  participate,  and  a reduced  atomic  number  should  be  used  to 
account  for  this  fact. 


\yt7.  i * ttTLz, t r 


(e  v) 


(A -2 7) 


s 


Figure  A-8.  Model  for  Electronic  Stopping,  Indicating 
Firsov  Quasimolecule  that  Forms  as  Elec- 
tronic Orbitals  Overlap  (Source:  Ref. 25) 
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Projected  Range  and  Distribution 

Knowing  and  Sg,  Equation  A-l  can  be  used  to  calculate  the 

range  or  total  path  length  traveled  by  the  ion.  However,  for 

practical  applications  this  is  not  as  relevant  a statistic  as  the 

actual  penetration  depth  into  the  substrate.  To  obtain  the  Gaussian 

distribution  shown  in  Figure  A-l  R and  AR  are  the  calculations 

P P 

of  interest.  Lindhard  and  his  coworkers  have  done  extensive  work 
on  this,  formulating  the  LSS  theory  for  atomic  stopping.  By  noting 
that  integration  of  the  probability  that  an  ion  has  a range  between 
R and  R + dR,  integration  from  0 to  » is  1,  analogous  integral 
equations  can  be  written  for  R^  and  other  lengths  indicated  in 
Figure  A-9,  which  can  be  converted  to  differential  equations.  Thus 
the  integral  equation  for  R^  is 


I 


. nf dr  [*T  <*>  - * ^'T>  cos 


( A-29) 


where 

da  = Differential  cross  section  for  energy  transfer  T 
0^  = Scattering  angle  in  laboratory  system 
N = Atomic  density. 

After  the  application  of  several  useful  approximations,  this 
equation  has  as  its  solution  in  terms  of  stopping  powers, 


e 

" f /, 


^ e*f  J J 


(A-30) 
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INCIDENT 
DIRECTION 
OF  AN  ION 


Figure  A-9.  Coordinate  System  for  Distances  R , R , and  Ri  for 
Ion  Entering  Substrate  at  0 and  cBming  to  Rest  at 
O'  (Source:  Ref .29) 


aL(E)  = M2  N Sn(E) 


M,  2 E 


(&)  - * 


<5V 

Ce)-frn1^ 


( A-31) 


(A-32) 


(A-33) 


2 2 2 2 

Under  the  same  energy  restriction,  (R^  = Rp  - 1/2R ± ) and 

2 2 2 2 

R (R  = R + Ri  ) have  solutions: 


i y \ . 2-Rr  Ce  ) d £ 

Rc  ^ / /V/>„  te')  + s<z  (-*03 


where 


= 


aj^CE) 


(A-35) 


( A-  3 6 ) 


and  x A*1 » o 1 /£r\ 

Thus  the  straggling  can  be  found  using 

CA*fp  *pGog&&-i*ty 


(A-37) 


( A-38 ~) 


Numerical  techniques  involving  the  use  of  a Thomas-Fermi  potential 


or  stopping  powers  with 


(A-39) 


where  S = constant  total  stopping  power  approximation,  enable 
n 
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calculation  of  the  straggling.  Table  A-l  indicates  the  results  of 
computer  calculations  of  and  for  various  ions  and  energies 
(Refs. 24  and 30) . 

As  a rough  extimate  for  intermediate  energies,  the  approximation 
of  Nielsen  is  sometimes  convenient  to  use  (Ref. 31) 

Cob  ^ rfe,%+z*y*- 


( A-40) 


where  g = density  (gm/cm  ): 


130  B C**  V)  I 

% ^ 


(A-41) 


For  ions  heavier  than  the  substrate  atoms,  an  approximate  value  for 
Rp  is  given  by 


1 + M2(3M1) 


(A-42) 


2 

/-l  ^ . 


(a-43) 


(A- 44) 


As  shown  in  Figure  A- 10,  the  more  energetic  an  ion,  the  deeper 
it  penetrates  the  substrate.  Similarly,  lighter  ions  have  longer 
ranges . 

The  maximum  in  the  ion  distribution  profile  thus  occurs  at 
distance  R^  into  the  substrate,  as  shown  in  Figure^-1,  with  a 


1 
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Table  A-l 


PROJECTED  RANGE  (Rp)  AND  STANDARD  DEVIATION  IN  Rp  (ARp)  IN  X 
FOR  COMMON  IONS  IMPLANTED  IN  SILICON  AT  VARIOUS  ENERGIES 

(Source:  Ref.  24) 


Projected  Range, 

R 

Energy 

(keV) 

ARp  (X) 

20 

40 

60 

80 

100 

120 

140 

160 

180 

200 

Boron 

wm 

1413 

2074 

2695 

3275 

4289 

4745 

5177 

5588 

mm 

443 

562 

653 

726 

855 

910 

959 

1004 

Nitrogen 

491 

961 

1414 

1847 

2260 

2655 

3034 

3391 

3728 

4046 

191 

312 

406 

479 

540 

590 

633 

672 

710 

745 

Aluminum 

289 

564 

859 

1141 

ESI 

■tit! 

2036 

2335 

2633 

2929 

107 

192 

271 

344 

■si 

mm 

535 

591 

644 

693 

Phosphorus 

255 

488 

729 

976 

1228 

1483 

1740 

1998 

2256 

2514 

90 

161 

228 

291 

350 

405 

459 

509 

557 

603 

Gallium 

155 

272 

383 

492 

602 

712 

823 

936 

1049 

1163 

37 

64 

88 

111 

133 

155 

176 

197 

218 

238 

Arsenic 

151 

263 

368 

471 

574 

677 

781 

885 

991 

1097 

34 

59 

81 

101 

122 

141 

161 

180 

198 

217 

Indium 

133 

223 

304 

381 

456 

529 

601 

673 

744 

815 

23 

38 

51 

63 

75 

86 

97 

108 

119 

129 

Antimony 

132 

221 

300 

376 

448 

519 

590 

659 

728 

797 

22 

36 

49 

60 

71 

82 

92 

102 

112 

122 
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Figure  A-10  Ion  Range  Profiles  Indicating  Variation  of  Shape  and  Depth 
of  Distribution  with  Energy  (Source:  Ref.  32) 


r 


concentration 


(atoms/cm^) 

2toRp) 


(A-45) 


where  = ion  fluence  (atoms /cm  ) and  /\R^  is  given  in  centimeters 
(Ref. 24  ). 

Lateral  Spreading 

The  lateral  spread  of  implanted  ions  becomes  of  particular 
concern  as  interest  in  very  well  defined  doped  regions  grows.  In 
addition  to  the  standard  deviation  of  the  projected  range,  the 

lateral  standard  deviation  can  be  calculated  similarly.  Noting  from 


Figure  A -9  that 


2 2 2 
R - X + Y 
P 

~~2  ~~2 
X « 


for  an  amorphous  solid  substrate  the  lateral  standard  deviation  is 


(A-46) 


Assuming  a Gaussian  distribution  in  three  dimensions,  the  spatial 
probability  distribution  for  the  point  where  an  implanted  ion  comes 


to  rest  is  given  by  (Ref .32) 


. . . f /Vy1- 


(Ref.  32)  (A-47) 


Using  the  LSS  computer  calculation  techniques  for  boton  implanted  in 
silicon,  for  example,  the  curves  of  Figure  A-ll  are  generated. 
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INCIDENT  ENERGY 


Figure  A-ll.  Standard  Deviation  ^ R in  Projected  Range 

and  Lateral  Standard  ^Deviation  A±  as  Functions 
of  Ion  Energy  (Source:  Ref. 29). 


Particularly  for  low- implantation  energies,  the  predicted  lateral 
scatter  appears  small. 

The  few  experimental  studies  of  lateral  spreading  support 
standard  deviations  predicted  in  this  manner.  Figure  A-12  presents 
the  results  observed  for  boron  implanted  in  silicon,  and  Table  A-2 
indicates  both  theoretical  and  experimental  work  for  Kr  and  Ar  into 
Si. 

Channeling 

The  preceding  predictions  of  the  implanted  ion's  ranges  and 
straggling  were  based  on  the  premise  that  the  substrate  lattice 
appears  amorphous  to  the  incoming  ion.  However,  in  a single-crystal 
lattice  resulting  in  deep  ion  penetration  and  less  radiation  damage. 
Providing  the  ion  beam  enters  the  crystal  within  a critical  angle 
of  a crystallographic  axis,  the  repulsive  force  between  the  screened 
Coulomb  fields  of  the  ions  and  substrate  atoms  acts  to  steer  the 
ions  between  the  atomic  rows  so  that  thev  lose  little  energy  to 
nuclear  interactions  (Figure  A-13) . 

This  critical  angle,  *c,  can  be  estimated  by  considering  the 
ion's  incident  energy,  E,  to  be  composed  of  two  components,  E,j  and 
E , parallel  and  perpendicular  to  the  atomic  rows  (Figure  A-14). 

If  E is  less  than  the  repulsive  potential  energy  from  the  substrate 

i 

atoms,  the  ion  undergoes  channeling.  Approximating  the  configuration 
to  a line  of  charge,  the  shielded  potential  can  be  written  using  the 
screening  function  of  Equation  A-13.  Setting  U (rm^n  ) = E = E sin 
h allows  the  minimum  approach  distance  to  be  calculated.  Because 
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Figure  A-12.  Lateral  Spread  of  Boron  Implanted  in  Silicon 
as  a Function  of  Energy  (Source:  Ref.  33) 


I 


-2 


STANDARD  DEVIATION  AND  LATERAL  SPREAD  VALUES 
FROM  THEORY  AND  EXPERIMENTS 
(Source:  Ref.  29) 


js  'ai 


Figure  A-13.  Steering  of  Implanted  Ion  Along  Channel  Formed  by 

Rows  of  Substrate  Atoms  with  Spacing  d(Source:  Ref . 24) 


Figure  A-14.  Division  of  Energy  into  Components  for  Calculation 
of  Critical  Angle  h , which  Determines  Whether 
Channeling  Occurs  ( Source:  Ref . 24) 
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It  (<■)  = 2.  Zf Jl 

UL) 


- ■ J 


(A-48) 


an  atom  could  penetrate  the  channel  walls,  since  they  are  not  really 
continuous  walls  of  charge,  a requirement  that  the  ion  pass  many 
substrate  ions  while  traveling  near  the  middle  of  a channel  is 
imposed  as  an  additional  condition  for  channeling: 


min  > d 


(A-49) 


where  v„  and  v are  parallel  and  perpendicular  ion  velocity  compon- 
M 1 


ents . Thus : 


>»  - ^ TT-  ' ^ ~tanf 


(A-50) 


for  small  tlj . 


Using  this  r , , the  critical  angle  can  be  estimated, 
min 


e > ed  = (eS) 

° 60  a 


(A-51) 


then : 


with 


(Refs  .24  and  27)  (A-52) 


(A-53) 


However,  if  E<:Eo  (Refs . 24  and27) 
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K. 


(A- 54) 


Table  A-3  indicates  critical  angles  for  channeling  for  several 
common  ions  implanted  at  different  energies  along  the  major  crys- 
tallographic axes.  The  fraction  that  are  channeled  can  be  estimated, 
as  well,  by  the  model  that  ions  approaching  within  r ^ of  a sub- 
strate atom  are  scattered  out  of  the  channel.  Assuming  a surface 
density,  Nd,  in  the  substrate,  the  fraction  of  the  random  beam 
(portion  scattered  out)  is 

9 -9 

(A-55) 


f , - t (r  . )^N  d(cm  ^) 

random  v min' 


However, 


r . 
mxn 


d ii/C 


(A-56) 


for  small  i!jCj  and  thus  the  channeled  fraction  is  (1  - f , ). 

random 

Although  the  range  distribution  for  the  random  fraction  may  be 
calculated,  the  path  length  traversed  by  a channeled  ion  is  not  as 
predictable.  Due  to  variable  influences  of  increasing  dose  and 
surface  disorder  (making  the  lattice  more  amorphous)  and  temperature 
(large  vibrational  amplitudes  of  the  substrate  ions  reduce  channel 
size),  only  the  maximum  penetration,  R[iax,  that  a perfectly  chan- 
neled ion  travels  is  a reproducible  calculable  number  (Figure  A-15). 
Thus,  considering  the  ion  beam  to  separate  into  independent  random 
and  channeled  components,  the  resulting  distribution  profile  of 


~1 
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Table  A- 3 

CRITICAL  ANGLES  FOR  CHANNELING  NEAR 
MAJOR  CRYSTALLOGRAPHIC  AXES  OF  SILICON 
FOR  SEVERAL  ENERGIES 
(Source:  Ref.  27) 


Ion 

Energy 

(keV) 

Channel  Direction  (degrees) 

<110? 

<111? 

<100? 

Boron 

30 

4.2 

3.5 

3.3 

50 

3.7 

3.2 

2.9 

Nitrogen 

30 

4.5 

3.8 

3.5 

50 

4.0 

3.4 

3.0 

Phosphorus 

30 

5.2 

4.3 

4.0 

50 

4.5 

3.8 

3.5 

Arsenic 

30 

5.9 

5.0 

4.5 

50 

5.2 

: 

4.4 

4.0 
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DEPTH  (MICRONS  OF  SILICON) 


Figure  A-15 


Distribution  Profiles  of  p+  Ions  Implanted  in  Silicon 
at  40  keV  Along  <110>  Axis  showing  Variation  with: 

(a)  Substrate  Temperature  and  (b)  Dose  (Source:  Ref.  34) 
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implanted  ions  is  a superposition  of  statistics  for  each.  This 
results  in  a double  peaked  distribution,  as  shown  in  Figure  A_i6, 


with  a fraction,  f^^om,  stopped  near  and  (1  - francjon,)  stop- 
ping near  Rmx,  with  the  area  between  peaks  filled  in  by  other 
dechanneled  ions  (Refs. 24  and  36). 

Perfectly  channeled  ion  energy  loss  through  nuclear  processes 
is  negligible;  therefore,  energy  loss  may  be  described  by 

dE  « S (E)  (A-57) 

dx  c 


Using  the  Firsov  formula 
T » 


a Sjr 

— ft  / V ft- 

‘t.S-fifO  (X, 


fj  + 3,/*  /o'7  (z , tZ+yj  pjr 

for  energy  T lost  in  each  collision  and  substituting  d*c  for  p and 

reduce  atomic  numbers  Z, ’ and  Z S can  be  estimated.  However, 

1 Z e 

due  to  oscillations  in  with  Z^  (Figure  A -17),  this  method  is 
difficult.  A more  accurate  technique  utilizes  measured  S e for  the 
desired  ion  target  and  orientation  so  that : 

J E . i r*  1 /—A  _ ^ er 


■ • % (.*& 

(A-58) 

(A-59) 

N-s.'  (.?.*) 

l£o*)^ 

(A-60) 

where 


can  be  found  from  measurements  performed  at  energy  Eq  . Thus  the 

maximum  range  of  a channeled  ion  can  be  fond  using 

2 (^o+)  /Bo  Y^- 
° ^ - — ' — (A-61) 


M5ft' 
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Concentration 


Figure  A-16  Distribution  of  Implanted  Ions  in  Single-Crystal  with  Ion  Beam 
Aligned  with  Major  Crystallographic  Axis  (Source:  Ref.  35) 


eV  cm 


Figure  A- 18  gives  Rmflv  for  various  implantation  conditions: 
however,  the  influence  of  crystallographic  axis  direction  (varying 
channel  size)  is  more  graphically  displayed  in  Figure  A-19.  Even 
\rtien  the  target  crystal  is  deliberately  misaligned  to  reduce  the 
amount  of  channeling,  some  channeling  still  may  occur,  causing  the 
formation  of  tails  in  the  range  distribution  as  ions  entering  chan- 
nels travel  further  than  R up  to  R 

p max 

Radiation  Damage 

A limitation  to  the  value  of  ion  implantation  as  a doping 
technique  is  the  problem  of  lattice  disorder  resulting  from  atomic 
displacements  due  predominately  to  elastic  collisions  during  the 
implantation  process.  For  a crystalline  target,  the  substrate 
atoms  are  bound  to  their  lattice  sites  by  a displacement  energy. 

As  the  implanted  ion  travels  through  the  crystal,  it  makes  many 
violent  collisions  with  atoms  in  the  substrate.  A displaced  atom 
may  act  as  a secondary  projectile  if  struck  by  an  ion  with  energy 
greater  than  its  displacement  energy.  Such  secondary  projectiles 
come  to  rest  through  elastic  and  inelastic  collisions,  producing 
another  generation  of  projectiles  if  they  are  sufficiently  energetic, 
and  so  on.  Thus  the  damage  production  is  a cascade  process,  with 
the  incident  ion  the  primary  disorder  causing  projectile.  The  extent 
and  character  of  the  disorder  zone  or  cluster  are  governed  by  the 
mass  of  implanted  ion  and  its  energy,  the  mass  of  the  substrate  atom, 
the  ion  fluence,  and  the  substrate  temperature. 


suojojuj  xv|ft,y 


Figure 


Energy  (keV) 


-18  Experimental  Measurements  of  R x for  P Ions  in 

Silicon  as  a Function  of  Energy  (After:  Ref.  38) 


For  relatively  heavy  ions  such  as  phosphorus  implanted  in 
silicon,  a large  fraction  of  the  energy  goes  into  elastic  collisions, 
resulting  in  more  lattice  damage.  There  is  a greater  tendency  for 
the  disorder  to  occur  within  a damage  core,  surrounded  by  point 
defects  (Figure  A-20).  These  point  defects  can  diffuse  away  from 
the  core.  A considerable  portion  of  the  radiation  damage  stable 
at  room  temperature  is  believed  to  result  not  directly  from  the 
incident  ion  impinging,  but  instead  from  this  migration  of  defects 
to  nucleation  centers. 

Lighter  ions  such  as  boron  lose  a large  part  of  their  energy 
into  inelastic,  electronic  processes.  Most  of  the  displaced  atoms 
form  isolated  defects  along  the  ion  track  (Figure  A-21),  rather  than 
amorphous  zones.  Of  course  the  substrate  atomic  mass  is  also  an 
important  parameter,  determining  the  degree  of  elasticity  in  the 
collisions,  depending  on  whether  it  / le'.  s or  more  than  the  incom- 
ing ion's  mass. 

As  the  fluence  or  dose  of  ions  is  increased,  the  number  of 
damage  zones  increases  linearly  until  they  apparently  begin  to  over- 
lap, resulting  in  a continuous  amorphous  layer  (Figure  A-22) . How- 
ever, the  substrate  temperature  has  a marked  effect  on  the  amount  of 
radiation  damage  produced  during  the  implant  process  (Figure  A-23). 

At  low  temperatures,  the  damage  is  frozen  in  as  it  occurs,  and  anneal- 
ing of  defects  proceeds  very  slowly.  With  elevated  temperatures, 
individual  disordered  zones  can  anneal,  preventing  the  formation  of 
disorder  and  resulting  in  less  overall  damage. 
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Figure  A-20.  Damage  Cluster  Formed  by  Heavy  Ion  Implantation 
(Source:  Ref.  39) 


Figure  A-21  Damage  Track  Produced  by  Light  Ion  Implantation 
(Source:  Ref.  39) 


DISORDER  (RELATIVE  UNITS) 


10,z  1 o13  1014  10'5  10'6  1 017 


Sb  DOSE  (ions/cm2) 


Figure  Ik-22. 


Increase  in  Lattic^  Disorder  as  a Function  of 
Dose  for  40-keV  Sb  Ions  Implanted  in  Silicon 
(Source:  Ref.  40) 


-290- 


i 


A variety  of  defects  may  be  generated  in  the  implantation  proc- 
ess. These  defects  have  been  studied  using  tools  such  as  EPR,  IR 
absorption,  luminescence,  channeling  and  backscattering,  transmission 
electron  microscopy,  and  electrical  measurements.  Through  such 
experiments,  the  nature  and  behavior  of  defects  can  be  characterized 
to  some  extent,  allowing  tentative  predictions  of  the  percentage 
of  electrically  active  implanted  sites  and  annealing  needed  to 
maximize  the  desired  characteristics. 

A Frenkel  defect  is  formed  when  a silicon  atom  is  displaced 
from  its  lattice  site  to  an  interstitial  position.  The  resulting 
vacancy  may  occur  in  various  charge  states  (+,  0,  and  =)  at 
different  energy  levels  in  the  bandgap.  Free  vacancies  anneal  in 
the  range  from  70°  to  150°K,  depending  on  the  charge  state.  Because 
they  are  highly  mobile  in  silicon,  vacancies  can  migrate  through  the 
crystal  and  become  trapped  by  other  vacancies,  impurities,  or  defect 
complexes  (Figure  A-24) . 

Interstitial  silicon  atoms,  displaced  from  their  sites  in  the 
substrate  lattice,  are  also  believed  to  migrate  at  very  low  temper- 
atures, as  low  as  4°K.  Moving  through  the  lattice,  they  may  dis- 
place selected  impurities  from  their  substitutional  sites,  resulting 
in  the  formation  of  an  interstitial  impurity  and  a substitutional 
silicon,  a conversion  commonly  called  the  "Watkins  Replacement 
Mechanism"  (Figure  A-25)(Ref.  43).  The  interstitial  silicon  may 
also  be  active  in  forming  complexes  with  impurities  and  other  defects 
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Figure  A-24.  Annealing  Curves  for  Vacancy  Complexes  in 

Radiation  Damaged  Silicon  (Source:  Ref.  42) 


mation  of  Vacancy-Interstitial  Pair  (Frenkel  Defe 
sequent  Watkins  Replacement  Mechanism.  (Source: 


The  divacancy  is  particularly  important  in  room  temperature 
implanted  silicon,  because  it  is  stable  up  to  about  550°K.  Divac- 
ancies may  be  formed  directly  when  an  energetic  ion  displaces  ad- 
joining silicon  atoms  or  when  vacancies  of  the  same  charge  state 
combine  in  pairs.  Like  vacancies,  divacancies  appear  in  several 
charge  states . 

In  addition,  a number  of  more  complicated  defects  are  present 
in  implanted  material.  Many  of  these  defects  have  been  studied 
through  EPR  and  other  measurements  and  identified  as  clustered 
defects  and  multiple  vacancy  complexes. 

The  number  of  point  defects  introduced  by  the  implantation 
process  can  be  estimated  by  modeling  the  situation  after  hard 
sphere  collisions,  providing  that  the  incoming  ion  is  not  moving 
too  fact,  the  limiting  energy,  EA,  for  this  assumption  being 


£*  * 2-  % z,  z,_  Cz, 


(A-62) 


where  E^  = 13.6  eV.  This  is  the  Coulomb  interaction  potential  of 
two  screened  nuclei.  If  the  incident  ion  has  E<EA>  then  an  estimate 
for  the  number  of  substrate  atoms  displaced  by  primary  recoils  is 
given  by 


(A-63) 


where 

E = Energy  of  primary  recoil 
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E,  = Displacement  energy  of  substrate  atom  (about  15  eV  for 
silicon) 


If  E>Ea,  not  all  of  the  nuclear  interactions  cause  atomic  dis- 
placements. For  such  an  energy.  Eg  or  greater 


5.-  Y 


+ y-EQ. 


(A-64) 


where  half  of  the  nuclear  collision  energy  goes  into  lattice  vib- 
rations instead  of  displacements,  the  number  of  displaced  atoms  is 
approximate ly 


N. 


(A-65) 


In  many  situations  the  electronic  stopping  power  dominates, 
E>E^ . Defining  Ec  as  the  energy  for  which  the  electronic  energy 
loss  is  greater  than  the  nuclear  loss  and  Eq,  the  ion  energy,  is 
»E„  and  ED,  the  total  number  of  displaced  atoms  can  be  estimated 

L.  D 


by 


(A-66) 


where  P ~ 10 


and 


b = 1/2  Eg  > Ec  (A-67) 

b = 1/4  Eg  < Ec  (Ref.  45)  (A-68) 

and  E„  is  the  energy  for  S = S . 

C n e 

Such  calculations,  for  example,  indicate  that  a heavy  ion  such 
as  antimony  will  at  100  keV  produce  about  3000  displacements.  A 
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much  lighter  boron  ion  of  the  same  energy,  which  loses  most  of  its 
energy  to  electronic  scattering  processes,  produces  only  700  dis- 
placements, most  of  which  occur  near  the  end  of  the  ion  track. 

The  distribution  of  defects  penetrates,  in  general,  slightly 
deeper  than  the  actual  implanted  ion  distribution  (Figure  A-26) . 
Using  an  approximation  to  the  Thomas-Fermi  potential,  numerical 
values  for  the  spatial  distribution  of  radiation  damage  can  be  cal- 
culated, giving  results  such  as  those  of  Figure  A-27.  These  resits 
exhibit  the  ratio  between  the  projected  range  and  depth  of  the 
damage  and  the  straggling  in  these  in  relation  to  the  ratio  of  the 
ion  to  the  substrate  mass.  For  relatively  light  ions,  which  can 
undergo  wide  angle  scattering,  the  peaks  for  the  range  and  defect 
distributions  occur  at  nearly  the  same  depth  and  the  distribution 
shapes  are  similar,  except  the  damage  exhibits  less  straggling.  If 
the  ion  is  heavy  and  moving  slowly,  the  ion's  penetration  depth  can 
be  half  the  disorder  depth,  since  little  wide  angle  scattering  can 
occur  and  the  range  straggling  is  slight.  As  shown  in  Figure  A-20, 
the  damage  for  this  case  is  spread  along  the  extent  of  the  ion  path. 

These  defects  in  the  crystalline  structure  can  be  healed  through 
annealing  for  selected  times  and  temperatures.  Such  annealing  may 
also  have  the  added  benefit  of  increasing  the  number  of  electrically 
active  implanted  sites  in  the  substrate.  While  some  defects  anneal 
at  or  near  room  temperature  or  below,  in  general  the  silicon  lattice 
damaged  by  ion  implantation  starts  to  reorder  above  200°C.  If  the 
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substrate  has  been  driven  amorphous,  temperatures  above  550°C  are 
needed  for  reordering  to  begin  (Figure  A-28). 

The  amorphous  layer  apparently  undergoes  epitaxial  regrowth  on 
the  silicon  substrate  underneath,  producing  a highly  ordered  elec- 
trically active  layer  after  annealing.  To  achieve  the  same  degree 
of  crystalling  perfection,  a lightly  damaged  substrate  would  require 
considerably  higher  temperatures  (~  900  C) . With  high-temperature 
annealing  dislocation  loops  and  lines  are  observed,  which  increase 
in  size  with  temperature.  The  precise  annealing  schedule  followed 
is  dependent  on  the  ion  species  and  substrate  being  used,  in  order 
to  obtain  maximum  electrical  characteristics  from  a particular  im- 
plant. 

Ion  Sites  and  Electrical  Activity 

Most  implantations  in  semiconductors  are  directed  toward  placing 
dopant  ions  at  particular  sites  in  the  substrate  crystal  lattice, 
usually  substitutional  for  silicon.  Using  the  channeling  phenomena 
as  a tool,  it  is  possible  to  estimate  the  numbers  of  implanted  ions 
occupying  substitutional  or  interstitial  sites  by  studying  the 
directional  character  of  very  light  ions  scattered  back  by  heavier 
ions  already  implanted  in  the  material,  as  shown  in  Figure  A-29. 

The  proportions  of  substitutionals  and  interstitials  are  highly 
dependent  on  a variety  of  parameters,  such  as  ion  species,  implant 
temperature,  dose,  and  annealing.  Detailed  information  allowing 
accurate  estimations,  however,  is  not  available  in  many  cases.  Table 
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Figure  A- 29.  Configuration  Utilizing  Channeling  to  Determine 
Impurity  atom  Locations  (Source;  Ref.  35) 


A-4  indicates  some  results  for  hot  implants  in  silicon. 

Apparently  some  correlation  exists  between  the  number  of  sub- 
stitutional sites  and  the  chemical  as  well  as  electrical  behavior 
of  the  implanted  ions,  particularly  if  defect  impurity  interactions 
are  influential.  The  perfect  substitutional  species,  in  general, 
could  be  expected  to  exhibit  a radius  near  that  of  a substrate  atom, 
a slight  variation  in  electronegativity  from  the  host,  and  a valence 
near  that  of  the  target  atom.  An  ion  that  is  stopped  in  the  substrate 
has  to  compete  for  lattice  sites.  If  the  energy  needfed  to  form  an 
interstitial  impurity  is  greater  than  that  needed  for  an  inter- 
stitial substrate  atom,  the  ion  can  move  onto  an  interstitial  site. 

Due  to  the  lesser  valence  electron  ionization  of  group  III  atoms, 
silicon  interstitials  may  compete  for  and  win  the  available  substit- 
utional sites  over  some  impurities.  As  the  ionic  radius  increases, 
this  phenomenon  is  less  marked.  Because  group  V atoms  need  more 
ionization  than  group  IV  substrate  atoms,  these  impurities  tend  to 
take  up  substitutional  sites  (Ref.  47  ).  However,  over  a particular 
temperature  range,  impurities  may  move  on  or  off  sites  and  then 
sometimes  return,  depending  on  the  particular  species. 

Using  Hall  effect  measurements,  a more  exact  determination  of  the 
number  of  active  ions  in  the  implanted  layer  than  estimations  by 
impurity  locations  can  be  made.  The  activity  level  depends  not 
only  on  the  ion  species  and  substrate  temperature  during  implantat- 
ion, but  also  on  the  anneal  time,  dose,  and  initial  substrate  resis- 
tivity. As  the  dose  is  increased,  higher  annealing  temperatures  are 
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LATTICE  LOCATIONS  OF  IONS  IMPLANTED  AT  350°C  INTO  SILICON 
(I  = Interstitial,  S = Substitutional  (Source:  Ref.  49)) 
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needed,  in  general,  to  achieve  maximum  activity,  unless  the  sub- 
strate is  driven  amorphous.  As  the  amorphous  layer  regrows  through 
an  epitaxial  process,  high  electrical  activity  is  possible  as  the 
damage  anneals  at  lower  temperatures.  Increasing  the  implantation 
temperature  often  increases  the  temperature  requirements  for  maxim- 
izing activity  as  well  as  reducing  radiation  damage. 

Implanted  layers  usually  show  increased  electrical  activity  for 
anneals  above  300°C  as  the  defects  around  the  ion  diffuse.  However, 
annealing  is  not  always  necessary,  because  boron,  for  example, 
produces  type  conversion  immediately  after  room  temperature  implan- 
tation. For  most  ion  species,  900°C  annealing  maximizes  the  carrier 
density.  As  more  implantation  phenomena  are  studied,  the  mechanisms 
involved  in  determining  the  activity  level  of  implanted  ions  will 
be  better  understood,  and  patterns  in  maximizing  the  number  of 


carriers  can  be  discerned. 


Appendix  B 

ELECTRON  PENETRATION  INTO  SOLIDS 
(Summary  of  Experimental  and  Theoretical  Understanding) 

This  appendix  provides  a brief  overview  of  the  present  under- 
standing of  electron  penetration  into  solids  (and  gases),  with 
emphasis  on  lateral  scattering  and  energy  deposition  profiles  at 
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beam  energies  less  than  10  keV  in  low-z  materials. 

Historically,  much  of  the  early  theoretical  work  in  this  area 
was  motivated  by  beta  ray  experiments,  where  the  electron  energies 
are  approximately  1 millen  eV,  and  in  more  recent  times  by  trans- 
mission electron  microscopy  in  the  energy  range  of  30  to  100  keV. 

Interest  in  energies  below  10  keV  was  primarily  motivated  by  a 
desire  to  check  certain  aspects  of  the  higher  energy  theories  over 
a wider  energy  range.  More  recently,  interest  in  electron  beam 
lithography  has  generated  renewed  interest  in  this  energy  range. 

Most  experimental  measurements  have  concentrated  on  secondary 
emission,  backscattering,  and  transmission  through  thin  foils.  Measure- 
ments of  lateral  scattering  are  rare,  because  of  their  obvious  diffic- 
ulty. 

When  electrons  penetrate  a solid  target,  they  can  be  scattered 
either  elastically  or  inelastically.  Elastic  collisions  are  primar- 
ily due  to  the  screened  nuclear  field  and  can  result  in  large  angle 
(>  90  degrees)  deflections,  although  most  are  smaller  angle  deflec- 
tions. Inelastic  collisions  are  primarily  electron-electron  events 
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that  result  in  a substantial  energy  loss  (~100  eV)  but  smaller 
angular  deflection,  because  of  the  small  electron  mass.  Hence 
the  electron  penetration  process  can  be  modeled  as  independent, 
large  angle,  elastic  scattering  and  zero  angle  inelastic  scattering. 

Inelastic  collisions  are  generally  treated  in  the  continuous 
energy  loss  approximation  (Ref. 50).  It  was  shown  that  the  rate  of 
energy  loss  per  unit  of  actual  path  length,  s,  is  given  by  (Ref.  51): 

_ (ae\ 

Ts'  LfiJ  ll/(B-l) 


This  formula  is  essentially  exact  for  hydrogen,  where  the  elec- 


tron wave  functions  are  known.  For  more  complex  atoms,  it  is  only 
a generalization.  NA  is  Avagadro’s  number,  o is  the  mass  density, 

Z the  atomic  number,  and  A is  the  atomic  weight.  I is  the  mean 
excitation  energy  ( I ~11.5  Z eV) . The  total  path  length,  S, 
traversed  by  an  electron,  slowing  down  from  energy  to  E2  is  then 


given  by: 
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( B-2) 


A useful  approximation  to  s is  the  Webster  law  (Ref. 52); 
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where  0 is  given  in  g/cm  , E in  keV,  and  s in  cm.  Slight  variations 
on  this  relation  are  common  to  most  theoretical  treatments. 


The  elastic  scattering  is  generally  treated  as  Rutherford 
scattering  from  a screened  Coulomb  field.  An  example  of  such  a 
calculation  is  due  to  Wentzel  (Ref . 53) ; in  which  Cosslett  and  Thomas 
assumed  an  exponential  screening  field.  They  found  the  differential 
cross  section/solid  angle  to  be: 


( B-4) 

where0o  = 3.69  x 10‘2Zl/3E~1//2 , with  E in  keV.  For  Z = 10,  E = 10 
keV  and@Q  = 0.07  radian. 

Many  theoretical  models  for  electron  penetration  have  been  built 
on  variations  of  these  two  equations.  Among  them  are  References  50, 
52 , and  55  through  69  • 

Most  of  these  studies  have  attempted  to  describe  x-ray  production 
as  a function  of  depth,  secondary  emission,  and  backscattered  elec- 
trons. No  attempt  has  been  made  to  predict  lateral  scattering, 
although  many  models  contain  lateral  scattering  information.  In 
these  types  of  models,  the  electrons  penetrate  depth  x^  into  the 
material,  losing  energy  according  to  the  Bet he  law  and  scattering 
according  to  a screened  Coulomb  cross  section.  At  depth  x^,  the 
scattering  has  essentially  randomized  the  directions  of  the  elec- 
trons. From  this  depth  the  electrons  diffuse  into  a sphere,  still 
losing  energy  according  to  the  Bethe  law.  The  treatment  of  the 
diffusion  is , of  course,  difficult  because  the  scattering  events 
are  not  random  in  angle,  so  a Boltzmann  equation  is  required.  This 
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type  of  approach  has  been  used  by  Bethe,  Rose,  and  Smith  (Ref.  60), 
Brown  and  Ogilvie  (Ref.  67),  and  Bennett  and  Roth  (Ref.  68)  to  study 
dE/dX,  X-ray  production,  and  secondary  electron  emission.  In  all 
cases,  only  the  variation  with  penetration  depth  was  studied;  how- 
ever, the  models  contain  lateral  scattering  information  if  it  can 
be  solved.  The  division  by  Archard  into  penetration  to  diffusion 
depth  x^  and  then  radial  diffusion  from  this  point  seems  especially 
promising,  because  then  the  Boltzmann  equation  has  radial  symmetry. 

More  sophisticated  recent  attempts  to  improve  on  this  type  of 
approach  are  the  plural  scattering  model  of  Greeneich  and  Van  Duzer 
(Ref.  70)  and  the  multiple  scattering  model  of  Hawryluk  et^  al_.  (Ref. 
71) . These  models  have  been  solved  for  lateral  scattering  information. 
The  work  to  date  has  been  restricted  to  beam  energies  at  or  above 
10  keV. 

Monte  Carlo  calculations  have  also  been  applied  with  some  success. 
Among  these  are  the  calculations  of  Lifshin  and  Bolon  (Ref.  72)  of 
Corporate  Research  and  Development  and  Hawryluk  et  al.  (Ref.  71) . 
Again  the  problem  is  that  these  calculations  have  not  been  applied 
below  10  keV.  Bolon  feels  that  the  Monte  Carlo  programs  as 
presently  written  may  not  be  valid  below  10  or  15  keV,  because  of 
the  simplified  scattering  laws  used. 

Experimental  Knowledge:  Energy  Deposition  Profiles  and  Electron  Range 

The  dE/dx  is  defined  as  the  energy  deposited  per  unit  of 
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penetration  (x)  into  the  material  by  a focused  or  defocused 
electron  beam  incident  on  a flat,  thick  target.  A typical  plot 
of  dE/dx  versus  x is  shown  in  Figure  B-l.  The  linear  portion  of 
the  dE/dx  curve  is  characteristic  of  electron  penetration,  and  its 
extrapolation  defines  the  so  called  Grun  range,  R„(E).  This  is  not 
the  definition  of  extrapolated  range  usually  applied  in  electron 
beam  studies. 

Typically  measurements  are  made  on  thin  foils,  and  the  extra- 
polated range  is  defined  as  the  extrapolated  thickness  that  reduces 
the  most  probably  transmitted  electron  energy  to  zero  (Thomson- 
Whiddington  range  R^) , the  mean  energy  to  zero  (R^),  or  the  extra- 
polated beam  energy,  Ec,  which  reduces  the  number  of  transmitted 
electrons  to  zero  R^Eq).  Ranges  based  on  energy,  such  as  Rg,  are 
more  relevant  for  devices  that  are  sensitive  to  energy,  but  because 
many  published  results  are  in  terms  of  number  ranges,  they  will  also 
need  to  be  referred  to. 

In  Figure  B-2,  ranges  are  shown  as  measured  by  several  inves- 
tigators using  various  methods.  Because  the  definition  of  "range" 
is  not  the  same  in  all  cases,  exact  numerical  agreement  is  not  to  be 
expected.  However,  all  ranges  scale  approximately  as  E^  This 

is  confirmation  of  the  integral,  S,  of  the  Bethe  energy  loss  relation, 
because  those  electrons  that  penetrate  farthest  into  the  target  suf- 
fer few  large  angle  scatterings,  so  S = range.  Note  also  that  the 
ranges  are  essentially  independent  of  the  atomic  number,  because  the 
Bethe  energy  loss  mechanism  depends  primarily  on  the  electron  density 
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Figure  B-l  General  Shape  of  dE/dx  Curves  and  Definition 
of  Grun  Range 
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Figure  B-2  Electron  Ranges  in  Solids  and  Gases 


of  the  substance  and  is  independent  of  Z (Z/A  ~ 2).  Feldman  (Ref. 

73)  did  find  a strong  variation  in  oR  with  Z,  but  his  measurement 
technique  was  shown  to  be  incorrect  by  Kanter  and  Sternglass  (Ref. 

74) . 

Measurements  of  dE/dx  (x)  are  less  common.  The  earliest  are 
those  of  Grun  Range  (Ref.  75),  who  measured  the  fluorescent  radia- 
tion from  nitrogen  gas  at  pressures  of  approximately  1/2  atm  in 
the  energy  range  of  5 to  50  keV.  Cohn  and  Caledonia  (Ref.  76) 
repeated  these  measurements  between  2 and  5 keV,  with  o approximately 
2 torr,  and  found  very  similar  results,  even  though  the  dominant 
fluorescent  energy  mechanism  is  different  at  the  reduced  pressure. 
Everhart  and  Hoff  (Ref.  77  ) measured  dE/dx  by  an  induced  conductiv- 
ity method  in  metal-SiC^-Si  capacitors  between  5 and  25  keV.  Possin 
and  Norton  (Ref.  78)  measured  dE/dx  using  enhanced  etching  in  poly- 
methyl methacrylate  (PMM) , between  5 and  10  keV.  All  four  of  these 
measurements  on  low-Z  solids  found  essentially  the  same  universal 
function  for  dE/dx  when  x was  normalized  in  the  appropriate  fashion 
for  changes  in  beam  energy  and  mass  density  of  the  target  material. 
The  data  of  Possin  and  Norton,  and  Cohn  and  Caledonia,  also  contain 
information  about  the  lateral  spread  of  the  ionization.  Figure  B-3 

shows  this  universal  function,  A (z/R  ),  as  measured  by  Everhart  and 

G 

Hoff  and  Grun. 

Function  dE/dx  is  given  by: 

(B-5) 
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Figure  B-3  Empirical  Universal  Function  A (Z/R) 
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where  Grun  range  = 3Eq  and  = beam  energy  in  keV.  ( B-6) 

Constant  £,  as  determined  by  the  various  investigators,  varies  bet- 

-6  2 

ween  4.0  and  4.57  x 10  g/cm  as  summarized  in  Table  B-l.  The 
materials  studied  were  ^ gas,  PMMM,  and  Si02*  No  dependence  on 
atomic  number  was  found,  as  indicated  by  the  form  of  Equations 
B -5  and  B-6;  - dE/dx  scales  as  the  mass  density. 

Schumacher  and  Mitra  (Ref.  79)  have  used  the  universal  Grun 
fromulation  for  dE/dx,  to  measure  the  thickness  and  composition  of 
thin  films  with  good  success,  and  it  has  been  used  at  the  General 
Electric  Research  and  Development  Center  to  predict  BEAMOS  target 
gain,  again  with  good  success. 

Lateral  Scattering  Measurements 

One  would  like  to  measure  the  ionization  function,  G(r,z), 

for  a focused  beam  entering  a material  normally  at  r = z = 0. 

G is  defined  as  the  ionization  density  per  unit  volume  at  (r,z). 

Experimental  measurements  of  these  profiles  are  rare.  Grun 

shows  one  set  of  profiles  for  nitrogen  gas  at  32  keV.  Ehrenberg 

and  King  (Ref.  80)  show  several  intensity  profiles  that  are  integrals 

of  G along  the  viewing  direction.  These  profiles  are,  unfortunately, 

only  sketches. 

% 

Cohn  and  Caledonia  measured  G by  electron  beam  induced  fluor- 
escence in  nitrogen.  This  experiment  is  shown  schematically  in 
Figure  B-4. 

A focused  electron  beam  enters  nitrogren  gas  at  a low  pressure. 
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approximately  2 torr.  The  size  of  the  luminescent  cloud  is  approx- 
imately 20  centimenters.  Using  a photomultiplier,  Cohn  and  Caledonia 
scanned  the  cloud  in  y and  z,  obtaining  the  sidewise  intensity  pro- 


files, I(y,z): 


1 0,  ■*•)  - f <s 

J . V J- 


-I*-/')*' 


( B-7) 


By  Abel  inversion,  I can  be  inverted  to  give  G(r,z): 


G(S>*)  - ft  J J Cyx) 

The  crosswise  intensity  profiles,  T(z): 

lC*)  uC«> 

XL*)  --ZV  ft*)  rJr-2  frty,z)Jr 


(B-8) 


(B-9) 


are  proportional  to  the  dE/dz  depth  profiles  measured  by  Grun, 
Everhard  and  Hoff,  and  Possin  and  Norton,  and  they  agree  well  with 
these  other  data. 

The  Gn  (r,z)  profiles  that  they  found  are  reproduced  in  Figure 
B -5.  They  found  that  they  obtained  the  same  profiles  (for  electron 
beam  energies  between  2 and  5 keV)  if  all  distances  were  scaled  to 
the  electron  range,  as  defined  by  Grun.  Unfortunately  some  prob- 
lems exist  in  their  analysis.  The  cutoff  in  their  inversion  inte- 
grals, L(z),  was  apparently  determined  from  the  scope  traces  where 
the  intensity,  I(y, z),  fell  to  zero.  However,  an  oscilloscope  trace 
only  resolves  zero  to  about  2 percent  of  the  full  scale, unless  scale 
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changes  were  used.  Cohn  and  Caledonia  claim  to  have  determined 
variations  of  G of  at  least  x 200.  Hence  it  would  be  necessary 
to  measure  changes  in  I(y,z)  of  at  least  x2000,  because  of  the  2irr 
factor  in  the  integrals. 

To  see  this  problem  more  clearly,  consider  calculating  the  Grub 
T(z)  a X(z/R)  from  their  data.  To  do  this  graphically,  one  first 
plots  F = G^(r/R, z)r/R  and  then  integrates: 

« ( B -10) 

Plots  of  GN(r/R,z)r/R  for  several  values  of  z are  shown  in  Figure 
B -6.  The  curves  are  cut  off  at  the  last  data  point  given  by  Cohn 
and  Caledonia.  Clearly,  from  Figure  B-6,  the  data  given  in  the 
paper  do  not  extend  to  large  enough  values  of  r/L  to  perform  the 
normalization  integral  or  the  inversion  integral.  However,  some 
of  the  reduced  data  in  the  paper  (particularly  their  Figures  B -7 
and  8)  indicate  that  their  data  extend  further  than  indicated  by 
Figures  B -5  and  B-6,  so  their  data  may  be  more  valid  than  indicated 
by  this  critique.  We  will  return  to  this  later. 

Possin  and  Norton  (Ref.  78)  of  General  Electric  Corporate  Res- 
earch and  Development  measured  the  lateral  scattering  of  5-  and  10- 
keV  electron  beams  by  the  enhanced  etching  method  in  PMM.  This 
approach  was  originally  proposed  by  Hatzakis;  the  method  is  described 
in  Reference  78.  The  results  are  shown  in  Figures  B-7  through  B-9 
for  the  line  ionization  distribution  I(x,z)  at  5 and  10  keV.  Constant 
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Constant  C «*  1.49  x 10  eV/cm  . Selected  results  are  compared  to 
Cohn  and  Caledonia  in  Figures  B-10  and  B-ll.  These  data  are 
normalized  as  follows:  all  dimensions  are  normalized  to  Grun  range 
RG  = 4.1  x Eq1  75  x 10*6g/cm2  and  IN(x/RG,  z/Rg)  = EQ/RG2I(x,  z) . 

This  is  analogous  to  the  normalization  of  the  Grun  depth  dose  func- 
tion, dE/dx.  Note  that  the  PMM  data  corroborates  the  Cohn  and  Cale- 
donia data  very  well  over  the  range  given  in  their  paper.  The  devia- 
tions for  small  x are  probably  due  to  a failure  of  the  FMMM  method. 
The  deviations  at  large  distances  are  probably  due  to  the  data 
range  problem  of  Cohn  and  Caledonia.  The  poorer  agreement  at  5 keV 
is  probably  a partial  failure  of  the  PMM  method,  again  because  of 
the  small  distances  involved. 

In  general  it  appears  that  the  Cohn  and  Caledonia  data  fairly 
well  describe  the  lateral  scattering  of  bw  (2  to  10  keV)  electron 
beams  in  low  Z solids.  These  data  and  the  Possin  and  Norton  data 
can  be  used  to  describe  the  generation  of  electron  hole  pairs  in  the 
archival  target. 

In  many  calculations  it  is  important  only  that  most  of  the  ioniz- 
ation energy  is  contained  within  a hemisphere  of  depth  0.9  RG  and 
diameter  0.9  R^,.  As  in  the  section  on  resolution  of  the  surface 

U 

diode  target,  I(x,z)  or  G(r,z)  can  be  integrated  over  an  assumed 
electron  beam  spot  profile  to  yield  the  ionization  distribution  in 
the  target  at  various  depths,  Z. 


-323- 


IN  (ELECTRONS) 


T 


Figure  B-10  Selected  Results  Compared  to  Cohn  and 
Caledonia  (5  keV) 
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Figure  B-ll  Selected  Results  Compared  to  Cohn  and 
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